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Abstract. Recent progress has been made in the understanding of the physical 
properties of chromatin - the dense complex of DNA and histone proteins that occupies 
the nuclei of plant and animal cells. Here I will focus on the two lowest levels of 
the hierarchy of DNA folding into the chromatin complex: (i) the nucleosome, the 
chromatin repeating unit consisting of a globular aggregate of eight histone proteins 
with the DNA wrapped around: its overcharging, the DNA unwrapping transition, the 
"sliding" of the octamer along the DNA. {ii) The 30nm chromatin fiber, the necklace- 
like structure of nucleosomes connected via linker DNA: its geometry, its mechanical 
properties under stretching and its response to changing ionic conditions. I will stress 
that chromatin combines two seemingly contradictory features: (1) high compaction 
of DNA within the nuclear envelope and at the same time (2) accessibility to genes, 
promoter regions and gene regulatory sequences. 
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1. Introduction 

Higher developed organisms face the problem to store and retrieve a huge amount of 
genetic information - and this in each cell separately. For instance, the human genom 
corresponds to 3 billion base pairs (bp) of the DNA double helix, two copies of which 
make up two meters of DNA chains that have to be stored within the tiny micron-sized 
nucleus of each cell pp. These two meters are composed of 46 shorter DNA pieces, 
each of which, if not condensed, would form a swollen coil of roughly lOO/im diameter 
- clearly much too large to fit into the nucleus. On the other hand, the densely packed 
genom would form a ball of just ~ 2 /im diameter due to the huge aspect ratio of contour 
length - 2 meters - versus diameter - 20 A - of the DNA material. Hence, the DNA 
indeed fits into the nucleus but a suitable compaction mechanism is required. 

DNA is a highly charged macromolecule carrying two negative elementary charges 
per 3.4 A. In the presence of multivalent counterions DNA condenses into dense - often 
toroidal - aggregates |2] , resembling the DNA packaged in virial capsids |S1 111 Ej • In 
viruses the DNA has just to be stored whereas such a simple way of DNA compaction 
cannot work for the long DNA chains in eucaryotic cells (cells of fungi, plants and 
animals) where many portions of the DNA have to be accessible to a large number 
of proteins (gene regulatory proteins, transcription factors, RNA polymerases etc.). 
Therefore the substrate that these proteins interact with must be much more versatile to 
allow access to certain regions of the DNA and to hide (i.e. to silence) other parts. That 
way each cell can regulate the expression of its genes separately according to its state 
in the cell cycle, the amount of nutrients present etc. Furthermore, the differentiation 
of cells into the various types that make up a multicellular organism relies to a large 
extent on the way the DNA - which is identical in all the cells - is packaged. 

The substrate that combines all these features is chromatin, a complex of DNA 
and so-called histone proteins. In 1974 it has been realized that the fundamental unit 
of chromatin is the nudeosome [HI E] : roughly 200 bp of DNA are associated with one 
globular octameric aggregate of eight histone proteins consisting of two molecules each 
of the four core histones H2A, H2B, H3 and H4. A stretch of 147 bp DNA is wrapped 
in a 1- and- 3/4 left-handed superhelical turn around the octamer and is connected via 
a stretch of linker DNA to the next such protein spool. Each octamer together with the 
wrapped DNA forms a nudeosome core particle with a radius of ~ 5 nm and a height 
of ~ 6 nm which carries a large negative electric charge |H1 E] • 

While the structure of individual core particles is now documented in great detail 
mainly on the basis of high-resolution X-ray analyses ^Ql 1^ , much less is known about 
the higher-order structures to which they give rise. When the fiber is swollen - as this 
is the case for low ionic strength - it has the appearance of "beads-on-a-string." It is 
sometimes referred to as the " 10-nm fiber" since its "beads" have ~ 10 nm diameter |12j . 
With increasing salt concentration, heading towards physiological conditions (roughly 
100 mM), the fiber becomes denser and thicker, attaining a diameter of ~ 30 nm 
[T^ . Longstanding controversy surrounds the structure of this so-called 30-nm fiber 




''lOnm-fiba'' 



Figure 1. Steps of the DNA compaction into chromatin. The DNA molecule of 
length ~ 1 cm is compacted with the help of 10^ histone octamers leading to 10000- 
fold reduction of its original length (see text for details) 

|14[ IT^ [TBI EI- In the solenoid models ITHl IT^ it is assumed that the chain of 
nucleosomes forms a hehcal structure with the hnker DNA bent in between whereas 
the zig-zag or crossed linker models pUl 1^ 122 12^1 121] postulate straight hnkers that 
connect nucleosomes that are located on opposite sites of the fiber. The higher-order 
folding of the 30-nm fiber into structures on scales up to microns is yet to be elucidated. 
In Fig. 1 I sketch the steps of the DNA folding starting with a DNA chain of length 
~ 1 cm and ~ 10^ octamers and ending up at the highly condensed chromosome. This 
highly condensed structure occurs before cell division and contains the chain and its 
copy neatly packaged for the distribution into the two daughter cells. The size of the 
chromosome is ~ 10000 times smaller than the contour length of the original chain. 

It is instructive to draw a comparison between the structure and function of 
chromatin to that of a daily-life example: the library. As the nucleus stores the long one- 
dimensional string of basepairs, so does the library contain the huge one-dimensional 
string of letters, the text written down in all of its books. A book like Ref. jl] contains 
~ 10 km text, a library with 10000 books stores roughly 100000 km text! How can the 
user find and retrieve the little piece of information of interest? The way this is handled 
is that the text is folded in a hierarchical fashion in lines, pages, books and shelves. This 
makes it is relatively easy - with the help of a few markers - to find the corresponding 
text passage. Furthermore, all the text is stored in a dense fashion but the book of 
interest can be taken out of the shelve and opened at the appropriate page without 
perturbing the rest of the library. Apparently, the result of this hierarchical structure 
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is a relative high efficiency in storing a huge amount of information in a relatively small 
space and, at the same time, having a high accessibility to it. 

The similarities between hierarchies in the library and in chromatin are pretty 
obvious. What is less obvious and in many respects still an open question is how 
the dense chromatin structure can be opened locally to allow access to its genes. 
As mentioned above for one nucleosomal repeat length, typically 200 bp, 147 bp are 
wrapped around the octamer, i.e., roughly 75 percent of the DNA chain is tightly 
associated to the histone aggregates. It is known that many essential proteins that 
interact with DNA do not have access to DNA when it is wrapped (reviewed in Ref. |25j). 
Moreover, even the unwrapped sections - the linker DNA - are somewhat buried inside 
the dense 30-nm fiber. Therefore, it is necessary for the cell to have mechanisms at hand 
to open - unfold - the fiber and then, somehow, to unwrap the DNA from the protein 
spools or to temporarily remove them from the DNA piece of interest. 

This leads to the problem of how the chromatin structure changes its shape with 
time. As this structure involves length scales of many orders in magnitude (from 
Angstroms to microns) so do its dynamical processes take place in a wide range of 
time scales, beginning with fiuctuations of the nucleosome structure in the micro- to 
millisecond timescale [2^1 up to large scale variations in the condensation degree of the 
chromosome that follows the cell cycle with a typical period of hours to days pp. 

It has been shown through competitive protein binding to nucleosomal DNA p7|l^ 
that thermal fiuctuations might lead to partial unwrapping of the DNA from the 
nucleosome. This mechanism provides intermittent access to nucleosomal DNA. Not 
only unwrapping but even "sliding'^ of nucleosomes along DNA seems to be facihtated 
by thermal fiuctuations as it has been demonstrated in well-defined in vitro experiments 
j^niEOlEIl- Whereas this kind of repositioning dynamics is quite slow - even at elevated 
temperatures the timescale is of the order of minutes to hours - nucleosome repositioning 
appears also to be of great importance in vivo where it is aided by chromatin-remodeling 
complexes, large multi-protein complexes that use energy by burning ATP (reviewed in 
Refs. |Zll321ISnilSlllSSl)- These complexes might catalyze and direct the displacement of 
nucleosomes out of regions where direct access to DNA has to be granted (like promoter 
regions of transcriptionally active genes). 

Based on a range of experiments it has also been speculated that RNA polymerase, 
the protein complex that transcribes - copies - genes, can itself act "through" 
nucleosomes without having to disrupt the structure completely jHEl IHZl IHHl IBHl 1^ - 
An appealing picture is the idea that the polymerase gets around the nucleosome in a 
loop. Such a mechanism is especially important since genes are typically distributed 
over a DNA portion of the length of lOO's to lOOOO's of bp which means that there are 
a few up to lOOO's of nucleosomes with which a polymerase has somehow to deal with 
during transcription. 

Concerning the unfolding of the 30-nm fiber, I mentioned already above that a 
decrease in the ionic strength leads to a swelling of the fiber. Obviously, in vivo such a 
global swelling is not possible within the tiny space available within the nucleus. Hence 
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only the unfolding of local regions should be expected as the above given analogy of a 
library suggests. Experimental observations seem to indicate that the swelling degree 
can indeed be tuned locally by the acetylation of the lysine-rich (i.e. cationic) tails 
of the eight core histones that appear to be long flexible polyelectrolyte chains [H] 
that extend out of the globular part. Furthermore, transcriptionally active chromatin 
portions show a depletion in the linker histone HI, a cationic proteins that is believed 
to act close to the entry-exit region of the DNA at the nucleosome. As long as HI is 
present the fiber is relatively dense and the individual nucleosomes are inert against 
thermal fluctuations (no unwrapping, sliding or transcription through the octamer). If 
HI is missing the chromatin fiber appears to be much more open and the nucleosomes 
become "transparent" and mobile due to thermal fluctuations |31j . 

Undoubtedly chromatin lies in the heart of many essential biological processes 
ranging from gene expression to cell division. Most of these processes are controlled 
by a huge amount of specific proteins. Their investigation clearly belongs to the 
realm of biologists and is beyond the scope of a physicist. Nevertheless, many insights 
gained in this field were achieved through in vitro experiments under relative well- 
controlled conditions, in some cases essentially only involving DNA and histone proteins. 
Furthermore, new physical methods like micromanipulation experiments allowed to gain 
access to certain physical properties of chromatin. The purpose of this paper is to discuss 
some of these results and, especially, to review the physical theories they gave rise to. 

The interest of physicists into single nucleosomes was mainly sparked by the above 
mentioned fact that the core particle carries a large negative net charge. This is due 
to the fact that much more negatively charged DNA is wrapped around the cationic 
octamer than necessary for its neutralization. Beginning around 1998 a considerable 
activity started among several research groups to explain on which physical facts 
overcharging is based (recently reviewed in [42j), a phenomenon not only occurring 
in chromatin but also in DNA-lipid complexes, multilayer adsorption of polyelectrolytes 
etc. To gain insight into this the nucleosome was translated into different types of toy 
models, usually consisting of one charged chain and an oppositely charged sphere and 
the amount of wrapping was calculated |13lllligiinill3llHlliniEDlEIlEaE3ElESl 

(cf. also the related studies [HEiisziinEinniinniiniiinaEHiiEiiMiEiiisziE^ 

It turned out that overcharging is a fairly robust phenomenon that occurs even if certain 
mechanisms are neglected (like counterion release upon adsorption). It was also shown 
that multi-sphere-complexation can lead to undercharged systems jlHl EI] • 

The relevance of these toy models for "real" nucleosomes might be questionable, 
especially since the binding sites of the DNA to the octamer are relatively specific jlO|ITT] 
and since under physiological conditions (100 mM) the screening length is fairly short 
- 10 nm, i.e., half the diameter of the DNA double helix. Nevertheless, these models 
might give some insight into the unwrapping transition that takes place when the DNA 
which is fairly rigid on the nucleosome length scale unwraps from the nucleosome due 
to a decrease of the adsorption energy. This can be achieved by a change in the salt 
concentration and has indeed been observed experimentally [SI [121 • ^ simple approach to 
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this problem comparing adsorption energies between a "sticky" spool and a semiflexible 
chain and its bending energy has been given in Ref. [73] and led to the prediction of an 
unwrapping transition in an all-or-nothing fashion, i.e., the cylinder is fully wrapped by 
the chain or not wrapped at all. A more complex picture has been found in Refs. [^ ITF] 
where the unwrapping has been calculated for complexes with short chains (" nucleosome 
core particles"). Here the chain showed different degrees of wrapping (dependent on the 
ionic strength) which seems to agree fairly well with the corresponding experimental 
observations on core particles jHl IZ2] • It was shown in Refs. [SSI IZ1| that besides the 
wrapped and unwrapped structures there is for longer chains the realm of open multi-loop 
complexes (so-called rosettes) that have now also been observed in computer simulations 

The nucleosome repositioning, mentioned above, is another single-nucleosome 
problem that has been investigated theoretically jZSl HEl (ZZl • It is suggested that the 
nucleosome " sliding" is based on the formation of loops at the ends of the nucleosome 
which then diffuse as defects around the spool, leading to the repositioning. This process 
somewhat resembles polymer reptation in the confining tube of the surrounding medium 
(a polymer network or melt), cf. Ref. [ZHl (ZHl ISOl • A different "channel" for repositioning 
might be a cork screw motion of the DNA helix around the octamer induced by twist 
defects which has now also been studied theoretically [77j . 

On the level of the 30-nm fiber recent progress has been made in the visualization 
of these fibers via electron cyromicroscopy [23 • The micrographs reveal for lower ionic 
strength structures that resemble the crossed linker model mentioned above. However, 
for increasing ionic strength the fibers become so dense that their structure still remains 
obscure. An alternative approach was achieved via the micromanipulation of single 
30-nm fibers [HH IH21 jEH] • The stretching of the fibers showed interesting mechanical 
properties, namely a very low stretching modulus for small tension, a force plateau 
around 5 pN and, at much higher tensions, saw-tooth-type patterns. 

These experimental results led to a revival of interest in 30-nm fiber models. 
Theoretical studies [23 IHH IHS| as well as computer simulations [SSI IHZ| attempted to 
explain the mechanical properties of the fiber. What all these models have in common 
is that they assume straight linkers in accordance with the experimental observations, 
at least found for lower ionic strength [2ni- The low stretching modulus of the fibers 
is then attributed to the bending and twisting of the linker DNA which is induced by 
the externally applied tension. Indeed, all the models find a relatively good agreement 
with the experimental data. The bpN force plateau is usually interpreted as a reversible 
condensation-decondensation transition of the fiber which can be attributed to a small 
attractive interaction between the nucleosomes [211IHS1- This is also in good agreement 
with recent studies on single nucleosome core particles where such a weak attraction has 
been observed and attributed to a tail-bridging effect [HHllSni- Finally, the sawtooth- 
pattern which leads to a non-reversible lenghtening of the fiber probably reflects the 
"evaporation" of histones from the DNA I9()j. 

The geometry of the crossed-linker models suggest that the density of nucleosomes 
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in the fiber depends to a large extend on the entry-exit angle of the DNA at the 
nucleosomes. It is known that in the presence of the linker histone the entering 
and exiting strands are forced together in a "stem" region j^. A recent study pT] 
investigated how the electrostatic repulsion between the two strands dictates this entry- 
exit angle. Especially, it was shown how this angle can be controlled in vitro via a change 
in the salt concentration. It also has been speculated that via biochemical mechanisms 
that control the charges in the entry-exit region (like the acetylation of certain histone 
tails) the cell can locally induce a swelling of the fiber [T7j . 

Before discussing in the following all the above mentioned issues concerning the 
physics of chromatin, I note that there are also important studies using a bottom-down 
approach by studying the physical properties of whole chromosomes that have been 
extracted from nuclei preparing for cell division. Via micropipette manipulation of 
these mitotic chromosomes it was demonstrated that they are extremely deformable by 
an externally applied tension 123 IMj and that a change in ionic strength induces a 
hypercondensation or decondensation of the chromosome, respectively jOH]- Meiotic and 
mitotic chromosomes were compared to simple polymer systems like brushes and gels 
|96j . A problem with this bottom-down approach is still the lack of knowledge of the 
chromosome structure at this level and of the proteins that cause them. I will therefore 
dispense with giving a discussion on this subject. 

This review is organized as follows. In the next section I give a discussion of 
single nucleosome problems. After providing some experimental facts on the structure 
of the nucleosome (Section 2.1), I discuss simple model systems (Section 2.2), the 
unwrapping transition (Section 2.3) and the nucleosome repositioning (Section 2.4). 
Section 3 features the next level of folding, the 30-nm fiber. I briefly review some of the 
proposed models (Section 3.1) and then give a systematic account on the crossed-linker 
model (Sections 3.2 and 3.3) and its mechanical response to stretching, bending and 
twisting (Section 3.4). Then fiber swelling (Section 3.5) is discussed. Section 4 gives a 
conclusion and outlook. 

2. Single nucleosome 

2.1. Experimental facts on the core particle 

The structure of the nucleosome core particle is known in exquisite detail from X-ray 
crystallography: the octamer in absence of the DNA was resolved at 3.1 A resolution 
pT7] and the crystal structure of the complete core particle at 2.8 A resolution [10] and 
recently at 1.9 A ^I]. The octamer is composed of two molecules each of the four 
core histone proteins H2A, H2B, H3 and H4. Each histone protein has in common 
a similar central domain composed of three a-helices connected by two loops. These 
central domains associate pairwise in form of a characteristic "handshake" motif which 
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Figure 2. Schematic views of the nucleosome core particle. At the top picture the 
upper half of the 8 core histories and the nucleosomal DNA are depicted. At the bottom 
a simplified model is displayed where the octamer is replaced by a cylinder and the 
DNA by a wormlike chain. Also indicated are the dyad axis and the DNA superhelix 
axis. 

leads to crescent-shaped heterodimers, namely H3 with H4 and H2A with H2B|. These 
four " histone-fold" dimers are put together in such a way that they form a cylinder with 
~ 65 A diameter and ~60 A height. With grooves, ridges and relatively specific binding 
sites they define the wrapping path of the DNA, a left-handed helical ramp of 1 and 
3/4 turns, 147 bp length and a ~ 28 A pitch. In fact, the dimers themselves are placed 
in the octamer in such a way that they follow this solenoid: They form a tetrapartite, 
left-handed superhelix, a spiral of the four heterodimers (H2A-H2B)^, (H3-H4)^, (H3- 
H4)^ and (H2A-H2B)^. This aggregate has a two-fold axis of symmetry (the dyad axis) 
that goes through the (H3-H4)2-tetramer apex and is perpendicular to the superhelix 
axis. A schematic view of the nucleosome core particle is given in Fig. 2. 

There are fourteen regions where the wrapped DNA contacts the octamer surface, 
documented in great detail in Refs. jTHI and ^l]. These regions are located where the 
minor grooves of the right-handed DNA double helix face inwards towards the surface of 
the octamer. Each crescent-shaped heterodimer has three contact points, two at its tips 
and one in the middle, altogether making up 12 of the 14 contacts. Furthermore, at each 
end of the wrapped section (the termini of the superhelix) there is a helical extension 

X At physiological conditions stable oligomeric aggregates of the core histones are the H3-H4 tetramer 
(an aggregate of two H3 and two H4 proteins) and the H2A-H2B dimer. The octamer is stable if it is 
associated with DNA or at higher ionic strengths. 
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of the nearby H3 histone making contact to a minor groove. At each contact region 
there are several direct hydrogen bonds between histone proteins and the DNA sugar- 
phosphate backbone ^0] as well as bridging water molecules ^I]. Furthermore, there is 
also always a (cationic) arginine side chain extending into the DNA minor groove. The 
free energies of binding at each sticking point is different which can be concluded from 
the fact that for each binding site there is a different number of hydrogen bonds located 
at different positions; this is also reflected in the fluctuations of the DNA phosphate 
groups in the nucleosome crystal [Ej. However, a reliable quantitative estimate of these 
energies is still missing. 

An indirect method to estimate the adsorption energies at the sticking points is 
based on studies of competitive protein binding to nucleosomal DNA 12HI HH] • Many 
proteins are not able to bind to DNA when it is wrapped on the histone spool due to 
steric hindrance from the octamer surface. However, thermal fluctuations temporarily 
expose portions of the nucleosomal DNA via the unwrapping from either end of the 
superhelix. It was demonstrated that sites which are cut by certain restriction enzymes 
showed - compared to naked DNA - an increased resistance to digestion by these 
enzymes when they are associated to the octamer. Furthermore, the farer these sites 
are apart from the termini of the superhelix the less frequent they become exposed for 
cutting. The rate for digestion is reduced roughly by a factor 10~^ for sites at the 
superhelical termini and ~ 10~^ — 10~^ for sites close to the center of the nucleosomal 
DNA portion (i.e., close to the dyad axis). From these findings one can estimate that 
the adsorption energy per sticking point is of order ~ 1.5 — . 

It is important to note that the 1.5 — do not represent the pure adsorption 

energy but instead the net gain in energy which is left after the DNA has been bent 
around the octamer to make contact to the sticking point. A rough estimate of this 
deformation energy can be given by describing the DNA as a semiflexible chain with 
a persistence length Ip of ~500 A [99j. Then the elastic energy |l()()j required to bend 
the 127 bp of DNA around the octamer (10 bp at each terminus are essentially straight 
[10]) is given by 

Eelastic/kBT = lpl/2R^ (l) 

Here / is the bent part of the wrapped DNA, ~ 127 x 3.4 A= 431 A and Rq is the 
radius of curvature of the centerline of the wrapped DNA (cf. Fig. 2) which is roughly 
43 A JU]. Hence, the bending energy is of order 58/cb^- This number, however, has to 
be taken with caution. First of all it is not clear if Eq. Q is still a good approximation 
for such strong curvatures. Then it is known that the DNA does not bend uniformly 
around the octamer [lOj and finally the DNA might show modified elastic properties 
due to its contacts with the octamer. Nevertheless, using this number one is led to the 
conclusion that the bending energy per ten basepairs, i.e., per sticking site, is of order 
mkBT/lA ~ AkpT. 

Together with the observation that the net gain per sticking point is ~ 2A;bT this 
means that the pure adsorption energy is on average ~ GksT per binding site. Note that 
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the huge pure adsorption energy of ~ Ofc^T x 14 ~ SbksT per nucleosome is cancelled 
to a large extend by the ~ 58k bT from the DNA bending, a fact that has important 
consequences for the unwrapping transition discussed in Section 2.3 and especially, for 
the nucleosome repositioning reviewed in Section 2.4. 

Of great importance are also flexible, irregular tail regions of the core histones which 
make up ~28% of their sequences jH]. Each histone has a highly positively charged, 
flexible tail (which is the N-end of the polypeptide chain IJ) that extends out of the 
nucleosome structure. Some of them exit between the two turns of the wrapped DNA, 
others on the top or bottom of the octameric cylinder. These N-tails are extremely 
basic due to a high amount of lysine and arginine aminoacids (aa's). They are sites 
of posttranslational modification and are crucial for chromatin regulation. Especially, 
the tails have a strong influence on the structure of the 30nm chromatin fiber as I will 
discuss in more detail in Section 3. X-ray scattering data on core particles [HEj suggest 
that the tails are adsorbed on the complex for small ionic strengths and extended at 
high salt concentrations (cf. Fig. 8 in that paper). If the octamer is associated with a 
longer DNA piece the N-tails desorb at higher ionic strength. 

Finally, let me mention the amount of charges found on the nucleosome core particle 
[S]. The histone octamer contains 220 basic side chains (arginine and lysine). From 
these are about 103 located in the flexible histone tails mentioned above. The rest, 
117 residues, are in the globular part of the octamer, of which 31 are exposed to the 
solvent, the rest being involved in intra- and interprotein ionic interactions. On the 
other hand, one has 147 bp of DNA wrapped around the octamer, each contributing 
two phosphate groups. Hence there are 294 negative charges from the DNA versus 220 
positive charges of the octamer (or even less, 134 - not counting the charges buried 
inside the octamer), i.e., the nucleosomal complex is overcharged by the DNA. At first 
sight this is a surprising fact and indeed led to the development of simplified toy models 
containing charged spheres and oppositely charged chains which I will discuss in the 
following section. 

2.2. Poly electrolyte- charged sphere complexes as model systems for the nucleosome 

The complexation of (semi-)flexible polyelectrolytes and oppositely charged macroions 
is an important ingredient in biological processes. For instance, the non-specific part of 
the interaction between DNA and proteins is governed by electrostatics |l()lj . In fact, 
the nucleosome is an example of this kind of interaction. A number of experimental 
[ma CUSl imi CUS] and theoretical studies UHl CUni EUl CUHl have demonstrated 
that the complexation of highly charged macroions (e.g. DNA) is governed by an unusual 
electrostatics mechanism: counterion release. The free energy of complexation is then 
dominated by the entropy increase of the released counterions that had been condensed 
before complexation. This electrostatic contribution to the free energy has to compete 
with the energy cost of deforming one or both macromolecules to bring them in close 
contact. 
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In the next section I discuss how the counterion release leads to an overcharged 
sphere-chain complex which might be considered as a simplified model system of the 
nucleosome. I also consider the case of a polyelectrolyte chain placed in a solution of 
oppositely charged spheres (Section 2.2.2). Both sections follow closely the treatments 
given in Refs. ■ Afterwards I review studies where it was found that overcharging 

occurs also in weakly charged systems due to "standard" electrostatics (Section 2.2.3). 
Section 2.2.4 is devoted to physiological conditions (strong screening) and also to the 
question whether such toy models can "explain" the large net charge of nucleosomes. 

2.2.1. Single-sphere complex (highly charged case) Consider a single sphere of radius 
Rq with its charge eZ homogeneously smeared out over the surface (the "octamer") 
and a fiexible rod with a charge per unit length — e/6, persistence length /p, contour 
length L S> -Rq and radius r (the "DNA chain"). Both are placed in a salt solution 
characterized by a Bjerrum length Ib = /eksT (e: dielectric constant of the solvent; 
in water e = 80 and Ib = 7 A at room temperature) and a Debye screening length 

= (SttCs/b)^^^^; furthermore the solvent is treated as a continuum. Clearly, such a 
model neglects all the intricate features of the nucleosome discussed in Section 2.1. It 
is, however, indispensable to start from such a simple model system to identify general 
features that occur in polyelectrolyte-macroion complexes. 

I will focus in this section on salt concentrations Cg that are sufficiently small such 
that is large compared to the sphere radius, kRq <^ 1. The persistence length 
is assumed large compared with Rq. The chain is here highly charged which means 
that the so-called Manning parameter ^ = Ib/^ is required to be much larger than one. 
In this case (1 — ^~^) L/h ~ L/h counterions are condensed on the chain reducing the 
effective line charge density to the value —c/Ib The entropic cost to "confine" 

a counterion close to the chain is QkBT with Q = 2ln (4^/t~^/r) |KT |lll2j . This leads to 
the following entropic electrostatic charging free energy of the isolated chain in the salt 
solution: 

Fchain (L) ^ f2L 

ksT ~ b ^ ' 

On the other hand the corresponding electrostatic charging free energy of the spherical 
macroion of charge Z is given by 

Fsphere (■^) ^ J '^Rq' 1"^! ^ ^max /q\ 

ksT ~\ \Z\n (Z) for \Z\ > Z^ax ^ ' 

where VL{Z) = 2\n{\Z\lBK^^/Rl) |llUj and Z^ax ~ ^Rq/^b (see below). For weakly 
charged spheres {\Z\ < Zmax) Fsphere is the usual electrostatic charging energy. In the 
highly charged case, \Z\ ^ Z^axi most of the counterions are localized close to the 
sphere with an entropic cost VL (Z) fc^T per counterion leading to Eq. Olfor |Z| Z^ax- 
Only the small fraction Z^ax/Z of counterions is still free leading to an effective sphere 
charge Z^ax- The value of Z^ax follows from the balance of electrostatic charging energy 
IbZ^^^/2Ro and counterion entropy —Q (Z) Z^ax [HHj- 
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The total free energy of the sphere-chain complex can be determined as follows. 
Assume that a length I of the chain has been wrapped around the sphere. Divide the 
sphere-chain complex in two parts: the sphere with the wrapped part of length / of the 
chain and the remaining chain of length L — l. The first part - which I will refer to as the 
"complex" - carries a net charge Z (/) = Z — l/b. The electrostatic free energy Fcompi (0 
of the complex is then estimated to be equal to Fgphere {Z (/)) (neglecting higher-order 
multipole contributions). There is a special length /j^o = bZ, the isoelectric wrapping 
length, at which Z (liso) = 0. The usual principle of charge neutrality would lead one to 
expect that the total free energy is minimized at this point. 

The total free energy of the sphere-chain complex is approximately given by the 
following terms j51j : 

(^) Fcompi (0 ~^ ^chain (-^ ~^ ^compl— chain (0 ~^ ^elastic (0 (^) 

The first two terms have already been discussed. The third term is the electrostatic free 
energy of the interaction between the complex and the remainder of the chain which is 
of the order 



Fcompi— chain (0 



Z* (/) In {kRq) (5) 



ksT 

where Z* (/) is the effective charge of the complex (the smaller value of Z (/) and Zmax)- 
The final term in Eq. 0] describes the elastic energy of the wrapped portion of the chain 
that has a typical curvature I/-R0 and is given already above, Eq. [T] 

Following Ref. jHI] the two cases \Z < Zmax and \Z > Zmax have to 
be treated separately. The first case applies for wrapping lengths / between Imin = 
kso - hZraax and Imax = kso + bZ^ax- The free energy intakes then the form 



Ib (r, , CI 



, ^ — 1 - , I + -; — I- const. (6) 
keT 2Rq\ h) h ^ ' 

where C = lpb/2Rl — In (kRq) — fl. For the second case, when \Z (/)| > Zmax, one finds 

— + const (7) 

with = lpb/2Rl — Q^Q. The "— " sign refers to the case / < Imin (i-e. Z (/) > Zmax) 
when for every segment b of adsorbed length a negative counterion of the sphere and a 
positive counterion of the chain are released while the " +" sign refers to the case / > Imax 
(equivalently, Z (/) < —Zmax) when for every adsorbed segment a positive counterion is 
transferred from the chain to the sphere leading to a change kp — of its entropy. 
The three different cases are depicted in Fig. 3. 

Using Eqs.|Bland[7|one can describe the complexation as a function of chain stiffness. 
For large Ip, B~ > and there is no wrapping; the free energy is minimized for / = 0. 
There is, however, still the possibility of more open complexes with many point contacts 
between the sphere and the chain which we shall discuss in Section 2.3. As is reduced, 

changes sign which marks the onset of wrapping. For B^ < and C > the 
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Figure 3. Schematic view of the complex between a highly charged chain and an 
oppositely charged sphere. Depicted are three scenarios. From top to bottom: For 
short wrapping lengths complexation is driven by release of counterions from sphere 
and chain (top picture). For intermediate values (around the isoelectric wrapping 
length) all counterions of the sphere have been released, further complexation leads, 
however, still to release of counterions from the chain. For even larger wrapping lengths 
there is no further counterion release. 



minimum of Fi {I) lies in between Imin and Uso- According to Eq. IHlthe position of the 
free energy minimum /* is given by 

/* = kso - CR^/i (8) 

a result first given by Park, Bruinsma and Gelbart jJHl- Further reduction of Ip leads 
to increasing /* until the complex reaches the isoelectric point at C = 0. For smaller Ip, 
C < and according to Eq.|Sl/* > /j^o so that the complex is overcharged. Consequently, 
for a fully flexible chain with /p = 0, the complex is always overcharged. The critical 
persistence length below which complexes are overcharged is /p = 2 (fi + In {kRq)) Rl/b. 

It can be seen clearly from this line of arguments that it is the release of counterions 
from the chain that drives the overcharging. What opposes this effect is the charging 
energy of the complex, the repulsion between the chain and the overcharged complex 
and, most importantly, the bending stiffness of the chain. 

2.2.2. Multi-sphere complex (highly charged case) In the previous section it was 
discussed how the release of counterions from the chain upon adsorption causes 
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overcharging. Even though the focus of this chapter is on single-nucleosome properties 
it is instructive to consider next the case of a highly charged chain placed in a solution 
of oppositely, highly charged spherical macroions. This case has been investigated by 
myself, Bruinsma and Gelbart jM]. The solution is represented as a reservoir with 
concentration Cm of uncomplexed spheres. The chemical potential is the sum of the 
usual ideal solution term and the electrostatic free energy of a spherical macroion with 
charge Z > Z^ax (cf. Eq. E)): 

= In (c^Rt) + Zn^Zn (9) 

The number of spheres that complex with the chain are determined by requiring this 
chemical potential to equal that of the complexed spheres. In Ref. we assumed a 
beads-on-a-string configuration, with a mean spacing D between spheres. The Euclidean 
distance S between the beginning and the end of this configuration is related to the 
number of complexed spheres via N = S/D. The wrapping length I per sphere follows 
from S and L to be L ~ A^/ + — 2NRo. The Gibbs free energy of the bead-on-a-string 
configuration is given by 

G {N, I) = NF, (/) + {N, I) - flsphereN (10) 

where Fi is the above given single-sphere complex free energy, Eq. EJ and Fint is the 
interaction between the complexed spheres. For sphere-sphere spacings D [N, I) = 
S{N,l) /N with 2Rq < D < k,~^, the complexed spheres feel a mutual electrostatic 
repulsion approximately given by (for |2' (/)| < Z^ax) 

keT ~ D{N,l) ^ ^ 

with A ~ 2 In {k^-^ N / L) . For D ^ 2Rq adjacent spheres interact via a strong excluded 
volume interaction, for D ^ the electrostatic interaction is screened. 

G {N, I) has to be minimized with respect to both N and /. We showed in Ref. [HI] 
that due to the large chemical potential of the spheres (last term in Eq. ITTH) it is 
energetically favorable to keep adding spheres to the chain up to the point when 
D ^ 2Rq. At this point the hard-core repulsion terminates complexation and the chain 
is completely "decorated" with spheres. It follows then that the number A^ of spheres 
and the wrapping length I per sphere are related via N c::^ L/l, i.e., essentially the whole 
chain is in the wrapped state. This argument holds for any Z^m < I < Imax- Because 
of this relation (namely N = L/l) the Gibbs free energy depends only the number N of 
complexed spheres: 

^ - {(^ + 1) 2^ ('•■» - ^/^)^ - ^) + (^^' 

Clearly, the first term of Eq. El favors the isoelectric configuration L/N = Uso- 
However, because of the second term, we can lower the free energy further by increasing 
A^ beyond L/liso- This is not a small effect since ^sphere/ ^bT is of order Z ^ Z^ax while 
the first term of Eq. [T21 the capacitive energy, is of order (Ib/Rq) Z^iax ~ Zmax (since 
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Zmax ~ Rq/Ib)- The spheres in the many-sphere complex are thus undercharged. The 
optimal wrapping length is as follows 



Physically, this effect can be illustrated by first setting L/N = Uso- In that case the 
complex is isoelectric. Now add one more sphere. By equally redistributing the chain 
length between the + 1 spheres, one has an individual wrapping length I = L/ [N + 1) 
close to the isoelectric one. Therefore the previously condensed counterions of the added 
sphere are released and increase their entropy. By adding more and more spheres - while 
reducing I = L/N - more and more counterions are liberated. 

In both cases, the single-sphere case of the previous section and the multi-sphere 
case discussed here, it is the counterion release that is the driving force which brings 
oppositely, highly charged macroions together. In the first case the release of the cations 
of the chain is responsible for bringing more monomers to the complex than necessary 
for its neutralization; in the second case it is the release of the anions of the spheres that 
attracts more spheres to the chain than "optimal" and the spheres are undercharged. 
Comparing the similarities between the two cases it might be more appropriate to say 
that in the latter case the spheres overcharge the chain. 

Finally, there is also the possibility to have a solution of chains and spheres in a 
certain stoichiometric ratio such that there are N < L/l* (with /* being the single-sphere 
wrapping length, Eq. |HJ. Then essentially all spheres will complex with the chains (due 
to the large contribution from the counterions to the chemical potential, Eq. IHl). Since in 
this case there is enough chain available, each sphere will be overcharged by the chain. 
Indeed, using Eq. HIl with ^sphere = we found in Ref . [SI] 



as the optimal wrapping length. This is the single chain wrapping length, Eq. |H1 with 
a slightly reduced deviation from the isoelectric point due to the electrostatic repulsion 
between the complexed spheres, Eq. [TT] 

2.2.3. Weakly charged case The first theoretical models on sphere-chain complexation 
were presented in 1999, each of which used quite a different approach to this problem 
^3 E3 ESj • Park, Bruinsma and Gelbart 43ij considered a semifiexible and highly 
charged chain and showed that counterion release leads to overcharging, as discussed in 
Section 2.2.1, cf. Eq. |H1 The other studies jUllinillE] considered weakly charged chains 
and arrived at the conclusion that also in this case overcharging should be a common 
phenomenon. 

Gurovitch and Sens studied a point like central charge (the "sphere") and a 
connected chain of charges (the flexible chain). On the basis of a variational approach 
(self-consistent field theory using an analogy to quantum theory [78|) they came to the 
conclusion that the chain collapses on the central charge even if the total charge of the 




(13) 




(14) 
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resulting complex becomes "overcharged." The critical polymer charge up to which this 
collapse occurs is 15/6 times the central charge. In a subsequent discussion |114j it 
became clear that this number has to be taken with caution and that other effects, like 
the formation of tails, loops etc. were not included in the class of trial function which 
were used in that study. 

Mateescu, Jeppesen and Pincus jH] used a purely geometrical approach in order 
to calculate the zero-temperature state of a complex of a sphere and a perfectly flexible 
chain in the absence of any small ions (no salt, no counterions) . They divided the chain 
into two regions, one straight tail (or two tails on opposite sites of the complex) and 
a spherical shell around the macroion. The only approximation in that study was to 
uniformly smear out the monomer charges within the spherical shell. Starting from a 
point-like sphere and then gradually increasing its radius Rq they found the following 
typical scenario (cf. Fig. 1 in that paper): For very small Rq the complex is slightly 
undercharged and shows two tails. With increasing Rq more and more chain wraps 
around the sphere leading to an overcharging of the complex. For sufficiently large 
sphere radius the whole chain is adsorbed. Before this point is reached there are - for 
sufficiently long chains - two jump-like transitions: one from the two-tail to the one-tail 
configuration and then one from the one-tail case to the completely wrapped state. 

Finally, Netz and Joanny [IHI considered the complexation between a semiflexible 
chain and a sphere. For simplicity, they considered a two-dimensional geometry and 
calculated - using a perturbative approach - the length of the wrapped section and 
the shape of the two-tails for different salt concentrations. That study focuses on 
the wrapping transition and its discussion (together with that of subsequent studies, 
Ref. HaiMI) will be relegated to Section 2.3. 

The four studies mentioned above agreed in that respect that overcharging should 
be a robust phenomenon occurring in these systems but there was still a transparent 
argument missing that would clarify the nature of the underlying mechanism that 
leads to overcharging. Nguyen and Shklovskii |1H] bridged this gap by showing that 
correlations between the charged monomers induced by the repulsion between the turns 
of the wrapped chain can be considered as the basis for this effect. They studied again 
a fully flexible chain and neglected the entropy of the chain conflgurations. The chain 
is assumed to be in the one-tail conflguration with the tail radially extending from the 
sphere. Then - as it is the case in Ref. - the energy of the chain-sphere complex is 
completely given by the electrostatic interactions between the different parts: 



I use here the same symbols as in the previous sections (cf. beginning of Section 2.2.1; 
liso = hZ denotes again the isoelectric wrapping length) . The flrst term in Eq. ^1 is 
the charging energy of the complex (the sphere plus the wrapped chain of length /), the 




(15) 
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second term is the self energy of the wrapped chain portion and will be discussed in 
detail below. The third term is the self-energy of the tail of length L — I, and the forth 
term accounts for the interaction between the complex and the tail. 

I discuss now the second term in Eq. following the arguments given in Ref . jlH] • 
The length A denotes the typical distance between neighboring turns of the wrapped 
chain, i.e. A ^ Rq/^- Consider an isoelectric complex, / = Zj^o, and assume that 
A ^ i?o (multiple turns). Pick an arbitrary charged monomer on the wrapped chain. It 
"feels" the presence of other neighboring charged monomers up to a typical distance A 
beyond which the chain charges are screened by the oppositely charged background of the 
sphere. Hence the wrapped portion of the chain can be " divided" into fractions of length 
A that behave essentially like rods of that length and of radius r having a self energy 
~ /fi (A/6^) In (A/r). One has //A such portions leading indeed to Ib {l/b^)lia.{A/r). 
Another interpretation of this term can be given as follows (again following Ref. |48j ) : 
One can consider the formation of the complex as a two-step process. First one brings in 
sections of length Rq and places them on the sphere in random positions and orientations. 
This leads to a self energy ~ Ib {Rq/IP) In (Ro/r) whereas the interaction of each segment 
with the random background charge can be neglected. Then, as the second step, 
one reorients and shifts these pieces on the ball in order to minimize their mutual 
electrostatic repulsion, i.e., one forms something like an equidistant coil with distance A 
between the turns. Now there is an additional contribution stemming from the attraction 
between each chain piece of length Rq with a stripe on the sphere of length Rq and width 
A leading to a gain in the electrostatic energy, scahng as —Is {Ro/b"^) In {Rq/ A). This 
contribution from all these i?o-sections {I/Rq pieces leading to — (//6^) In (i^o/^^)) 
constitutes the correlation energy of the wrapped chain. Together with the self energy 
of these pieces {I/Rq times Ib {Rq/IP') In {Rq/t)) the correlation leads to Ib {l/b^) In {^/f) 
which is indeed the second term of Eq. [T31 

What is now the prediction of Eq. WbV Minimization of E with respect to / leads 
to the following condition jlHl 

{l-kso)i— ^ ]=^^{- ^1 +2^1nf— 1 (16) 

On the right hand side the argument of the logarithm was simplified assuming L—l ^ Rq 
(long tail), / 3> -Ro (many turns of the wrapped chain) and / ~ kso- Eq. [121 can be 
interpreted as follows: The left hand side describes the cost (if / > liso) of bringing in 
a chain segment from the tip of the tail to the surface, the simplified term on the right 
hand side is the gain in correlation energy (cf. Eq. ITH|l . This leads to the following 
optimal wrapping length: 

r -hso + Roln(^-^^ (17) 

which demonstrates that the correlations induce indeed an overcharging of the complex. 

Note that Eq. [T7| gives the asymptotic value of /* for long chains, L ^ I*, the case 
where most of the monomers are located in the tail. As discussed in detail in Ref. |48J one 
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encounters a discontinuous transition from the one-tail configuration to tlie completely 
collapsed state when one decreases L to such a value that the length of the tail is 
just of order Rq (up to a logarithmic factor). This collapse is similar to the collapse 
discussed in Ref. jlS] but occurs for much shorter chains (namely for chains of the length 
liso + -Ro In (...))• The authors also considered the two-tail configuration and showed 
that it it is formed, again in a discontinuous fashion, for very large chains of length 
L > Ifg^/RQ. This might be - contrary to the claim in Ref. [iB| - in qualitative agreement 
with the prediction of Ref. [S] who found - for sufficiently large spheres - with increasing 
L two discontinuous jumps, from the collapsed to the one-tail configuration and from 
the one-tail to the two-tail state (cf. Fig. 1 in that paper). 

It is also worth mentioning that this theory can be easily extended to semifiexible 
chains. One has to add Eq. Q to the free energy El and finds (following the steps that 
led to Eq. IT7|l for the optimal wrapping length 



i.e., the mechanical resistance of the chain against bending decreases the wrapped 
amount. 

In Refs. jini and [SIl Nguyen and Shklovskii also considered the many-sphere 
case for weakly charged components. Similar to the case considered in the previous 
section they found in the case of an abundance of spheres undercharged complexes (a 
phenomenon which they call "polyelectrolyte charge inversion" as opposed to "sphere 
charge inversion"). It was shown that in this case the imbalance is mainly caused by 
the reduction of the self-energy of each complexed sphere (cf. Ref. jlH] for details). 

Nguyen and Shklovskii argue that the correlation effect is the basis of all the 
phenomena discussed in Section 2.2. The approximation given in Ref. [H], for instance, 
is to homogeneously smear out the charges of the wrapped chain and in this respect 
it overestimates the gain in electrostatic energy upon complexation, i.e., the second 
term in Eq. El is neglected. They call this approximation the " metallization approach" 
|52j . an approximation that obviously holds for sufficiently tight wrapping only. Also 
in that reference they argue that the overcharging via counterion release as discussed 
in Refs. jlH] (cf- Section 2.2.1) is ultimately based on the correlation effect. Their 
argument is as follows: Consider an isoelectric single-chain complex and assume that 
the chain is wrapped in a random fashion around the sphere. Then the electrical field 
close to the complex is essentially vanishing. If more chain would be wrapped around 
the complex no counterions would be released. It is only the fact that the chain will be 
adsorbed in an orderly fashion due to its self-repulsion that each section is surrounded 
by a correlation hole that leads to counterion release even beyond the isoelectric point. 

Sphere-chain complexes have also been considered in several computer simulations 



ISiinniinilinainilMlinSlinniinilinHlinnilZnilZIl- WaUln and Ilnse ISH studied the effect 



of chain flexibility on the geometry of a complex of a single sphere with a polyelectrolyte 
- we will come back to this problem in Section 2.3. The same authors also varied the 
line charge density of the polyelectrolyte ^]] and the radius of the sphere jH!]; flnally 
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they considered the case when there are many chains present [HI]- Chodanowski and 
Stoll [U^ considered the complexation of a flexible chain on a sphere (assuming Debye- 
Hiickel interaction) and found good agreement with Ref. concerning overcharging 
and the discontinuous transition to the one-tail configuration for longer chains. The 
case of multisphere adsorption was studied by Jonsson and Linse, having flexible [Hlj or 
semiflexible chains and taking explicitly into account the counterions of the spheres 
and the chain. Their findings show the same qualitative features as discussed in this 
section above and in Section 2.2.2. A recent study by Akinchina and Linse jZO] focused 
systematically on the role of chain flexibility on the structure of the complex; the results 
will be discussed below in Section 2.3. Messina, Holm and Kremer [HZlinH] demonstrated 
that in the case of strong electrostatic coupling (large values of Ib) it is even possible 
that a polyelectrolyte chain forms a complex with a sphere that carries a charge of the 
same sign - a process which is made possible by correlation effects making use of the 
neutralizing counterions. Most recently Dzubiella, Moreira and Pincus [7T] studied the 
polarizibility of overall neutral chain-sphere complexes in electrical fields as well as the 
interaction between two such complexes. 



2.2.4- Physiological conditions Up to here I discussed only sphere-chain models for the 
case of weak screening, kRq < 1. At physiological condition, however, the screening 
length is roughly 10 nm and hence ten times smaller than the overall diameter of the 
nucleosome. This section is devoted to this case {kRq ^ 1). 

I will mainly focus here on weakly charged chains and spheres where linear Debye- 
Hiickel theory can be applied. For strong screening, kRq 3> 1 the potential (psphere close 
to the ball looks essentially like that of a charged plane with charge density Z/ (AnRl): 
G(psphere{h) /ksT ~ {IbZ /2kE^) e~'^^ {h\ height above surface). Neighboring turns of 
the adsorbed chain have locally the geometry of (weakly) charged rods for which it has 
been predicted that they form a lamellar phase |115t I116j . The lamellar spacing A 
follows from the competition between the chain-sphere attraction and the chain-chain 
repulsion. The chain-sphere attraction leads to the following adsorption energy per area: 

f chain— sphere ^^C\\ 

assuming that the chain is so thin that its adsorbed charged monomers feel an unscreened 
attraction to the surface, Kr <^ 1. To calculate the rod-rod repulsion one starts from the 
potential around a single rod: ecpchain (R) /ksT = —2lBb~^Ko (kR) {R: radial distance 
from rod axis); Kq denotes the modified Bessel function that has the asymptotics 

1 /2 

Ko{x) ~ — Inx for a; <^ 1 and Kq{x) ^ {n/2x) exp (— x) for x ^ 1. This leads 
to the following free energy density of the chain-chain repulsion: 
f 97 °° 

T^^J^T.JM'^'^^) (20) 

To proceed further one might consider two limiting cases. If the lamellar spacing is 
much smaller than the screening length, kA <^ 1, the sum in Eq. EDI can be replaced by 



CONTENTS 



21 



an integral |115t llltij : 



f chain _ , 2/b r _ (21) 



The free energy density / = f chain-sphere + f chain is then minimized for the isoelectric 
lamellar spacing: A = 47ri?o/ (bZ) leading to the wrapping length / ~ AttRq/A = 
bZ = liso- Note, however, that going to the continuous limit means to smear out 
the charges, neglecting the correlation energy discussed after Eq. ^1 This is similar 
to the approximation used by Mateescu, Jeppesen and Pincus jH] who studied the 
unscreened case (discussed above in Section 2.2.3). There they considered, however, the 
self-energy of the tails, a contribution driving more chain monomers to the sphere and 
hence leading to overcharging. This overcharging was overestimated since the correlation 
effects (included in the theory of Nguyen and Shklovskii jlH], cf. also Eq. IT^ were 
washed out. Here, on the other hand, for strong screening the self-energy of a chain 
section remains the same, whether it is adsorbed or not (on a length scale it always 
looks straight). Hence it is here appropriate not to include the tail contribution. 

This might lead one to expect that there is no overcharging for the case of strong 
screening. However, as mentioned above, there is an approximation involved when 
going from Eq. 1211 to |^ This approximation is only good for A ~ r. A more careful 
calculation leads also here to the prediction of overcharging. To see this one has to 
realize that / (knA) dk — J2k (knA) ^ Jq Kq (kA) d/c ~ — In (kA). Taking this 
into account one can replace Eq. |^ by 

fchain ^ T^^B ( ^ , 2 

Minimizing the free energy density with this additional contribution (coming from 
correlation effects) leads to a slightly smaller lamellar spacing 

A ^ ^ (23) 

and to a wrapping length that is larger than the isoelectric one (overcharging): 

r.lfi.,„.«^.ln(y^) (24) 

Lowering the ionic strength leads to smaller ^-values and hence to a reduction of the 
degree of overcharging. When Rq one recovers Eq. El the result presented by 

Nguyen and Shklovskii 0H| for the case of weak screening. 

Netz and Joanny [IH] also considered the case of spheres with an even smaller charge 
density Z/AtcRq where A > k^-^. In that case only the interactions with the two next 
neighboring turns count. From Eq. [201 follows 



2 \ 

l + -/tAln(KA) (22) 

TT / 



fchain \'27iIb _kA (0K\ 

1^ - 62^1/2^3/2^ ^^^^ 

Minimizing / = f chain-sphere + fchain leads then approximately to 

A [in + l) (26) 
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i.e., the spacing is of the order of the screening length (up to logarithmic corrections). 
The overcharging can then become very large. Clearly the term "overcharging" becomes 
quite questionable when there is such a strong screening that charges in the complex 
interact only very locally over length scales of order -C Rq. 

So far the bending energy was not accounted for, i.e., the chain was assumed to be 
perfectly flexible. Bending leads to an additional energy (per area): fbend — Ipl (2-Ro^). 
This contribution scales with the lamellar spacing as 1/A as also the chain-sphere 
attraction, Eq. El does. One can therefore interpret the bending to renormalize the 
sphere charge to a smaller value Z = Z — Iphn/ls- In fact, liso in Eq. |^ has to be 
replaced by Uso — Iph'^K/ls- 

When going to highly charged systems one encounters nonlinear screening 
(counterion condensation). This case has been extensively discussed by Nguyen, 
Grosberg and Shklovskii |117j . They showed that the case A > corresponds 
essentially to the above described case of strongly screened lamellas, Eq. UHl but with 
an effective rod line charge density 1/Ib (instead of the bare value 1/b). There is nearly 
no counterion release in this case. The other limit A ^ is highly non-linear and 
quite complicated. For details I refer the reader to Ref. |117j (section VI in that paper). 

Concluding, what can be learned from the sphere-chain systems with regard to the 
nucleosome structure? It is clear that many of the assumptions entering the models 
do not agree with the details of the nucleosome conformation. Most of the studies 
discussed in Sections 2.2.1 to 2.2.3 assume a weak screening which does not correspond 
to physiological conditions. Many studies neglect the bending energy which is a major 
energetic penalty for DNA wrapping around the octamer. Hydrogen bonds, solvent 
effects are not included. Most importantly, these studies do not account for the fact 
that the wrapping path is clearly prescribed by more or less specific binding patches 
on the octamer surface (cf. Section 2.1). And the octamer is, of course, not a sphere. 
Even when going to the strong screening case, discussed above, the assumed lamellar 
arrangement of the wrapped chain has only a vague resemblance to the 1 and 3/4 turns 
of wrapped nucleosomal DNA. Here one might at least say that the ~ 28A pitch of the 
superhelical ramp is in rough agreement with the prediction A ^ given in Eq. 1^ 

Nevertheless, the improved understanding of the chain-sphere complexes that has 
been achieved in the recent years is in my opinion very helpful. Many of the investigated 
model systems resemble closely complexes between macroions (colloids |118| I119j . 
dendrimers |12()j . charged micelles |121j etc.) and synthetic polymers which is of 
technological relevance as a means of modifying macroion solution behavior. But also for 
the nucleosome the sphere-chain systems help to understand its behavior under changing 
ionic conditions. Especially when the ionic strength is lowered the electrostatics becomes 
long ranged. In this case many local details become overruled by the electrostatics. The 
behavior of the nucleosome at low and high ionic strength will be the subject of the 
next section. 
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2.3. Unwrapping transition 

The nucleosomal complex is only stable at intermediate salt concentrations. In Section 
2.3.1 I will give an overview over experimental results that shed some light on the 
instabilities of the nucleosome core particle at low and high ionic strengths. I also report 
on how these instabilities can be understood in the framework of a sphere-chain model 
(for short chain length). Then I provide a thorough discussion of the instabilities of 
sphere-chain complexes at high ionic strength in Section 2.3.2 and at low ionic strength 
in Section 2.3.3 - both, for chains of arbitrary chain length. 

2.3.1. Instabilities of the nucleosome core particle at low and at high ionic strength 
Yager, McMurray and van Holde characterized the stability of the nucleosome 
core particle as a function of the salt concentration (NaCl) and the concentration of 
core particles (measured via the 260 nm absorbance). Using a variety of experimental 
methods (velocity sedimentation, gel exclusion chromatography and gel electrophoresis) 
they arrived at the following main conclusions (cf. also the schematic phase diagram. 
Fig. 4). For not too low concentrations core particles are stable for ionic strengths 
ranging from 2 mM to 750 mM (called region 1 in that paper; this includes physiological 
relevant salt concentrations ~ 100 mM). For slightly higher salt concentrations (region 
3 jZ2]) or low concentration of core particles (region 2) the DNA is partially dissociated; 
an equilibrium between histones, free DNA and core particles is observed. Especially, 
in region 3 there is also the occurrence of (H2A-H2B)-depleted particles. At salt 
concentrations beyond 1.5 M the core particle is completely dissociated into histone 
oligomers (the (H2A-H2B)-dimer and the (H3-H4)2-tetramer) and free DNA (region 4). 
On the other end, for very low salt concentration < 1 mM (called region IE in j72]) one 
finds an "expanded" state of the nucleosome. 

Using a different experimental approach (measurement of the fluorescence of the 
aa tyromisine in the histone proteins) Khrapunov et al. [H] came to similar conclusions 
(cf. Fig. 4): For ionic strengths between 5 and 600 mM the core particle is intact 
but one finds different degrees of contact between the core histones (the resulting two 
forms are called Ni and N2 in :8.J). At larger ionic strength (^ 1.2 M) the terminal 
regions of the DNA unwrap and the (H2A-H2B)-dimers are dissociated (N4) and at an 
even larger value (~ 1.5 M) the (H3-H4)2-tetramer leaves the DNA. Finally, at low salt 
concentations one encounters an open state (called the Ns-form): the dimers break their 
contact with the tetramer and the DNA termini unwrap (with the dimers still attached 
to them). 

The key features of the behavior of core particle DNA (neglecting the substructure 
of the octamer) are indeed recovered in the framework of the sphere-chain models. 
Most clearly this has been demonstrated in a study by Kunze and Netz |47| IH^ (cf. 
also ^3B"). They considered the complexation of a charged, semiflexible chain with an 
oppositely charged sphere. All charges in their system interact via a standard Debye- 
Hiickel potential. The optimal DNA configuration r (s) (0 < s < L) follows then from 
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Figure 4. The effects of salt on the conformation and stability of nucleosome core 
particles as a function of the salt concentration Cs and the DNA concentration measured 
by the 260 nm adsorbance (adapted from Ref. I^). The notations of Ref. for the 
different states are given as well as the ones of Ref. [S] (in brackets). 



the minimization of the free energy functional 

J^{r(s)} IP ^ ( ^ V IbZ rL e-^{\r{s)\-Ro) 




ksT 2 Jo \R{s) J b{l + KRo)Jo |r(s) 

g-«(|r(s)|-|r(s')l) 



(27) 

|r (sj — r (s'jl 

The first term on the right hand side is the bending energy of the chain where R{s) 
denotes its curvature at point s along the contour. The second and third term account 
for the electrostatic attraction between monomers and the sphere and the monomer- 
monomer repulsion, respectively. The symbols are the same as used above (cf. Section 
2.2.1). Kunze and Netz chose the parameters [Ip, L, b and i?o) such as to mimic the 
values of the core particle and varied Z and as "free" parameters. The optimal shape 
r (s) was found by numerical minimization of Eq. |2ZI and characterized by two order 
parameters, a rotational and a torsional one. They found the following overall picture: 
For reasonable values of the sphere charge Z one finds for vanishing ionic strength 
(k^^ oo) an open, planar configuration where only a small fraction of the chain is 
wrapped whereas the two tails (of equal length) are extended into roughly opposite 
directions. This is reminiscent of the open structures reported in the experimental 



CONTENTS 



25 



studies (region IE in [72], Na-form in jH]). Upon addition of salt the structure stays first 
extended but looses at some point its rotational geometry (in form of a transition from 
a two- to a one-tail configuration) and then at even higher ionic strength the chain goes 
from a planar to a non-planar configuration. It begins then to wrap more and more 
and - at some point - regains its rotational symmetry (the wrapping path resembles 
then some kind of "tennisball seam pattern"). At that point already (which is well 
below physiological ionic strengths) the chain is almost completely wrapped. It stays in 
this wrapped state up to very high salt concentrations. Only then the chain unwraps 
in a discontinuous fashion because the chain-sphere attraction is sufficiently screened. 
Also these features of the complex (the wrapped compact state in a wide range around 
physiological conditions and the unwrapping at high salt content) reflect findings in 
Refs. jSj and jES!- The behavior of such complexes for chains of arbitrary length is the 
subject of the next two sections. 

2.3.2. The rosette state at high ionic strength For physiological conditions (or for even 
higher salt concentrations) the electrostatic interaction between the DNA chain and the 
octamer can be considered as short-ranged. It is then usually sufficient to assume some 
attractive short-range interaction between the chain and the octamer with a range of 
interaction <^ Rq. In this spirit Marky and Manning [ZS] considered the wrapping 
of a semiflexible chain around a cylinder. They came up with the picture of a simple 
"all or none" unwrapping transition. Denote with A the adsorption energy per length 
of the chain on the cylinder (in units of fc^T). Then for each additional wrapped length 
AZ, one gains AEads/ksT = —XAl. On the other hand, in order to wrap the chain it 
has to be bent with a curvature Rq^ with Rq being the radius of the octamer; this leads 
to an energetic cost /S.E elastic /ksT = lpAl/2Rl, cf. Eq. [T] From this follows that if the 
adhesion energy A is larger than 

Ac = (28) 
' 2Rl ^ ' 

then more and more chain will wrap around the histone spool (up to a point when 
there are no adsorption sites available anymore). On the hand, for A < Ac, the chain 
unwraps completely. This unwrapping transition can also be induced by increasing the 
persistence length of the chain, cf. Eq. |2H1 

This, however, is not the complete picture. Already the numerical study by Wallin 
and Linse jSH] indicated that with increasing chain stiffness one encounters a gradual 
change of the conformations of the complexed chain towards more extended structures. 
In Ref. [73] myself, Rudnick, Bruinsma and Gelbart showed in a systematic analytical 
study that there is indeed a wide range of parameters in which more open, multi- 
leafed ("rosette") states occur in this system. The results of this model study together 
with some additional material will be presented in the current section. Following that 
reference I will first discuss the ground-state configurations of the sphere-chain complex 
and then account for thermal fluctuations. Then I will present the general phase diagram 
that includes wrapped (A > Ac) and open (A < Ac) structures. 
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In Ref. we started from the popular worm-like chain model (WLC) which 
provides a good description of the mechanical properties of DNA (for reviews cf. 
Refs. |122[ ll2Hj ). The chain molecule is represented by a semiflexible tube of radius 
r characterized by two elastic moduli, the bending and torsional stiffnesses. The elastic 
energy of a WLC of length L can be expressed as 

A 



1 fL 

Eelastic ~ 2 




(29) 



Here A is the bending stiffness and 1/R{s) the curvature of the chain at point s along its 
contour. The stiffness is usually expressed as A = ksTlp where Ip is the orientational 
persistence length of the chain - as given above in Section 2.1, cf. Eq. ^ The torsional 
angle of the chain is and the torsional stiffness is C, for DNA C ~ kpT x 750A ^99j. 
In addition to this bending contribution there is a short-ranged attraction between the 
chain and the sphere (or cylinder) with a range of interaction 5 <^r. 

As mentioned above for strong attraction A > Ac the chain is wrapped around the 
sphere. On the other hand, if A < Ac, the WLC can only make point contacts to the 
sphere. The energy of a contact point (in units of kpT) 



fi ~ X^/r^ (30) 

follows from the length \/Rq5 of chain portion around that point contact that is located 
within the distance 5 from the sphere. In Ref. ^S] I gave the quantities A and 5 in 
terms of strongly screened electrostatics. In that case 5 = k^^. The adsorption energy 
per length can be estimated from the Debye-Hiickel electrostatic potential close to the 
surface (cf. beginning of Section 2.2.4). From Eq . IT^ follows 

and hence the unwrapping into the rosette takes place when A reaches the critical value 
given by Eq. |2H1 i-e., when the persistence length 

(32, 

is reached. That this is the upper bound for Ip for having a stable wrapped complex 
has been predicted by Netz and Joanny jlHI (cf. Eq. (35) in that paper). 

For A > Ac the optimal number M* of point contacts as well as the preferred 
configuration of the chain is obtained by a minimization of the energy 

E = Eelastic - keTfiM (33) 

with Eelastic givcu by Eq. 1201 In Ref. [SS] we searched first for the minima of the 
elastic energy (zero temperature conformations), neglecting thermal fluctuations that 
were included in a second step (see below). 

The search can be performed systematically by applying the Kirchhoff analogy 
[124j 1125] which relates stationary points of the WLC energy, Eq. 1221 to the well- 
studied classical mechanics problem of the trajectory of a supported, symmetric spinning 
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top in a gravity field. It can be easily shown that the action of the spinning top 
has precisely the same form as Eq. [221 with the time playing the role of the arc 
length, the orientation of the top corresponding to the tangent vector of the WLC, 
the gravitational force being a tension acting on the rod§ etc. This analogy has 
been repeatedly applied to DNA related problems during the last 20 years (e.g. see 
WTA 1123 CSni CSIl IIS2l CSa cm ESS)- For a mce visual review on the 
spinning top-elastic rod analogy the reader is referred to Ref . |125j . 

In Ref. jZl| we presented solutions of the corresponding Euler-Lagrange equation 
that combine the following features: (z) The WLC closes on itself ||; (m) it is possible 
to inscribe a sphere of radius Rq inside the WLC that touches the WLC at M points; 
{iii) there is no self-intersection of the WLC chain with itself if it is surrounded by 
a tube of radius r; and {iv) the solution is stable against small perturbations. The 
resulting rosette-type configurations can be characterized by the number of loops M. 
Fig. 5 shows such a rosette (computed numerically); as indicated in the figure it is 
indeed possible to inscribe a sphere (or a cylinder) in the central hole of the rosette. 
For each M, we adjusted the linking number of the loop to minimize the elastic energy. 
By varying the degree of spatial distortion (characterized by the so-called "writhe") a 
family of solutions was obtained, for given M, with different hole diameters. Solutions 
with the maximum amount of writhe have the smallest central hole diameter as well as 
the lowest elastic energy. The inset of Fig. 5 shows the elastic energy of a loop of length 
L, in dimensionless units, for the rosette state, as a function of the degree of writhe. 

These solutions are actually saddle points of the WLC energy, i.e., there is a finite 
subset of infinitesimal distortions that lower the elastic energy of the WLC; the rest 
raises the energy or leaves it unchanged. The role of the sphere is to stabilize this 
saddlepoint and turn it into a real maximum (for a detailed discussion cf. Ref. |133j V 
The energy of a minimum-hole rosette depends on the overall chain length L and the 
number of point contacts M (= number of leaves) as 



with X = 7.02 (see inset of Fig. 5). This result can be understood making use of results 
of the earlier work by Yamakawa and Stockmayer (YS) jl36j who showed that a loop of 
length /, formed by imposing common endpoints on a WLC strand, assumes the form 
of a lemniscate-shaped leaf with an 81-degrees apex angle. The elastic bending energy 
of a leaf is e (/) = 2xA/l and Emin{,M) given above just equals Me{L/M) plus the 
adhesion energy. We verified numerically that the leafs of the rosette indeed have apex 

§ The tension comes here from the "sticky" sphere that induces the rosette structure discussed below. 
Later, in Section 2.4.2, I wiU give an example where it becomes more obvious how the rod tension 
formally enters the Hamiltionian in form of a Lagrange multiplier T - leading to a term that resembles 
the potential energy of the spinning top. 

II This assumption is only for technical reasons. As discussed below, the behavior for an open chain is 
essentially identical. The solutions for closed WLC can be characterized in terms of the topologically 
conserved linking number. Loops with a non-zero linking number show a configuration that combines 
twist and spatial distortion (known as "writhing"). 



Emin (M) 2xlpM^ 
keT ~ L 



- nM 



(34) 
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Figure 5. Five-leafed rosette. This configuration corresponds to the minimal energy 
solution of the Euler-Lagrange equation with a maximum amount of writhe and the 
smallest central hole (cf. text for details) 



angles close to 81 degrees. The energy Emin (M) exhibits a minimum as a function of 
the number M of rosette leaves for 

We thereby obtained an open, multileafed structure controlled by the adhesion 
energy and the persistence length that competes with the wrapped state. Packing 

3 /2 

considerations imply an upper limit for the number of loops of order Mmax ~ (Ro/^) ! 
this is the maximal number of contacts, each excluding an area ~ ry/Riyr, that can be 
closely packed on the surface of the sphere. 

Note that the above results do not depend critically on the assumption that the 
chain forms a closed loop. In fact, for a chain with open ends we expect the same number 
of rosette leaves, each again having approximately the shape of a YS-loop with an 81° 
apex angle. This means that here neighboring leaves will have a relative orientation of 
~ 180° — 81° = 99°. In addition, the leaves will be slightly twisted with respect to each 
other (like propeller blades) to account for the mutually excluded volume. I give an 
overview over the results in a "phase diagram" (ground state configurations) in Fig. 6 
(adapted from Ref. jSS]; an alternative presentation with a different choice of axes is 
presented in Ref. |7i]). 

We next studied in Ref. [H] the stability of the rosette against thermal fluctuations. 
We started from a single, large loop of length L and constructed the rosette step by 
step, by attaching to the sphere lemniscate-shaped leafs of variable length of the kind 
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Figure 6. The sphere-chain complex in the case of short-ranged attraction (for 
instance, at high ionic strength). Depicted is the diagram of states as a function 
of the total length L of the chain and its persistence length Ip divided by the point 
contact energy ^ (both axes in units of Uso)- The thick vertical line indicates the sharp 
unwrapping transition from the wrapped to the rosette-type complexes. 



examined by YS |136j . The finite-temperature free energy cost F {keaf, (p) of introducing 
a single leaf of length keaf and apex angle into a large strand was computed by YS. 
Using path-integral methods they found |136j 



The function W (0) has a minimum when the apex angle of the leaf is approximately 
81°. The logarithmic contribution to f{lieaf,(p) ^ associated with the configurational 
entropy of the loop neglecting excluded volume interaction {6 solvent) - imposes a free 
energy penalty for large leaves, keaf ^ ip (corresponding to the entropy of a closed 
random walk in 3 dimensions). The enthalpic l/keaf contribution imposes an energy 
penalty for small leaves. Hence, for given 0, / {heafy 0) as a function of lieaf has a 
shallow minimum, its value being close to that of the persistence length Ip. The total 
free energy cost Fm ({/,;}) of introducing into a large loop an M-leafed rosette for a fixed 
distribution {/j} of leaf lengths is then given by 




(36) 



Fm {{h}) 
keT 



M 



fih) 



E 



(37) 



i=l 



Here 




(38) 
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(^i)"" (39) 



constitutes an interpolation formula between the large and small / limits of / (/, 81°) 
as given by Eq. EHl The partition function Zm for an M-leafed rosette follows by 
integration over all possible leaf distributions 

dh 

p^cr^ — 

Assuming /i 3> 1 (strong sticking points) one can assume the number of leaves always to 
be maximal, M = M^ax (the index is dropped here and in the following for simplicity 
of notation). In Eq. EHla Lagrangian multiplier P is introduced in order to satisfy the 
constraint Y^iLi h = L. Physically, P is the overall tension of the loops induced by their 
adhesion to the sphere. The free energy follows to be 



r ^ dk 

/ 11 — exp 

^0 fJi r 




G (P) = -kBTlnZM = MksT { 2^ + In | ) } (40) 

It is interesting to note that G (P) is mathematically identical to the free energy 
of a one- dimensional many-body system of M particles under a "pressure" P confined 
to a circular track of length L. The particles are interacting via a "nearest-neighbor 
pair-potential" / (/) while ksT^ is the "chemical potential" of the particles. According 
to Eq. EHl the effective pair potential / (/) is concave (i.e., d? f /dt^ < 0) for inter-particle 
spacings exceeding a spinodal threshold spacing of order Ip. Experience with mean-field 
theory of many-body systems suggests that one should expect phase-decomposition if 
the average spacing between particles exceeds a spinodal threshold. 

Using L = dG/dP one finds the non- linear "tension-extension" curve P (L) = 
2xlpkBT [M / Lf' . Using Eqs. IHUlandHUl it is straightforward to compute the first and 
(reduced) second moments of the leaf size distribution: 

(/zea/) = -fcBT^lnZi = A (41) 

and 



\j {{heaf - ikeaf)) 



^ (42) 



ikeaf) 2 V xMlp 

Thus one encounters no phase-coexistence: the leaf size grows with chain length in 
the same manner as the "T = 0" solution. This does not mean, however, that the 
rosette structure is not altered when the mean leaf size exceeds Ip, because the reduced 
second moment starts to exceed one at that point (actually, at (keaf) exceeding 4x/p). 
One can therefore identify L/M Ip as, the onset point of heterogeneity of the leaf size 
distribution; the orderly, symmetric rosette is starting to "melt" . It must be emphasized 
though that G (P) is analytic and that there is no true thermodynamic singularity. 

This heterogeneous rosette state is very fragile. For instance, the loop size 
distribution will change drastically if one considers an open chain, i.e., if one allows 
for two free ends. Consider a chain of length L with M leaves of length U {i = 1, M) 
and the two chain ends of length Iq and hi+i- One requires Y^iLt^ h = L. The calculation 
of the partition function goes along similar lines as above. One has just to multiply Zm 
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Figure 7. Molten rosettes of closed (left) and open chains (right). These structures 
are extremely fragile: cutting one loop leads to two long tails at the expense of the 
other loops. Shown is also the corresponding problem of particles on a one-dimensional 
closed (or open) track that interact via a nearest-neighbor pair potential. 



in Eq. 1221 with /q°° dlodlM+i exp (— P (/q + hi+i) I^bT) jr^ which accounts for the free 
ends. This leads to G (P), Eq. HOI with an additional additive term 2kBT\ia {rP/kBT). 
The "pressure" has to be chosen such that L = dG/dP = MJ2xlpk^T/P + 2kBT/P, 



I.e. 



P 



2kBT xksTlpM^ M x^bTI 



+ 



L ■ 2L2 
The average loop size is now given by 

d 



+ 



L 



'^BT—ln Zi 
dP 



l2xkBTlp 



P 



{heaf) 

The second moment of the loop size distribution obeys 

1/4 



J! for ^ « M 

V¥FL for^»M 



(43) 
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y {{keaf - {heaf)y 
iheaf) 



L/M 



2 (xtp) 



1/4 



for ^ < M 

for ^ > M 

ip 



(45) 



Hence one recovers the result for the closed- loop case but only for LjM <^ Ml p. In the 
opposite case, L/M ^ Mlp, one finds a different scaling (lieaf) oc \/lpL (Eq. Hlj) which 
shows that the mean leaf size is now small compared to L/M. Most of the chain is part 
of the two free ends that emerge from the rosette; the average length of a free end is 
(Zo) = kpT/P which leads indeed to L/2 — My/xlpL/2. A schematic view of molten 
rosettes formed from closed and open chains and the corresponding analogy of a ID gas 
of particles is shown in Fig. 7. 

In Appendix A I present some new results that extend the above calculations to the 
case of rosette formation in d dimensional space. This will shed some light on the nature 
of the "phase coexistence" within molten rosette structures. In the next section (2.3.3) 
I shall show that rosette structures do also occur at low ionic strength. In that section 
I will also contrast the unwrapping transitions into the rosette at low and high ionic 
strength. Furthermore, I will speculate if rosette structures could occur in DNA-histone 
complexes. 



CONTENTS 



32 



2.3.3. The rosette state at low ionic strength The rosette configuration discussed in 
the last section is a way to bring at least a small fraction of the chain in close contact 
to the ball (in the form of point contacts). The majority of the monomers resides in the 
loops that do not "feel" the presence of the sphere but are needed to connect the point 
contacts via small curvature sections. At first sight one might thus expect the rosettes 
to be a special feature for chain-sphere complexes with a short ranged attraction. 

This is, however, not true. Rosettes are quite robust and occur also in systems with 
a much larger range of interaction. Recently This became clear in a Monte Carlo study 
by Akinchina and Linse They considered the complexation of a semifiexible charged 
chain with an oppositely charged ball that carries the same absolute charge as the 
chain (isoelectric complex). No small ions were present so that the charged monomers 
were attracted to the sphere via a long-ranged 1/r-interaction. The authors simulated 
systems with chains of different persistence lengths and linear charge densities as well 
as spheres with different radii. Depending on the choice of parameters they encountered 
a multitude of structures - ranging from collapsed structures with a "tennisball seam 
pattern" or solenoid arrangement of the wrapped chain |m 131] to open multi-leafed 
structures very much resembling the ones discussed in the previous section. The rosette 
structures occur for stiffer chains on smaller spheres. The example configurations in 
Ref. [Zni clearly show some rosette structures with one, two and three leaves (cf. Fig. 1, 
system II in that paper). That these are representative example configuration can most 
clearly be seen in the adsorption probability of monomers as a function of the monomer 
index (Fig. 3 in jTHp. 

To understand better why rosettes occur also in the long-ranged case I developed a 
scaling theory for this system in Ref. [HEj that I will outline in the following. Consider 
first sufficiently short chains L = bN < bZ = li^o where the chain charge is smaller 
than (or equals) the sphere charge. The energy of the rosette with M leaves is then 
approximately given by 

El^^ ^ ^_Lm^ - (46) 
keT Lb ^ ^ 

The first term is the bending energy of M leaves of length L/M and typical curvature 
~ M/L. This has, of course, the same form as the elastic contribution to the energy of 
the rosette at strong screening, first term of Eq. |^ (up to a numerical constant that I 
do not consider here). The second term is the attraction between the ball charge Z and 
the chain charge L/b over the typical distance L/M. Remarkably this term shows the 
same scaling with M as the second term of Eq. I^that describes the energy of the point 
contacts! One has just to identify the point contact energy /x (for the strong screening 
case) with the leaf-sphere attraction 

^='-f (47) 
for the unscreened case. The optimal leaf number is thus again (cf. Eq. 1^3)) 
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and the leaf size is 

keaf (49) 

The rosette state competes with the wrapped state that was aheady discussed in Section 
2.2. The rosette state is expected to transform continuously into the wrapped structure 
when L/M* ~ Rq; then the leaves become so small that they touch with their contour 
the surface of the sphere. Indeed, by setting M = L/ Rq in Eq. ESI one finds 

Ewrap _ IpL IbZL I Ip IbZ 



kpT Rl bRo \Rl bRoJ^ 

which can be considered as the free energy of the wrapped state: The first term is 
the bending energy, Eq. d and the second accounts for the electrostatic chain-sphere 
attraction. All other electrostatic contributions (as written down in Eq. ITH|l are much 
smaller and do not occur on this level of approximation. 

On the right hand side of Eq. EDI I arranged the terms in such way that one can 
deduce directly an unwrapping transition at Ip/R^ = IpZ/bRo, i.e., at 

Ip = liRo (51) 

At first sight one might expect that at that point the chain unwraps in a strongly 
discontinuous fashion similar to the cases discussed above (the highly charged case, 
Eq. 13 and the short-range case discussed at the beginning of Section 2.3.2). However, 
this "unwrapping point" corresponds just to the point L/M* ~ Rq when loops form on 
the sphere. This leads to a smooth transition as pointed out before Eq. EDI That the 
unwrapping transition occurs rather smoothly at low ionic strength and sharp at high 
ionic strength has been predicted by Netz and Joanny even though the authors did 
not allow in their study for rosette structures. 

To complete the picture one has also to study chains that are longer than the 
isoelectric length, L > Uso- Then at least three terms are needed to capture the essential 
physics of the rosette state: 

^12^ J^M^ -^^M +^^M (52) 
fcsT I b b^ ^ ' 

Here the monomers are "allowed" to distribute between the rosette of length / and a tail 

of length L — l. The first two terms are the same as above, Eq.^^l the last term describes 

the self-repulsion of the monomers that constitute the rosette (additional logarithmic 

terms accounting for the self-energy of the tail and its interaction with the rosette, cf. 

Eq. are smaller and neglected here). Minimization with respect to / leads to the 

optimal rosette length 



The free energy 1^21 with the optimal wrapping length /* is minimized for the following 
number of leaves: 

M* = 4^ (54) 
Ip 
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Figure 8. The sphere-chain complex at low ionic strength. The axes are chosen as in 
Fig. 6 with /i being now the leaf-sphere attraction, Eq. 1471 The unwrapping transition 
is smooth in this case. 



Hence - on this level of approximation - /* ~ h^lpM* /Ib — kso, i-e., the rosette 
monomers just compensate the central ball charge; the rest of the monomers extends 
away from the rosette in a tail of length L — Uao- Each leaf is of size 

keaf ^ Ip/fJ' (55) 

The rosette disappears at Uso/M* = Rq, i.e., when Eq. EUis fulfilled. It is then replaced 
by a wrapped chain of length Zj^o (plus additional correction terms such as the one given 
in Eq. ITHjl and a tail of length L — Uso- 

In Fig. 8 I depict the complete diagram of the sphere-stiff chain complexes to be 
expected in the long-range case. I again plot L vs. x = Ip/fi (in units of liso), which is 
for rosettes just the leaf size, Eq. 0^1 When one starts in this diagram at a large a;- value 
and goes towards smaller values (with some arbitrarily fixed value L < liso) then all 
leaves shrink and more and more leaves can form. At x = -Rq the maximal number of 
leaves (for that given value of L) is reached and at the same time the leaves disappear 
simultaneously in a continuous fashion. For x < Rq the chain wraps around the sphere. 
For L > liso the excess charges are accommodated in tails and all rosettes have the same 
lengh liso- The borderlines between different rosette ground states are then independent 
of the total length of the chain and thus appear as vertical lines. Desorption occurs when 
the free energy, Eq. |3nior|H2l equals the thermal energy kpT. This point is reached when 

h ^ i VlpL for L < liso j-gg-j 
fJ' [ ylpliso for L > liso 

An arrow in Fig. 8 indicates the direction in which desorption takes place. 

Before comparing these results with the Monte Carlo simulations by Akinchina and 
Linse jTU] and with the properties of the nucleosomal complex it is instructive to take a 
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Figure 9. The unwrapping transition for (a) short-ranged and (b) long-ranged 
interaction. In the former case the chain unwraps discontinuously into large-leafed 
Yamakawa-Stockmeyer loops, in the latter case the transition is continuous. 

closer look at the unwrapping transition and to contrast the short- and the long-ranged 
case. The former case, kRq ^ 1, is depicted in Fig. 9(a). As discussed in Section 2.3.2 
the unwrapping transition is expected to occur at Ac — lp/2Rl (Eq. l28|l which leads 
to Eq. 1221 At this point the structure jumps in a strongly discontinuous fashion into 
a /arg' e-leafed rosette with leaves of size lieaf = L/M*. Using Eqs. IHUl and 1^ one finds 
indeed 

keaf ^ Ro\fi^ > Rq (57) 

As discussed in the paragraph after Eq. EHl neighboring leaves have a relative orientation 
of ~ 99° which I also indicated in Fig. 9(a). 

The unwrapping at low ionic strength is depicted in Fig. 9(b) and goes as follows. 
When the chain becomes so stiff that Ip/Rq > IpZ/ (bRo) the wrapped state is not 
stable anymore (cf. Eq. I^UI) . At that time many small leaves (M* = L/Rq ones) form 
simultaneously in a continuous fashion. Their size scales as keaf — Ro, the precise 
prefactor being not accessible to our scaling argument. The typical opening angle 7 of 
the loop at the point of its formation scales as (L/Rq) /M* ^ 1, again with an unknown 
numerical value. A multi-leafed configuration slightly above the unwrapping point is 
depicted in Fig. 9(b). 

Additional insight can be gained by generalizing the attractive force between a 
given chain charge and the sphere by a power law —AZ/ with an arbitrary exponent 
a > 0. An integer value a = d — 2 with d = 3,4, ... can be interpreted as a charged 
chain that adsorbs on a d- dimensional ball in a d-dimensional space. The electrostatic 
term for the rosette in Eq. Intakes then the form —AZN"'/ {bL°'~^) and the one for the 
wrapped state scales as —AZL/BRq . Unwrapping takes place at Ip ~ AZ/ (BRq'"^^. At 
this critical value the energy of the rosette E^osette {Ip = ^p) has (as a function of M) 
a minimum at M* ~ L/Rq for a < 2 {d < 4), suggesting a rather smooth unwrapping 
transition similar to the one depicted in Fig. 9(b). This minimum turns into a maximum 
at a = 2 (li = 4). For larger values of a we find M* = 0, i.e., the unwrapping transition 



CONTENTS 



36 



is sharp, similar to the short-ranged case discussed in the previous section. 

Let me also shortly come back to the highly charged case reviewed in Section 2.2.1. 
In that case the dominant contribution to the complexation energy is the release of 
counterions which is a rather short-ranged interaction and consequently the unwrapping 
transition should be expected to be discontinuous - even at low ionic strength. To 
determine the unwrapping point one has to use Eq. [7| (with the " — " sign) from which 
follows that unwrapping into the rosette occurs at Ip ^ {VL + n) Rl/h. The point 

contact energy is of order {Vt + Vtj ^JR^XQc/h where the so-called Guoy-Chapman length 
Xgc — 1/ (ct/b) is the thickness of layer of condensed counterions around the sphere 
(Agc ^ Rq for strong counterion condensation [SIl)- The leaf size at the unwrapping 
point is then given by 

W-\/^^o>i?o (58) 

V ^GC 

which indeed indicates a sharp unwrapping transition for highly charged systems. 

We compare now the results of this section to the Monte Carlo simulations by 
Akinchina and Linse [70]. As mentioned above the simulated systems were always 
at the isoelectric point, i.e. L = hZ = Uso- Furthermore there was no screening. 
The simulation results have thus to be compared with Fig. 8. Four systems have been 
considered, each having a fixed set of parameters 6, Z and Rq but with 7 different values 
oilp. This means that for each case the systems were located on the dashed horizontal 
line at L = Zj^o in Fig. 8. In one system (called system II in [TUj) the continuous 
development from a wrapped to the rosette configurations has been seen very clearly. 
Example configurations are shown in Fig. 1, system II in that paper. For Ip = 7 A the 
chain is wrapped, at Zp = 60 A there is already a slight indication of very small loops 
(A^ = 4 or 5, cf. the small oscillations in Fig. 3, systems II, open squares). The next 
system depicted has already a much stiffer chain, Ip = 250 A, and shows very clearly 
three leaves, then two leaves at Ip = 500 A and one leaf for the stiffest chain, Ip = 1000 
A. In Fig. 8 I have chosen the parameters such that Amax — kso/ Ro equals 4 so that this 
corresponds roughly to system II in jTO]. To compare with the simulations one has to 
follow the L = /jso-line in Fig. 8: One starts with wrapped structures for x = Ip/ fi < Rq 
and finds then the continuous evolution of rosettes when the line x = Rq is crossed. The 
leaves grow at the expense of their number (first 4, then 3, 2 leaves), just as it has been 
observed in the simulations. The other three systems considered have different sets of 
parameters Rq and b. The observed behavior of these systems is also in good agreement 
with the theoretical predictions. For details I refer the reader to the last section of 
Ref. 1^. 

Finally let me speculate up to what extend the results of the previous and 
the current section apply to the nucleosome system under different ionic conditions. 
Especially let me ask if rosette structures could in principle occur on DNA-histone 
complexes. First consider the core particle (147 bp DNA). As already discussed in 
Section 2.2.4 at physiological salt concentrations (~ 100 mM) one has k^^ ~ 10 A so 
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the short range case of the previous section apphes. However, the estimate for A from 
Eq. I^is not rehable since r ~ 10 A (partial screening), since the binding sites between 
DNA and the histones are quite specific and since hnear Debye-Hiickel theory is not 
rehable for such highly charged components. As mentioned in Section 2.1 A can be 
derived instead experimentally from competitive protein binding to nucleosomal DNA 
P7] to be of order QksT per sticking point, i.e., A ~ (1/5) A~^. This is roughly 5 
times larger than what one would expect from Eq. |^ (but note also that contact is 
made mainly between DNA minor grooves which are 10 basepairs apart!). In any case, 
Eq. 1221 predicts an unwrapping for sufficiently small values of k,~^ but the numbers are 
not reliable. On the other hand, the unscreened long-range case of the current section 
applies when Cg < 1 mM {k~^ > 100 A). As mentioned in Section 2.3.1 the completely 
wrapped configuration is not stable at this point anymore. This is not surprising since 
the nucleosomal DNA overcharges the protein octamer by at least 74 negative charges 
(if not by 160 charges since 86 charged residues are inside the octamer, cf. Section 2.1). 
Hence it is expected that a considerable part of the terminal DNA unwraps and is part 
of one or two tails. In Fig. 8 this corresponds to the wrapped chain structures with tail 
that are found for small values oilp/ ^ < R and large values of L > Uso- 

It would be interesting to redo the experiments |H1 [Z2] mentioned in Section 2.3.1 
with complexes between histone octamers and DNA segments that are longer than 147 
basepairs. For sufficiently large salt concentrations I expect the formation of DNA 
rosettes (if there is no interference with the partial disintegration of the octamer which 
could be avoided by introducing covalent linkages between the core histones). For the 
other limit (low salt) it might be appropriate to use the argument for highly charged 
chains and spheres as given before Eq. EH The linear charge density of DNA is very high 
(2 phosphate groups per bp, i.e., per 3.4 A). Manning theory |lllj indeed predicts that 
counterion condensation reduces the linear charge density to —c/Ib with Ip = 7 A. Also 
the charge of the histone octamer is so high that counterion condensation is important. 
It was argued above that the unwrapping occurs around Ip = (jl + ^Ij Rl/b {^l and f2 
are numbers of order (but larger than) one). However, since the DNA persistence length 
is so small that Ip < Rl/b (50 A)^/1.7 A^ 1500 A I expect the wrapped state to be 
stable at low ionic strength and there should be no formation of rosettes in this limit. 

2.4- Nudeosome repositioning 

It has been shown that nucleosome repositioning occurs spontaneously via thermal 
fluctuations (under certain conditions). This autonomous repositioning is the subject 
of the current section. I flrst review the relevant experiments (2.4.1) and then discuss 
three theoretical models proposed to account for this effect: bulge diffusion (2.4.2), 
large loop repositioning (2.4.3) and twist diffusion (2.4.4). I contrast the three cases 
in Section 2.4.5 and speculate that similar modes might be catalyzed by remodeling 
complexes that use the energy of ATP hydrolysis. Nucleosome repositioning induced 
through transcription on short DNA segments is also briefly discussed in that section. 
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2.4-1 ■ Experiments An early study of uncatalyzed nucleosome repositioning was 
presented by Beard |137j . He constructed " chromatin-DNA-hybrids" where segments of 
radioactively labelled naked DNA were covalently joined to sections of chromatin fibers 
(derived from a simian virus). Remarkably the experiments suggested that nucleosomes 
spontaneously reposition themselves and invade the naked DNA. This was shown by 
several experimental methods: (z) the occurrence of a radioactive component in the 
~ 175 bp-band of a gel electrophoresis experiment that separated the products of 
a nuclease digestion, (m) Electron- micrographs showed an increased spacing of the 
nucleosomes on the chromatin pieces and {Hi) binding of radioactive components of the 
recut hybrids on nitrocellulose filters (to which naked DNA would not bind under the 
given conditions). All these methods indicated that something "happens" on the time 
scale of hours in a 150 mM NaCl solution at elevated temperatures (37°C). The methods 
used gave, however, not more quantitative information about the nucleosome migration 
rate. Furthermore, the conditions were not well-defined enough to exclude ATP-driven 
processes. Finally, the histone tails belonging to the virial DNA were strongly acetylated 
compared to the "normal" chromatin of its host cell (which - as shown much later - 
has not a big effects on nucleosome mobility, see below). 

In that paper Beard discussed also three different possible modes of nucleosome 
repositioning: (a) jumping where the octamer dissociates completely from the DNA 
and complexes at another position, (b) sliding or rolling in which the octamer moves 
along the DNA without dissociation and (c) displacement transfer where a naked 
segment of the DNA displaces a nucleosomal DNA section and takes over. Beard 
argues that two of the mechanism (jumping and displacement transfer) imply that 
in the presence of competing DNA nucleosomes would be transferred from one DNA 
chain to another, a fact that was not observed in his experiment when chromosomes 
and radioactively labelled, naked DNA were mixed together. Nucleosomes could only 
migrate onto the naked DNA when the chromosome and DNA pieces were covalently 
joined. Beard concluded that "sliding and rolling modes" should be responsible for 
octamer repositioning. 

Spadafora, Oudet and Chambon |138j showed via gel electrophoresis that 
nucleosomal rearrangement occurs in fragments of rat liver chromatin when certain 
conditions are fulfilled: Either one needs to go to high ionic strength above 600mM 
NaCl or the fiber has to be depleted of the linker histone HI; in that case rearrangement 
occurred also around physiological conditions. It was observed that when either of these 
conditions is fulfilled and if temperatures are elevated (again 37°C), then nucleosomes 
seem to move closer to each other (in the time range of hours), away from the natural 
200 bp repeat length towards a ~140 bp repeat length. That length corresponds to a 
close packing of nucleosomes, a fact that the authors assigned to an internucleosomal 
attraction. Furthermore, since HI is dissociated from the nucleosome at around 600 
mM salt content they concluded that one role of HI is to prevent nucleosome mobility. 

Similar conclusions were also drawn by Watkins and Smerdon |13T^ from 
corresponding experiments on human chromatin. An interesting additional experiment 
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Figure 10. Schematic view of the two-dimensional gel electrophoresis experiment by 
Pennings et al. i29j that allowed to demonstrate autonomous nucleosome repositioning, 
(a) Under conditions where no repositioning takes place the final products line up on 
a diagonal - each spot corresponding to a certain nucleosome position, (b) If the 
mononucleosomes are incubated the final products form a square of dots. 



studied the exchange of histone proteins between different DNA chains. This was shown 
by mixing ^^C labelled chromatin and labelled naked DNA. It became clear that at 
high ionic strength (600 mM) there is a formation of (probably intact) nucleosomes on 
the competing DNA whereas at physiological conditions this was not unambiguously 
detected. Moreover, the double-radioactive component showed the footprint of a ~146 
bp repeat length. This might indicate that a combination of histone protein transfer 
and subsequent nucleosome "sliding" into tight packing took place. 

An important series of experiments on nucleosome repositioning under rather well- 
defined conditions was performed by Pennings, Meersseman and Bradbury [SHI EDI IS]- 
Even though the original focus of this study was to understand better the positioning 
of nucleosomes on special natural sequences that have a high affinity to octamers |140j . 
the authors came up with elegant methods to monitor the nucleosome repositioning. It 
was found [21] that on tandem repeats of 5s rDNA positioning sequences (each of length 
207 bp) nucleosomes assemble in one dominant position surrounded by minor positions 
multiples of 10 bp apart. The most interesting observation was that there is a dynamic 
redistribution between these positions. This was shown by cutting the 207i8 chromatin 
into its repeating subunits and then studying the nucleosome dynamics on such 207 bp 
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fragments. The authors took advantage of the fact that different nucleosome positions 
on the chain give rise to different electrophoretic mobihties^ and that the motion of 
the nucleosome along the chain can be suppressed by subphysiological temperatures or 
ionic strengths, and by the presence of Mg^+. The 207 bp mononucleosomes were first 
separated by an initial dimension of the electrophoresis in conditions where mobility 
is suppressed. In this way they obtained different bands indicating a set of preferred 
positions. An entire track from such a gel was then incubated for some period of 
time in new conditions where mobility may occur, then changed back to physiological 
conditions where mobility is suppressed again, and run in a second, equivalent, dimension 
of gel electrophoresis. Essentially, the first dimension of electrophoresis created a 
nonequilibrium distribution. Depending on the conditions this distribution relaxed 
during the subsequent incubation, which in turn was detected as products moving off 
the diagonal in the second dimension of electrophoresis - cf. Fig. 10 for a schematic 
depiction of the two-dimensional separation technique. 

In Ref. |E| it was found that substantial redistribution took place when the sample 
was incubated for an 1 hour at 37°C but not at 4°C. The experiments were done in 
low ionic conditions in a tris-borate buffer (0.5xTBE; cf. Ref. |144j for a discussion 
of the effects of this buffer on naked DNA). There was a set of preferred positions, all 
multiples of 10 bp (the DNA helical pitch) apart. That means that the nucleosomes 
had all the same rotational positioning with respect to the DNA. Another feature that 
was observed is that the nucleosomes have a preference for a positioning at the ends of 
the DNA fragments, a typical feature for nucleosomes on short DNA that was recently 
discussed by Sakaue et al. jnEI- The 5s rDNA positioning sequence itself, however, is 
located more towards the middle. When gel-separating the mononucleosomes directly 
after they have been excised from the tandemly repeated nucleosomes this position led 
indeed to the strongest band. After incubation, however, the end positions showed the 
highest probability. 

The authors extended their study to head-to-tail dimers of 5S rDNA (2072) PH] - 
In a first dimension of a 2D gel electrophoresis the mononucleosomes were separated 
according to their position on the dimer. This was followed by an incubation at 4° 
or 37° C and a subsequent cutting of the dimer into its monomers. The resulting 
product was then separated through electrophoresis in a second dimension. Again for 
the sample incubated at elevated temperatures a repositioning of the nucleosomes was 
found. Interestingly, however, the study indicated that the repositioning took place only 
within a cluster of positions around each positioning sequence but not between them, a 
fact that was shown by radioactive labelling of one half. This finding indicates that there 
is no "long-range" repositioning at low ionic strength. Other systems studied in that 
paper were fragments of Hl-depleted native chromatin and nucleosomes reconstituted 

^ One reason is the bend induced by the octamer; it is known that a bend on a naked DNA fragment 
affects its mobihty in the gel in very much the same way |141| . Another reason is the inhomogeneity of 
the charge distribution along the chain that traps such a polymer in the gel in a U-like conformation 
(cf. Fig. 3 in Ref. [1^11^ 1 
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on Alu repeats; also in these repositioning was detected result of an 

elevated temperature incubation. The authors concluded that the repositioning "may 
be visualized as following a core screw movement within the superhelical path of the 
DNA" The same authors studied in Ref. [3T] the nucleosome mobility on the 20 72- 
dimer in the presence of linker histone HI (or its avian counterpart H5) and found that 
the mobility of nucleosomes was dramatically reduced. 

Ura et al. |145l I146j - following Ref. IHOj - studied nucleosome mobility on the 
20 72-dimer under varying conditions, namely in the presence of various chromosomal 
proteins and in the case when the core histones are acetylated, respectively. In the 
former case mobility was suppressed (depending on the type and concentration of the 
chromosomal protein), in the latter case the mobility was not changed much. 

Flaus et al. jl47j developed a different strategy to determine nucleosome positioning 
and repositioning. In their method they used a chemically modified H4 histone that 
induces - after addition of some chemical - a cut on the nucleosomal DNA 2 bp away 
from the dyad axis. Via gel electrophoresis of the resulting product they were able to 
determine the nucleosome position with base pair resolution. Using this method Flaus 
and Richmond |148j studied the nucleosome dynamics on a mouse mammary tumor 
sequence which revealed several features of repositioning more clearly. The longest 
fragment of this sequence studied was 438 bp long and had two positioning sequences 
where two nucleosomes assembled, each at a unique position. These positions were also 
found when mononucleosomes were assembled on shorter fragments that included only 
one of the two positioning sequences. The authors studied the degree of repositioning 
of the mononucleosomes on such shorter fragments (namely nucleosome A on a 242 
bp- and nucleosome B on a 219 bp-fragment) as a function of heating time (ranging 
from 20 to 80 minutes) and temperature (ranging from 0° to 50°C). It was found 
that the repositioning rates - as estimated from the occurrence and intensity of new 
bands - increase strongly with temperature but also depend on the positioning sequence 
(and/or length of the fragment). The difference of repositioning for the two sequences 
is remarkable: at 37°C one has to wait ~ 90 minutes for the A242 and more than 30 
hours for the B219 for having half of the material repositioned. Another feature found 
was again a preference for end positions (roughly 70 bp from the dyad axis, similar to 
the finding in Ref. P^). For nucleosome B which showed a slower repositioning the 
set of new positions were all multiples of 10 bp apart (namely at a 20, 30, 40, 50 bp- 
distance from the starting position), i.e., they all had the same rotational phase. On 
the other hand, nucleosome A did not show such a clear preference for the rotational 
positioning. It was argued that these differences reflect specific features of the underlying 
base pair sequences involved. Nucleosome B is complexed with a DNA sequence that 
has AA/AT/TA/TT dinucleotides that show a 10 bp periodicity inducing a bend on the 
DNA whereas nucleosome A is positioned via homonucleotide tracts. 

The authors speculate in Ref. |148j that the preference of end positions might be 
caused by one (or several) of the following mechanisms (1) direct histone interaction with 
a special structure at the DNA terminus, (2) relief of the repulsion between entering 
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and exiting strand and (3) entropy gain by having a long unbound DNA extension. 

Recently Hamiche et al. |149j demonstrated that the (uncatalyzed) nucleosome 
mobility along DNA depends on the presence of histone tails. Especially, in the absence 
of the N-tail of H2B that passes in between the two turns of the nucleosomal DNA 
[10] spontaneous repositioning of the nucleosomes was detected (here during a gel 
electrophoresis in a second direction). 

How does repositioning work? The studies seem to indicate that the octamer is not 
transferred to competing DNA (at least hardly under physiological conditions, cf. 
It is also clear that repositioning should preferably not involve the simultaneous 
dissociation of all the 14 bindings sites which would be too costly to be induced by 
thermal fluctuations. A mechanism that requires only the breakage of a few contacts is 
loop diffusion. Strictly speaking one has to distinguish here between two different kinds 
of loops: small loops or bulges (2.4.2) and big loops (2.4.3) that show a qualitative 
different energetics and lead also to a different picture of the overall nucleosome 
repositioning dynamics. Another possibility could be the above mentioned corkscrew 
motion as suggested in Ref. [301 (and even before by van Holde and Yager |15Uj ). This 
could be facilitated through the twist diffusion of small defects that will only require to 
break one or two contacts at a time, a mechanism discussed in Section 2.4.4. 

2.4-2. Bulge diffusion In Ref. [73] myself, Widom, Bruinsma and Gelbart argued that 
the repositioning of nucleosomes without dissociation from the DNA chain that wraps 
them might be possible through the diffusional motion of small intra-nucleosomal loops. 
This biological process is analogous to the familiar physical situation of reptation of 
" stored length" in polymer chains. Thirty years ago de Gennes [IE] discussed the motion 
of a flexible chain trapped in a gel, modelled by a matrix of fixed point-like obstacles that 
cannot be crossed by the polymer. Fig. 11 depicts schematically the mechanism whereby 
diffusion of these "defects" of stored length AL gives rise to overall translation of the 
chain. Specifically, when the loop moves through the monomer at B, this monomer 
is displaced by a distance AL. de Gennes wrote down a conservation equation for 
this motion of defects along the trapped chain and calculated its overall mobility, 
and thereby, in particular, the molecular-weight dependence of the overall translation 
diffusion coefficient. In our present situation the reptation dynamics do not arise from 
obstacles due to a host matrix (as in a gel) or to other chains (as in a melt), but rather 
to loops associated with unsaturated adsorption of the DNA on the protein complex. 
Similar physics arise in the lateral displacements of a linear polymer adsorbed on a bulk 
solid surface. Sukhishvili et al. |15H I152j , for example, have measured the translational 
motion of adsorbed polyethylene glycol (PEG) on functionalized (hydrophobic) silica, 
specifically the dependence of its center-of-mass diffusion constant on molecular weight. 
They find an unusual scaling behavior - but one that can be accounted for by "slack 
between sticking points" , so that lateral motion of the polymer proceeds via a caterpillar- 
like diffusion of chain loops. In the case of intra-nucleosomal loops considered in this 
section, one is essentially in the limit of infinite molecular weight, because of the chain 
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stored length AL 




Figure 11. According to de Genncs and Edwards translational diffusion of trapped 
drains can be envisaged by a reptation mechanism. When a defect of stored length 
AL passes through a monomer at B it is moved by that amount to a new position B' . 

length being large compared to the bead (solid substrate) diameter. Furthermore, one 
deals with a lower dimensional problem, since the DNA chain is wrapped (absorbed) 
on a ID path rather than a 2D surface. But the basic features of loop formation and 
diffusion, and subsequent motion of the overall chain - in particular the exclusive role 
of equilibrium fluctuations in driving these processes - are the same in both cases. 

As shown in earlier studies of competitive protein binding to nucleosomal DNA 
123 EHl EH] , thermal fluctuations lead to lengths of the chain becoming unwrapped at 
the ends of its adsorbed portion. If some length of linker is pulled in before the chain 
re-adsorbs, then an intra-nucleosomal loop is formed - see Fig. 12(b). In Ref. [Z^ we 
calculated first the equilibrium shape and length distribution of these loops, in terms of 
the chain bending stiffness A, adsorption energy per unit length A, and protein aggregate 
size Rq (to be more precise: the radius of curvature of the DNA centerline). We then 
considered the diffusion of these loops from one end of the nucleosome to the other. 
Finally, treating this motion as the elementary step in the diffusion of the nucleosome 
itself along the wrapping chain, we were able to make estimates of the nucleosome 
repositioning rates as a function of A, A, Rq, and solvent viscosity rj. 

I present here (in more detail than in our letter [ZS]) the calculation of the 
equilibrium statistical mechanical probability associated with the formation of a small 
intra-nucleosomal loop (large loops are considered in the next section). When no loops 
are present, a nucleosome consists of a length / of DNA chain wrapped continuously 
around the octamer, see Fig. 12(a). In reality, as discussed in Section 2.1, the 
configuration of the adsorbed chain is a left-handed superhelix (of contour length I) 
spanning the full height of a cylinder. One can proceed, however, without making any 
explicit assumptions about the shapes of either the histone octamer or the wrapped 
DNA. 

Consider a fluctuation in which some length of the chain becomes unwrapped (this 
can only happen at the end of the adsorbed portion of chain) and simultaneously some 
length, say AL, of linker (i.e., previously unadsorbed chain) is "pulled in" before the 
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Figure 12. (a) Top view of the defect-free nucleosome looking down the superhelical 
axis, (b) Geometry of an intra-nucleosomal bulge. 



chain re-adsorbs. The fluctuation has then produced a loop of contour length 

L = L* + AL (59) 

where L* is the "exposed" length of nucleosome associated with the loop; see Fig. 12(b). 
Note that AL shows a strong preference for values that are multiples of 10 bp lengths, 
since this is the periodicity of the DNA helical pitch (binding sites are located where 
minor grooves face inwards to the octamer, cf. Section 2.1); other values require a 
twisting of the loop DNA which is energetically costly. For the energy associated with 
forming a loop of this kind, one can write [75] : 

AU _ Eelastic ( Ip .\ J, _IP f ds ( Ip \ 

The first term is the bending energy of the loop with 1/R{s) being the local curvature 
at distance s along its contour, cf. also Eq. 1211 The second term accounts for the length 
L* which has been adsorbed and bent with curvature I/-R0 before loop formation. The 
mathematical details of the functional minimization of AU, Eq. EHl are presented in 
Appendix B. There it is shown that the formation energy of an optimal small loop for 
given extra length AL is approximately given by"*" 



In Ref. |75| we made one further approximation, namely - when inserting 9* into Eq. IB 121 - we 
neglected the second and third term that nearly cancel when one inserts the typical parameters of the 
nucleosome. This led us to Eq. (2b) in that paper. Here I present the full expression, Ea.lCT 
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Using Eq. one finds that the energy to form a loop of minimal size of 34 A one 
helical pitch) is roughly 27 ksT (assuming i?o = 43 A, A = 0.176 A-\ Ip = 500 A); 
larger bulges are even more expensive. Here seems to arise a conceptional problem of 
the small loop mechanism: the formation energy of a small loop is of the same order as 
the complexation energy of the nucleosome itself (ca. 30k bT, cf. Section 2.1)! However, 
note two points: (i) The worm-like chain model is not very reliable for such strong 
curvatures and might overestimate the actual bending energy, (ii) We assumed in the 
calculation leading to Eq. IHI] that TR^/A ~ n'^/O*'^ > 1 (cf. Appendix B). From 
Eq. IB13I follows, however, that 6* ~ 1.75 and hence ^ 3. So one is actually in 

the crossover region between small and large "tensions" T. 

To account for this fact, the calculation can be redone using a refined estimation 
for A±, namely A+ = 1/A_ = ^2 + ^TRl/A which is a much simpler expression than 
Eq.Elbut still is a rather good approximation; especially it shows the right asymptotic 
behavior, Eq. |B7l This leads to 6* ^ 1.3 and AU ^ 20A;bT for AL = 10 bp. 

Now the probability distribution for the formation of loops of size AL is then simply 
given by the corresponding Boltzmann factor, normalized such that the maximal number 
of loops (expected for the unphysical case AU = 0) is the geometrically possible one, 
i.e. l/L*: 

ue, (AL) ^ l-e-^u/'^-^ (62) 

Even the value 20kBT calculated above shows that the spontaneous formation of a small 
loop is a very rare event. I will check in the following if it occurs sufficiently often so 
that it could account for the experimentally observed autonomous repositioning rates 
discussed in the previous section. In order to proceed here the dynamics of the loops 
needs to be considered and that of the resulting nucleosome repositioning. 

The key idea [ZSj here is that diffusion of the histone octamer along the DNA is 
achieved by formation and annihilation of loops. Let D denote the diffusion constant 
relevant to this motion of the ball along the chain, and let w be the rate at which loops 
are formed (by incorporation of linker length) D = wAL'^. These loops "disappear" 
due to their diffusion "off" the ball, at a rate that is proportional to the instantaneous 
number of loops, i.e., at a rate C^n. Accordingly, the overall rate of formation of loops 
is given by w — CAn, which must vanish at equilibrium, implying Ueg = w /Ca- Since 
this number is much smaller than unity, we are justified in assuming that only one loop 
at a time needs to be considered in treating the diffusion of intra-nucleosomal loops. It 
follows from the Boltzmann expression for n^q (see l62|l that w, the rate of loop formation, 
is given by 

w ~ Ca {l/L*) exp {-AU/keT) (63) 
and D by wAL"^. 

It remains only to evaluate Ca-, characterizing the rate of diffusion of loops "off" the 
ball. Let denote the diffusion constant associated with this motion [D^ characterizes 
the diffusion of loops through a wrapped ball, as opposed to the coefficient D that 
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describes diffusion of tlie hall along tlie cliain). Since tlie distance wfiicfi the loop must 
move to leave the ball is /, the wrapping length, one can write C^^ ^ / . From 
the Stokes-Einstein relation one has furthermore that = ksT/C, where C, — t]L* is 
the friction coefficient of the loop, with r] the effective solution viscosity. L*, as before, 
is the exposed length of the octamer associated with the loop, and hence provides the 
loop size relevant to its diffusion along the (ID!) nucleosome path of the chain. This 
hydrodynamic description is justified by the fact that loop diffusion requires unbinding 
of only a single sticking site, whose binding energy is of order of a few ksT. Combining 
all of the results from this and the preceding paragraph then gives 

with L* given by Ea. lHTSl and AU by Eq. EH 

Recalling 9* = 1.3 and AU = 20/^^^ for AL ~ 34 A and taking reasonable 
estimates for rj (a centipoise), -Ro(43 A) and / (500 A) we find that D is of order 
10~^^cm^/s. Hence typical repositioning times are in the order of an hour; furthermore 
there is a strong dependence on the temperature. A closer comparison of these 
theoretical estimates with the experiments (discussed in Section 2.4.1) will be given 
in the discussion (Section 2.4.5) after I have also presented the theories for large loops 
and twist diffusion. 



2.4-3. Large loop repositioning The perturbation calculation presented in Ref. [73] and 
reviewed in the previous section allows only to study small loops that store an amount 
of excess length of 10 or 20 bp. To describe also large loops a different approach is 
necessary as it was presented by Kulic and myself in Ref. [ZHj- In that paper we made 
use of the Euler-Kirchhoff theory for the static equilibrium of rods which allowed us to 
describe loops of any given excess lengths. The outcome of that paper changed our view 
of how repositioning via loop formation should work; besides the local repositioning 
based on bulge diffusion there should also occur a long-range hopping via large loops - 
at least for the case of very low nucleosome line densities as they are often encountered 
in in vitro experiments. 

In Ref. [Zni we started again from the Hamiltionian given by Eq. IHOl The section 
of the DNA constituting the loop has a contour length L and is parametrized by its arc 
length s ranging from —L/2 to L/2. L is the sum of the exposed length L* = 29*Rq [29*: 
opening angle) and the excess length AL, see Eq. EHl In order to compute the ground 
state for a trapped intranucleosomal loop the total energvlUUlhas to be minimized under 
two constraints: (z) The excess length AL is prescribed so that the following relation 
between the opening angle 9* and the total loop length L has to be fulfilled 

AL = L- 29* Ro = const. (65) 

(a) At the two ends s = ±L/2 the rod has to be tangential on an inscribed circle of 
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Figure 13. The KirchhofF analogue for the case of a planar pendulum and a planar 
semiflexible rod under tension. The inset displays how to construct a intranucleomal 
loop by inscribing a circular disc representing the octamer. 



given radius (representing the nucleosome): 

y {L/2) 



Ro = = const. (66) 

Here x (s) and y (s) are the Cartesian coordinates of the rod axis as a function of the 
arc-length parameter s (cf. Fig. 13). The absolute value in the second constraint needs 
to be introduced formally for dealing with crossed rod solutions (which are considered 
later on) and can be omitted for simple uncrossed loops. 

For the analytical description of the loop geometry it is convenient to introduce the 
angle 6 = 6 (s) between the DNA tangent and the Y-axis (cf. Fig. 13) that describes the 
DNA centerline (note that this is the same angle 6 as the one introduced in the previous 
section, cf. Fig. 12). Integrating the sine (cosine) of 9 over the arc length parameter s 
yields the X (Y) Cartesian coordinate of any point along the rod, and the derivative 6' 
gives the rod curvature R~^. Furthermore the nucleosome opening angle 6* is simply 
related to 6 at the boundary, namely 6* = 9 {L/2) for simple loops and 6* = n — 6 {L/2) 
for crossed loops (see below). 

The two constraints Eqs. ESI and EEl can be rewritten in terms of 6 and then be 
introduced into the minimization by two Lagrange multipliers Ti and T2. We then 
arrive at the following functional 

j^{e{s)}=Aj {efds-l^-ksTxU* 
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+ Ti[L-{AL + L*)]+T2 



L/2 

COS Ods — Rn sin 9* 



(67) 



Here the first term is the bending energy, the second accounts for the exposed length 
L* = 26*Ro and the third and forth term are the imposed length and tangency 
constraints. Eq. can be rearranged in the more familiar form 

'^^^ 'a {O'f + Ta cos e) ds + h.t. (68) 



where b.t. denotes boundary terms (depending on 6* (L/2) only) that obviously do not 
contribute to the first variation inside the relevant s interval. The integral in Eq. IHHl is 
analogous to the action integral of the plane pendulum with A [6')^ corresponding to 
the kinetic and — T2 cos 9 to the potential energy. The latter analogy is nothing else but 
Kirchhoff 's kinetic mapping between deformed rods and the spinning top that contains 
the present problem as a simple special case (cf . the paragraph after Eq. ESI for a brief 
discussion of the Kirchhoff analogy). 

Kirchhoff 's analogy provides one directly with explicit expressions for DNA shapes 
subjected to twist, bending and various geometric/topological constraints. Here, for the 
case of planar untwisted rods, also called the Euler elastica, where the corresponding 
"spinning top" reduces to the plane pendulum, the rod conformations are most generally 
given by 



cos 6* (s) = 1 — 2msn^ ^-^ | mj (69) 
This can be integrated to obtain the general planar rod shapes in Cartesian coordinates: 



X [s) = 2v^Acn y— I mj (70) 

y{s) =2AEi^j\m}j-s (71) 

with sn, cn(. | m) (and later below dn) denoting the Jacobi elliptic functions with the 
parameter m and E [u \ m) being the incomplete elliptic integral of the second kind 
in its "practical" form |153j . The two parameters m > and A > in Eqs. 1701 and 
ITT] characterize the shape and the scale of the solution, respectively. These solutions 
are up to trivial plane rotations, translations, reflections and shifting of the contour 
parameter s — s + sq the most general solutions to the Euler elastica. For different 
parameters m one obtains different rod shapes corresponding to different solutions of 
the plane pendulum motion |125j . The case m = describes a pendulum at rest which 
corresponds to a straight rod. For < m < 1 one has strictly oscillating pendulums 
corresponding to point symmetric rod shapes where the turning points of the pendulum 
have their counterparts in points of inflection of the rod. For m < 0.72 the rod is 
free of self intersections like the one depicted in Fig. 13. For m larger than 0.72 the 
rods show varying complexity with a multitude of self-intersections and for m = 1 one 
has the homoclinic pendulum orbit corresponding to a rod solution with only one self 
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Figure 14. Diagram of solutions of the Euler elastica providing an overview of the 
possible loop shapes as a function of the shape parameter m and the contact parameter 
(7. Loops of constant excess length AL = 10 bp are located on the dashed lines. The 
solid lines separate regions with different geometrical characteristics: simple loops (" 1" 
to "3"), crossed loops ("4" to "6") and more exotic shapes ("7" to "10"). 



intersection that becomes asymptotically straight for s ±oo. For even higher values 
of m, i.e., for m > 1 one has revolving pendulum orbits corresponding to rods with self- 
intersections lacking point symmetry. Finally, the limiting case m — oo corresponds to 
the circular rod shape. 

In order to describe a trapped loop one needs to use Eqs. [TOl and |7T] imposing the 
constraints inSl and inni For details of this calculation I refer the reader to Ref. |7(ij . 
There we present explicit solutions for the scaling parameter A, the opening angle 9* 
and the excess length AL as functions of the "contact parameter" a = L/2A and the 
shape parameter m, i.e., A = A (a, m), 6* = 6* {cr,m) and AL = AL {a,m). Inserting 
A (cr, m) and 6* (cr, m) into Eq. 120 leads to the final expression for the loop formation 
energy 

4A 
Ro 



AU (a,m) = — 
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(72) 



with ± = sign {2E {a \m) — a). 

Now the problem of finding the ground state loop for given excess length AL reduces 
to a two variable (a, m) minimization of Eq. [Z21 under the constraint AL (a, m) = AL. 
This final step was performed numerically ^| . An overview over the different solutions 
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Figure 15. Density plot of tlie total loop energy, Eq. 221 as a function of m and a 
(same parameter range as in Fig. 14). The white lines denote lines of contant excess 
length (in multiples of 10 bp). 

can be obtained by inspecting some loop geometries in the resulting (a, m) parameter 
plane. Both parameter values a and m vary between and oo, though the loops of 
practical importance are all found within the range < m < 1 and < a < 5. In 
Fig. 14 this relevant section of the parameter space is depicted together with a few 
example loops. The dashed lines indicate parameter values which lead to constant 
excess length AL = 10 x 3.4 nm (corresponding to 100 bp). On these lines are located 
the loop shapes " 1" to "7" that are examples of such 100 bp-loops. The whole parameter 
plane is subdivided into regions of structurally different solutions that are separated by 
solid lines. The large region starting at cr = contains exclusively simple loops (like 
"1", "2" and "3") without self- intersections and nucleosome penetration. Above that 
simple-loop region there is a region that contains loops with a single self-intersection; it 
includes the branch of 100 bp-loops with the example configurations "4", "5" and "6". 
To the right there are nonphysical cases where the loops penetrate the nucleosome, like 
example "10". There are also three other regions with single and double crossing points 
("7", "8", "9") where the loop can be founds on the "wrong" side of the nucleosome 
like in "7" and "8". 

In Ref. we determined the energy minimizing loop for each value of excess 
length AL. Fig. 15 shows a contour plot of the loop energies, Eq. [Z21 in the (cr, m) plane 
for the same range of parameters as in Fig. 14. The nucleosome parameters chosen in 
Fig. 15 are the same as in Ref. [IE], namely i?o = 40 A, A = 0.23 A~^, Ip = 500 A; these 
values are close (but not identical) to the one used in the previous section. Shown in that 
picture are also the corresponding lines of constant AL (with AL = 1, 2, .., 50 x 3.4 nm). 
As already observed in Fig. 14 there are, for any given AL, different branches of {a, m) 
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Figure 16. The energy of the optimal loops for given excess length AL. The kink in 
the curve reflects the switch from simple to crossed loops. The dashed line gives the 
free energy for entropic loops that are much longer than their persistence length. 



values corresponding to uncrossed, simply crossed and other, more exotic, structures. 
For short excess lengths one finds that the loops with the smallest formation energy Af/, 
Eq. E21 belong to the simple, uncrossed kind. For example, the optimal loop geometry 
for a lOObp-loop is structure "2" in Fig. 14 which is located at the point where the lines 
of the lOObp-loops in Fig. 15 encounter the total energy minimum. Interestingly, the 
optimal loop shape switches from the simple to the crossed type when an excess length 
of ~ 500 A is reached. 

The ground state energy as a function of the excess length AL is given in Fig. 16. 
Consider first the simple loops that are energetically preferable for AL < 50 nm. 
Inspecting Fig. 16 one finds the remarkable fact that the loop formation energy is 
non- monotonous in that range. First it increases sharply (namely as At/ oc AL^/^, 
cf. Ref. |701, in accordance with the small loop behavior Eq. WiT\ . Then at some critical 
excess length AL = AL^rit (which is approximately ALcru ~ 2.2 x 3.4 nm for the above 
given parameters) the loop energy reaches a maximum AU{ALcrit) ~ 26kBT. Beyond 
that the energy decreases with increasing AL. 

In order to explain this behavior one might naively argue as follows: For excess 
lengths shorter than the DNA persistence length it is energetically unfavorable to store 
additional length into the loop because it requires increasing deformation of the loop 
DNA. On the other hand, for loops longer than Ip the bending energy contribution 
becomes very small; to add more length should even decrease this energy since the loop 
can lower its curvature. However, the occurrence of the maximum of AU at ~ 22 bp 
excess DNA lengths, a value that is considerably smaller than the persistence length, is 
surprising at first sight. 

The explanation for this small value is given in Ref. [ZO! • There it is shown that the 
condition for the critical excess length ALcru is given by a simple geometric distinction 
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Figure 17. Subcritical (a) and supercritical loops show a different behavior when 
additional length is added. In the former case the energy goes up, for the latter it goes 
down. This allows to understand why the maximum of AJ7 is reached already for a 
very short excess length, cf. Fig. 16. 

between two loop shapes: the subcritical loop (Fig. 17(a)) that has none of its tangents 
parallel to the X axis (i.e. 6 (s) 7^ 7r/2 for all s) and the supercritical loop (Fig. 17(b)) 
that has two or more tangents parallel to the X axis where 9 (s) = 7r/2. Now it can 
be easily envisaged that adding some extra length dL to a subcritical loop increases 
its energy At/ (cf. |76ij) whereas for supercritical loops additional length decreases its 
energy. In the latter case one might just cut the loop at the two points Pl and Pr in 
Fig. 17(b) and introduce there the additional length (this operation does not change 
the energy) and then relax the shape by letting it evolve to the new equilibrium while 
keeping 9* constant. It was demonstrated in Ref. [TBj that this condition of the parallel 
tangents indeed leads to the above given small value of ALcrit- 

The ground state of loops switches from simple uncrossed loops to crossed loops 
when one reaches an excess length of ~ 500 A . Here, however, arises an additional 
complication: As can be seen by inspecting loops "4", "5" and "6" in Fig. 14 these 
structures contain a self-intersection at the crossing point. Therefore in principle a 
planar theory cannot capture the geometry of crossed loops. One might thus leave the 
plane and describe the self-contacts of the rod with corresponding point-forces in 3d as 
done by Coleman et al. |134j in a general theory of rod self- contacts. However, since the 
self-avoiding crossed loops stay close to a plane in all cases of practical interest (namely 
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loops for which the self-contact point is not too close to the nucleosome) , it is here 
sufficient to treat the self- interaction as a perturbation (cf. Ref. jTH] for details). One 
finds an additional small contribution for crossed loops due to the out-of-plane bending 
caused by the self-contact. In Fig. 16 this contribution (a few /c^T) has already been 
added; as can be seen from this figure crossed loops are still favored for sufficiently long 
excess lengths AL > ALcross (here ~ 6OOA). This can be rationalized by the fact that 
for long enough loops the adsorption energy (proportional to 6*) starts to dominate 
over the bending energy so that loops with smaller 6*, namely crossed ones, become 
favorable. 

Increasing the length even further one leaves the energy- dominated regime in which 
entropic effects can be neglected due to short loop length L < Ip. For larger lengths 
entropic effects become important and one enters the entropic loop regime (cf. the 
discussion of large- leafed rosettes, Eq. IHEl and Ref. |136j ). In the large loop hmit 
where the loop is longer than several Ip the chain looses its "orientational memory" 
exponentially and behaves as a random walk which starts from and returns to the 
same point. The entropic cost for gluing the ends of this random walk together is 
approximately given by 

AU - 3/2kBT \n{AL/lp) + Eq (73) 

Here Eq ^ 6.5 k^T denotes the bending plus adsorption energy contributions of the 
overcrossing DNA segments that enter and leave the nucleosome in the large-loop limit 
AL ^00. As already discussed after Eq. EHlthe free energy minimum occurs at the 
crossover between the elastic (AL < Ip) and entropic (AL ^ Ip) region where the 
decreasing elastic energy is overtaken by the increasing entropic contribution - as can 
also be seen in Fig. 16. 

The free energy, Eq. EHl leads to an algebraically decaying probability w [AL) for 
the jump lengths scaling as w oc exp (-AU/kpT) oc (AL)"^/^ In general, power law 
distributions of the form w oc (AL)~^ with 7 > 1 lead to superdiffusive behavior of 
the random walker (here the nucleosome displacement along the DNA). According to 
Levy's limit theorem the probability distribution of the random walker (more precisely, 
the distribution of the sums of independent random variable drawn out from the same 
probability distribution w oc (AL)~^) converges to a stable Levy distribution of index 
7 — 1 |154| I155[ I156j . This so-called Levy-flight |157j differs in many respects from 
the usual diffusion process as for short time intervals big jumps are still available with 
significant probability. Moreover, all moments (besides possibly the first few ones) 
diverge. For the present case 7 = 3/2 even the first moment does not exist. Note that 
the value 3/2 is based on the assumption of an ideal chain (no excluded volume); in 
general the excluded volume leads to self-avoiding-walk statistics with a slightly larger 
value of 7 around 2.2 |156j (cf. also Ref. |158j ). In that case one has a finite value of 
the first moment, i.e., a finite average jump length. 

We presented in Ref. jTH] some numerical estimates of the dynamics of the 
nucleosome repositioning on DNA fragments of different lengths. The basic idea 
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is that the transition rate w (AL) for a jump of length AL is proportional to 
C^exp (— Af/ (AL)) with Af/ being the loop formation energy. The Arrhenius constant 

involved in the loop formation has in principle to be determined experimentally. 
A theoretical estimate [73] was reported in the previous section, cf. Eq. IHSl where it 
was shown that Ca corresponds roughly to the inverse life time of the loop. Hence 
^ pT]R/kBT ^ 10"^ - 10-S. In Ref. [70] we considered two DNA lengths: 
(147 + 90)bp (short segment) and (147 + 300)bp (intermediate length). For the short 
piece the octamer repositioning occurs on the time scale of hours (in accordance with 
the previous section); on the intermediate segment the repositioning times is of the 
order of seconds to minutes. Important in both cases is where the nucleosome initially 
starts. If the start position is at an end of the DNA fragment then the nucleosome jumps 
preferentially to the other end since large jumps are energetically favored, cf. Fig. 16. 
This leads to a fast relaxation of the initial position. Smaller jumps take also place but 
less frequent; these jumps, however, lead on long time scales to an equal distribution 
(in accordance to Boltzmann's law) of the octamer along the DNA fragment. On the 
other hand, if the initial position is chosen in the middle of the DNA piece then the 
relaxation process is slower since a smaller loop is initially required; this first jump is 
then preferentially to an end position. 

As mentioned in Section 2.4.1 the repositioning is often followed via a gel 
electrophoresis (cf. Fig. 10). It is therefore helpful to ask how the resulting band 
structure evolves with time j76j . Let us start again with an end positioned nucleosome. 
In that just mentioned - the octamer initially shows mainly jumps back and 

forth between the two ends. Since the mobility is symmetric with respect to the middle 
position this leads to the interesting conclusion that in a standard gel electrophoresis 
these jumps would not be detected at all! Only slowly shorter jumps will allow the 
nucleosome to inhabit positions away from the ends. But this slower process cannot be 
distinguished easily from a short-ranged diffusion (away from the initial end), a process 
like the one discussed in the previous section. On the other hand, when starting from 
the middle position the preference to jump to the end positions will lead to an initial 
population gap in the band structure between the fast band (end positions) and the 
slow one (middle position). This gap would not occur for short-ranged repositioning. 
For more details on the expected band structures, the reader is referred to Ref. [7Kj . 

2.4-4- Twist diffusion Let me now discuss twist diffusion which might be another 
possible mechanism for nucleosome repositioning. I will give here some first theoretical 
quantitative estimates for this mechanism (without any free parameters) based on recent 
calculations by Kulic and myself; the full presentation will be given elsewhere [77]. If 
a 1 bp twist defect (one missing or one extra bp) forms through thermal activation at 
one end and manages to get through to the other end, this results in a 1 bp step of the 
nucleosome along the DNA and at the same time in a rotational motion by ~ 36°, i.e., 
the nucleosome performs a short fraction of a cork screw motion. 

The possibility of twist defects was demonstrated as soon as the high resolution 
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crystal structure of the core particle was resolved [TU]. In that study the core particles 
were prepared from a palindromic 146 bp DNA and core histones assuming that the 
resulting complex would show a perfect two-fold symmetry. However, it turned out that 
one base pair is localized directly on the dyad axis so that one half of the nucleosomal 
DNA is of length 73 bp whereas the other is only 72 bp long. The missing base pair of 
the shorter half is, however, not localized at its terminus but instead at a 10 bp stretch 
close to the dyad axis (cf. Fig. 4d in ^U] that shows a superimposition of the two DNA 
halves). The reason is presumably the attraction between the DNA termini of adjacent 
particles in the crystal (cf. Fig. 4c in that paper) that try to come close to mimic a 
base pair step at the cost of forming a twist defect far inside the wrapped chain portion. 
In fact, crystals of core particles with 147bp DNA do not show this defect |llj . 

To proceed further we describe the DNA chain within a Frenkel-Kontorova model, 
i.e., we view it as a chain of particles connected by harmonic springs in a spatially 
periodic potential. The original Frenkel-Kontorova model was introduced more than 
sixty years ago in order to describe the motion of a dislocation in a crystal |159j . In 
the meantime variants of this model were applied to many different problems including 
charge density waves |16()j , sliding friction |16H I162j , ionic conductors |163| I164j , chains 
of coupled Josephson junctions |165j and adsorbed atomic monolayers |166| I167j . Here, 
in the context of DNA adsorbed on the protein octamer, the beads represent the 
basepairs. The springs connecting them have an equilibrium distance of 6 = 0.6 nm 
(which is here taken to be the distance along the DNA sugar-phosphate backbone, not 
the distance h = 0.34 nm along the fiber axis) and the spring constant is chosen such 
that it reflect the DNA elasticity. Specifically 



with Xn being the position of the the nth basepair measured along the helical 
backbone and Kb /2 ~ 70 — IOOAj^T accounts for the coupled twist-stretch elasticity 
|168| I169t I170j . Finally, the external potential comes from the contact points to the 
octamer. The distance between neighboring contact points is 10 bp which corresponds 
to 60 nm along the arclength of the minor groove. A contact point at position xq is here 
modelled by the following function 



with ^ (x) = 1 for X > and zero otherwise. The two parameters, the depth Uq of 
the potential and its width a, can be estimated as follows. Uq represents the pure 
adsorption energy per point contact which can be estimated from competitive protein 
binding j2Zll2H] to be of order GksT (cf. Section 2.1). The other parameter, a, can then 
be estimated from the fluctuations of the DNA in the crystal (measured by the B-factor, 
cf. Fig. lb in JHl)- The fluctuations of the DNA at the binding sites are much smaller 
than in the middle in between. Using a quadratic expansion of Eq. ESI one finds from 




(74) 




(75) 
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a straightforward normal mode analysis j77] that a ^ 6/2, i.e. the adsorption regions 
lead to a strong localization of the DNA. 

Now having all the numbers at hand we can answer the question whether the twist 
defects are localized between two contact points. The deformation energy of the defect 
localized along a lObp stretch is of order IksT. On the other hand, by distributing the 
defect homogenously over 20bp the elastic energy goes down by ~ 7/2/cbT at the cost 
of releasing the adsorption point in the middle (roughly 6k bT). The smearing out of 
the defect costs therefore ~ 3kBT. This already shows that the kink is not so strongly 
localized which points towards a high mobility of the twist defect. 

Let us assume now that a kink with one missing basepair is located between two 
binding sites. When this kink jumps to the neighboring lObp (say to the right) it has 
to cross a barrier of height ~ 3k bT. Using Eqs. [^andEHlwe obtain an explicit form 
for the barrier and together with the single bead friction ( ^ lQ~'^kBTs/nm? \\.7\\ 1172] 
are able to calculate the Kramer's escape rate for the kink from the given localized 10 
bp stretch to a neighboring one. This leads us to a typical time tgtep ~ 10ns for going 
from one stretch to the next one (cf. Ref. ^| for details). 

To determine the rate at which twist defects are formed at the entry/exit points 
of the DNA one can now use an argument similar to the one presented in Ref. [TSj (cf. 
also Section 2.4.2): the ratio of the life time tufe of a kink to the time interval 
between two kink injection events at the end of the wrapped DNA portion equals the 
probability to find a defect on the nucleosome, i.e. tufe/Unj — NgUee'^^^''^^ ~ 10~^. 
Here NsUe = 13 denotes the number of possible positions of the defect between the 14 
binding sites. 

How is the average life time tuje of a defect related to tstep, the typical time needed 
for one step? It is possible to calculate the mean first passage time for a defect that 
starts at one end (say the left one) and leaves the nucleosome at the same end, r^g/t, or 
at the other end, Tright- From Ref. |173j one finds r^g/t = (25/6) tgtep and Tright = '^Ststep- 
Furthermore, the probability to leave at the left end is p;e/t = 12/13 and at the right 
end Pright = 1/13 |17Hj which gives the life time as the weighted average tufe = 6t step- 
Only a fraction Pright of the defects reaches the other end and will lead to a repositioning 
step, i.e., the time of one diffusion step of the nucleosome along the DNA is of order 
T = tinj/pright ~ 10~'^s whcrc usc was made of the above presented relations between the 
time scales. From this follows then directly the diffusion constant D of the nucleosome 
along the DNA: D = b"^/ (2T) ^ 7 x lO^^^cm^/s (with b = 0.34nm). 

Therefore we find a diffusion constant that is much larger than the one expected 
for the repositioning via bulges {D ^ lO^^^cm^/s, cf. Eq. l64j) . Most importantly, it is 
also orders of magnitude larger than the diffusion constant observed in the experiments. 
How can this apparent inconsistency be resolved? 

Most likely, the diffusion is considerably slowed down due to the quenched disorder 
stored in the base pair sequence of the DNA. In fact, the bulk of the repositioning 
experiments has been made on DNA with rather strong positioning sequences leading 
to a strong rotational positioning of the nucleosome. Starting from a preferred position 
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the nucleosome would arrive already after 5 steps to the left (or to the right) on the 
"wrong" site of the nucleosome - forcing the DNA to be bent in an unfavorable direction. 
This means that the nucleosome needs to cross a barrier in order to reach a position 10 bp 
apart. For instance, in the case of the 5S rDNA sequence theoretical estimates indicate 
a barrier height of the order IQksT \\7A\ 117,^ . This leads to a strong reduction of the 
effective diffusion constant, namely Df.ff k, De~^^ ~ 10~^^cm^/s - a value comparable 
to the one found for bulge diffusion. In Ref. [77] we formally include these effects into 
the Frenkel-Kontorova framework by introducing an octamer-fixed bending field and by 
attributing "bending charges" to the beads. This allows us to give a rough quantitative 
treatment of the nucleosome mobility as a function of the underlying base pair sequence. 

2.4-5. Discussion: Bulge versus twist diffusion Comparing the different repositioning 
mechanisms presented in the previous sections one has to conclude that they lead to 
very different "sliding" scenarios. On short DNA fragments repositioning could in 
principle work via bulge (Section 2.4.2) or twist diffusion (Section 2.4.4); large loops 
cannot occur because there is not enough free DNA length available. Bulge diffusion is 
rather slow (timescale of hours) since the formation of a small loop is costly - mainly 
because the opening angle of a bulge is rather large leading to several open binding 
sites. The repositioning rates should show a strong temperature dependence as well as 
a strong dependence on the adsorption strength (i.e., a strong dependence on the ionic 
conditions). The preferred repositioning steps are multiples of ten base pairs. 

On the other hand repositioning via twist defects should be much faster (timescale 
of seconds). The nucleosome should slide in a cork screw motion along the DNA and 
should forget its initial position rather quickly. However, if the underlying DNA sequence 
induces a strong rotational positioning signal the timescale becomes comparable to that 
of small loop repositioning. Even more, due to the underlying base pair sequence one 
should expect a 10 bp spacing between the dominant positions that is, however, here not 
the result of 10 bp jumps but just reflects the relative Boltzmann weights of favorable and 
unfavorable positions. The estimates of the diffusion constants of these two mechanisms 
are too unreliable (activation energies appear in the exponent!) to allow one to predict 
which of the mechanisms should be favored. If on the other hand a rather homogenous 
DNA sequence is used, our prediction is that cork screw motion is the much faster and 
therefore predominant mechanism. The experiment by Flaus and Richmond |148j goes 
already in this direction; comparing their experimental results (cf. Section 2.4.1) with 
the theoretical pictures seems to point towards sliding motion for one of the positioning 
sequences (the one that has a homonucleotide tract) whereas it is not clear whether the 
nucleosome escapes from the rotational positioning trap via bulge or via twist diffusion. 

On long DNA fragments single nucleosomes could also be repositioned via large 
loops (Section 2.4.3). Our theoretical model suggests that large loop repositioning 
would be much faster than bulge diffusion. Also it should be expected that a similar 
mechanism allows the nucleosome to be transferred to competing naked DNA chains. 
As discussed in the experimental section above there is, however, not much evidence 
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for such processes. Only Watkins and Smerdon jl39j report such a nucleosome transfer 
to free DNA at higher ionic strength. This again allows to speculate that each of the 
different mechanisms might play a dominant role in a certain parameter range. To come 
to more definite conclusions more systematic experiments have to be made on short as 
well as long DNA molecules with and without positioning sequences under varying ionic 
conditions. 

Repositioning in vivo might be actively facilitated by chromatin remodelling 
complexes whose action is currently studied in vitro (reviewed in Refs. [3 1221 ESI OH ESI)- 
There are two major families: SWI/SNF and ISWI. They both burn ATP to enhance 
nucleosome dynamics but their underlying modes of action seem to be fundamentally 
different. The SWI/SNF class disrupts many of the DNA- nucleosome contacts making 
the nucleosomal DNA vulnerable to DNA digestion. It has been even observed that some 
of these complexes are capable to transfer the octamer to another DNA chain (176] . This 
might indicate that their mode of action is the creation of large loops (similar to the one 
discussed in Section 2.4.3) that could lead to large repositioning steps. In fact, Bazett- 
Jones et al. |177j observed that the SWI/SNF complex creates loops on naked DNA 
as well as on bead-on-string nucleosome fibers (cf. the electron spectroscopic images, 
Fig. 1 and 3, in that paper). On the other hand, the mode of action of the ISWI family 
seems not to interrupt the nucleosome-DNA contact on an appreciable level. Since these 
complexes induce nucleosome repositioning it has been speculated that they might work 
via twist or bulge diffusion. In the meantime the latter mechanism seems to be more 
likely since ISWI induced nucleosome sliding appears even if the DNA is nicked and 
hence a torsion cannot be transmitted between the complex and the nucleosome to be 
shifted HIH]. 

Another interesting and very prominent system known to mediate nucleosome 
repositioning is unexpectedly the RNA polymerase. It is found to be able to 
transcribe DNA through nucleosomes without disrupting their structure, yet moving 
them upstream the DNA template, i.e., in the opposite direction of transcription 
[Sni EH EHl Eni EOj ■ To rationalize this seemingly paradoxical finding Felsenfeld et al. [10] 
introduced a model which assumes that the polymerase crosses the nucleosome in a loop. 
This would indeed explain the backwards directionality of repositioning. Note that such 
a loop would have a different shape than the ones discussed above since polymerases 
induce a kink at the DNA with a preferential angle of ~ 100° |1791ll80lll81j . This means, 
however, that an RNA polymerase sitting in an intranucleosomal loop would soon get 
stuck since it transcribes the DNA in a cork screw fashion; this would complicate this 
mechanism [10] • It might well be that this effect only occurs on short DNA fragments 
as used in the experiments. If so, it would be an artefact that would not work in 
vivo. In that case another mechanism, namely induced cork screw motion of the bound 
DNA towards the polymerase and subsequently the recapturing of the nucleosome at its 
exposed binding sites by the other end of the DNA fragment might also be a possible 
scenario |182j . 
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3. 30-nm fiber 

3.1. Solenoid versus crossed- linker model 

Whereas the structure of the core particle has been resolved up to atomic resolution (TU], 
there is still considerable controversy about the nature of the higher-order structures 
to which they give rise. When stretched, the string of DNA/histone complexes has 
the appearance of "beads-on-a-string". This basic structure can be seen clearly when 
chromatin is exposed to very low salt concentrations, and is sometimes referred to as the 
"10-nm fiber" When the ionic strength is increased towards physiological values 

(100 mM), the fiber appears to thicken, attaining a diameter of 30 nm ^31- Linker 
histones (HI or H5) play an important role in this compaction mechanism: In their 
absence fibers form more open structures [T2j. These strongly cationic proteins act 
close to the entry-exit point of the DNA. They carry an overall positive charge and 
seem to bind the two strands together leading to a stem structure [l^; in fact, this stem 
is missing in the absence of linker histones. 

There is a longstanding controversial discussion concerning the structure of the 30- 
nm fiber ^1^1^21 EI- There are mainly two competing classes of models: the solenoid 
models [Hlllllin]; and the zig-zag or crossed-linker models [201 IHl 122 123 121 1113 ■ In 
the solenoid model (depicted in Fig. 18(a)) it is assumed that the chain of nucleosomes 
forms a helical structure with the axis of the core particles being perpendicular to the 
solenoid axis. The DNA entry-exit side faces inward towards the axis of the solenoid. 
The linker DNA (shown as a thick lines at the bottom of Fig. 18(a)) is required to 
bent in order to connect neighboring nucleosomes in the solenoid which in turn requires 
strong nucleosome-nucleosome interactions to hold this structure together. The other 
class of models posits straight linkers that connect nucleosomes located on opposite sides 
of the fiber. This results in a three-dimensional zig-zag-like pattern of the linker (cf. 
Fig. 18(b)). 

Images obtained by electron cryomicroscopy |23] should in principle be able to 
distinguish between the structural features predicted by the two different models. The 
micrographs show a zig-zag motif at lower salt concentrations and they indicate that 
the chromatin fiber becomes more and more compact when the ionic strength is raised 
towards the physiological value. A similar picture emerges also from atomic force 
microscopy I184j . However, neither method allows to identify the linker geometry 
at physiological ionic conditions so that one still cannot exclude the possibility that the 
fiber folds close to physiological conditions into a solenoid-like structure by a bending 
of its linkers. This is in fact the structure that is depicted in most of the standard 
textbooks on cell biology (e.g. Ref. PP). Also X-ray diffraction data that constituted 
the basis for many models lead to controversial interpretations, cf. Ref. ^Hj for a critical 
discussion. 

In view of this fact it is an important recent experimental achievement that single 
chromatin fibers can be stretched via micromanipulation techniques |HI1 IH21 EBj ■ The 
force-extension curves allow in principle to discern between the different structures. So 
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Figure 18. The two competing models for tlie 30-nm fiber: (a) tlie solenoid model 
and (b) the crossed linker model. Both are shown from the side and from the top; the 
latter view allows to distinguish the different linker geometries. Note that these are 
idealized models; real fibers are beheved to be less regular [T?)j . 



far computer simulations (HHl as well as analytical approaches |24t IM] to chromatin 
fiber stretching seem - when comparing their predictions to the experimental data - to 
support the crossed linker models. 

Another intriguing way that might allow to discriminate experimentally between 
these two types of structures is to measure the fiber orientation in strong magnetic fields 
- as it has been done already long ago |185j . Such a method has been used successfully 
to determine the persistence length of naked DNA |186| I187j . One makes use of the 
anisotropic magnetic susceptibility of the base pairs that cause the DNA double helix to 
orient its axis perpendicular to the field. As a consequence, single core particles orient 
their DNA superhelix axis parallel to the field |185j . Since the nucleosome axes in the 
two fiber models are oriented in different directions with respect to the fiber axis (cf. 
Fig. 18), an external field would induce orientations of the two fiber models in different 
directions. 

In the following I will focus on analytical models and computer modeling of the 
chromatin fiber which all belong to the class of crossed-linker models. In the next 
section the possible geometrical structures that follow from regular two-angle fibers (a 
"generalization" of crossed-linker models) are presented. Section 3.3 gives a speculation 
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Figure 19. Fraction of a two-angle fiber containing four nucleosomes. The two 
angles, the defection angle 9 and the dihedral angle are depicted together with 
the nucleosome diameter 2i?o and the "linker length" B. The arrows denote the 
nucleosomal superhelix axis, cf. Fig. 2. 

about the "optimal" fiber design from a biological point of view. Then in Section 
3.4 I will give a detailed account on the mechanical properties of the fiber comparing 
analytical results and computer models to recent stretching experiments. Finally, in 
Section 3.5 I report on a recent model that relates the degree of fiber swelling to the 
ionic strength. 

3.2. Structure diagram of the two- angle fiber 

To address the folding problem of DNA at the level of the 30-nm fiber myself, Gelbart 
and Bruinsma introduced a mathematical description for the different possible 
folding pathways which was based on Woodcock's crossed-linker model [20] (cf. also 
a related study on closed minichromosomes {188. ). At the simplest level, we assumed 
that the geometric structure of the 30-nm fiber can be obtained from the intrinsic, single- 
nucleosome structure. The specific roles of linker elastic energy, nucleosome-nucleosome 
interaction, preferred binding sites, HI involvement, etc. was then treated afterwards as 
"corrections" to this basic model [211 E]. To see how single-nucleosome properties can 
control the fiber geometry, consider the fact that DNA is wrapped a non-integral number 
of turns around the nucleosome, e.g., l-and-3/4 times (147 bp) in the case of no HI. 
This implies that the incoming and outgoing linker chains make an angle 9 with respect 
to each other - the entry-exit angle vr — ^ is nonzero. In the presence of the histone HI 
(or H5) the in- and outcoming DNA are glued together along a short section resulting 
in a stem-like structure [23] • While the precise value of the resulting exit-angle depends 
on salt concentration, presence or absence of linker histones, degree of acetylation of the 
histones, etc. (discussed in Section 3.4) one may nevertheless assume 6* to be a quantity 
that is determined purely at the single-nucleosome level. 

Next, there is a rotational (dihedral) angle between the axis of neighboring 
histone octamers along the necklace (see Fig. 19). Because nucleosomes are rotationally 
positioned along the DNA, i.e., adsorption of DNA always begins with the minor groove 
turned in towards the first histone binding site, the angle is a periodic function of the 
linker length 5, with the 10 bp repeat length of the helical twist of DNA as the period. 
There is experimental evidence that the linker length shows a preferential quantization 
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Figure 20. Diagram of geometries of two-angle fibers in the (6', (/i)-plane. Sliown 
are some example configurations with the arrows denoting their position in the plane. 
The lines give the boundaries to the forbidden structures due to short-range excluded 
volume (large 9 values) and long-range excluded volume (small 0- values). 

involving a set of values that are related by integral multiples of this helical twist |189j . 
i.e., there is a preferred value of 0. 

Treating the pair of angles {6,4'), together with the linker length B, as given 
physical properties (even though in vivo they are likely under biochemical control), the 
geometrical structure of the necklace is determined entirely by 6, and B. The model 
only describes linker geometry and does not account for excluded volume effects and 
other forms of nucleosome-nucleosome interaction; it assumes that the core particles are 
pointlike {Rq = 0) and that they are located at the joints of the linkers. The model also 
assumes that the linkers are straight. The {6, (f))-model is similar to the freely rotating 
chain model encountered in polymer physics literature (see, for instance, Ref. [EDI)- The 
main difference is that in the present case there is no free rotation around the linker 
and so torsion is transmitted (see also Ref. |19Uj ). 

Before giving a detailed discussion of two-angle geometries let me provide a short 
overview of the possible structures in the {6, (f))-space that is shown in Fig. 20. Both 
angles 6 and can each vary over the range to vr. At the edges of the diagram where one 
of the angles assumes an extremal value, the configurations are always planar. On the 
line = are located circles (see structure "2" in Fig. 20) and star-type polygons (that 
are closed for specific values of 6 like "5"). The planar zig-zag-structures are found on 
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the line (p = n (cf. "6" and "7"); for 6 = one has straight configurations ("1") and for 
(j) = TT "dimer" structures ("8"). If one moves from the fine = towards larger values 
of the circles and star-like polygons stretch out into the direction perpendicular to 
their plane, forming solenoids ("9") and fibers with crossed linkers ("11"), respectively. 
On the other hand, if one starts at the top of the diagram (0 = vr) and decreases the 
value of the planar zig-zag structure extends into the third dimension by becoming 
twisted ("11"). Various examples of two-angle fibers were displayed by Woodcock et 
al. |2ni in their Fig. 2, namely fibers with 6 = 150° and many different values of 0, 
corresponding to a vertical trajectory on the right-hand side of Fig. 20. Three different 
configurations with a fixed value of and different values of 6 are displayed in Fig. 3(c) 
in another paper by these authors [Z3j . 

An analytical decription of the structures can be achieved as follows (cf. Ref. |24j): 
It is possible to construct a spiral of radius R and pitch angle 7 such that the 
nucleosomes - but not necessarily the linker chain - are located on this spiral. The 
nucleosomes are placed along the spiral in such a way that successive nucleosomes have 
a fixed (Euclidean) distance B from one another. From straightforward geometrical 
considerations we derived in Ref. [23 analytical expressions that relate pitch angle 7 and 
radius R of the solenoid as well as Sq (defined as the vertical distance between successive 
"nucleosomes" along the helical axis) to the pair of angles 6, and linker length B. 
The corresponding relations B = B{'y,R,so), 6 = 6^(7,/?, sq) and = 0(7,i?, sq) 
are Eqs. (32) to (34) in Ref. [21]. I present here the reverse relations that have the 
advantage that they allow the direct calculation of the overall fiber geometry from the 
local geometry. Specifically, the radius R of the master solenoid is given by 
R = B sin {6/2) 

2-2cos2(^/2)cos2(0/2) ^ ' 

and its pitch angle ip by 

, tan {6/2) arccos (2 cos^ (6/2) cos^ (0/2) - 1) 

cot -0 = (77) 

2 sin (0/2) - cos2 {6/2) cos^ (0/2) 

Finally, the distance sq of neighboring nucleosomes along the fiber axis is obtained from 

B sin (0/2) 

So = . (78) 
y'sec^ {6/2) - cos2 (0/2) 

Using these relations, it is straightforward to construct a catalog of structures. 

If either one of the angles 6* or assumes the value or vr, then the resulting 
structure is planar. Consider first the line = 0. If one also has 6 = the fiber forms 
a straight line ("1" in Fig. 20). For small non- vanishing 6 the structure forms a circle 
of radius i? ~ 5/6' (as follows directly from Eq. [7B|) . For the special case 6 = 27r/n, 
with n an integer, the ring contains n monomers before it repeats itself and one obtains 
a regular polygon ("2"). The special case 6 = 'n/2 corresponds to the square ("3"). 
With increasing 6 the radius of the circle shrinks and approaches asymptotically the 
value B/2. For 6 = -n {n — \) /n with n being an odd integer one encounters a series of 
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closed star-like polygons with n tips. In particular, n = 3 corresponds to the equilateral 
triangle ("4") and n = 5 to the regular pentagram ("5"). 

Next consider the case = vr and 9 arbitrary. This case corresponds to 2D zig- 
zag-like structures, as shown by "6" and "7" at the top of Fig. 20. The length of a 
fiber consisting of N monomers is given by L = sqN = B cos [6/2) N (cf. Eq. I7H|) 
and the diameter is given by 2R = B sin [6/2) (cf. Eq. [7H|) . Note that the length of 
the fiber increases with decreasing 6. Finally, there are two remaining cases of planar 
structures: ^ = with an arbitrary value of leads to the straight line mentioned earlier 
("1"); 6 = 71 and arbitrary corresponds to linkers that go back and forth between two 
positions ("8"). 

Structures with 6^0 and 7^ form three-dimensional fibers. For small angles, 
9 <^ 1 and -C 1, structures resemble solenoids (see "9") where the linkers themselves 
follow closely a helical path corresponding to that of the master solenoid. For these 
structures one finds from Eqs. [70 and [THl the following limiting behavior of the fiber 
radius and length (for N monomers): 

BO BNcf) 

Furthermore, the pitch angle ip is given by 



R ~ L ~ f79l 

^ - ^2 I fl2 ' ^ - rur-r-m \ ' ^> 



^-7^ (81) 



cotV^~- (80) 

This suggests a classification of solenoids into dense helices with small pitch angle 
■0 ~ (j)/9 for 0^6' and open helices with large pitch angle ~ 7^/2 — O/cf) for cf) ^ 9. 
Other geometrical information can be obtained easily. For instance, the vertical distance 
d between two turns follows from d = 27cR/ cotip to be 

27r05 

Dense helices, -C ^, are characterized hj d <^ R and open ones by d ^ R. 

Structures where is still small but where the entry-exit angle 9 is large, i.e. 
n — 9 <^ TT, form fibers with crossed linkers. As discussed above for = one encounters 
star-shape polygons that are closed for 9 = tt {n — 1) /n with n odd. For non- vanishing 
<^ 1 the star-shaped polygons open up in an accordion- like manner into a three- 
dimensional fiber with the following radius and length (for monomers): 

R^ ^ . Lc^—^cot{9 2) 82 

2sm(6'/2) 2 

Assume now that 9n = Tf {n — 1) / n so that the projection of the fiber is a closed polygon 
(this is only strictly true for = but it is still a good approximation for 0^1). For 
this set of angles monomer i and i -\-n come very close in space; their distance d follows 
from the master solenoid that has n — 1 turns in between these two monomers: 



^^2n{n-l)R^n^ (83) 
cotV ' 
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Finally structures with a rotational angle close to vr, say = vr — 5 with 5^1, 
lead to twisted zig-zag structures - see " 11" . In this case monomer i + 1 is located nearly 
opposite to the iih monomer, but slightly twisted by an angle 5. Monomer i + 2 is then 
on the same side as monomer i but slightly twisted by an angle 25 and so on. The 
geometrical properties of the resulting fiber are the following 

— sin(0/2), Lc^ BN cos{9/2) (84) 

and show only a higher order dependence on that we gave explicitly in Ref. [23- For 
= 7r one recovers the planar zig-zag structure for which Eq. |M] becomes exact. 

If one takes into account the excluded volume of the core particles, then certain 
areas in that phase diagram are forbidden - reminiscent of the familiar Ramachandran 
plots used in the study of protein folding |191j . For simplicity we assume in the following 
that the core particles are spherical with a radius Rq and that their centers are located 
at the joints of two linkers, cf. Fig. 19. There are two different types of interactions. 
One is between monomers at position i and z it 2 (short range interaction), and leads to 
the requirement that the entry-angle must be sufficiently small: 

e <2 arccos {R^/B) (85) 

This condition excludes a vertical strip at the right side of the diagram, as indicated in 
Fig. 20 by a dashed line. 

There is also a long-range excluded volume interaction that comes into play for 
small values of 0. This is apparent for the case = where one finds planar structures 
that run into themselves. Starting with a circular structure one has to increase above 
some critical value so that the pitch angle of the resulting solenoid is large enough so 
that neighboring turns do not interact. This leads to the requirement d > 2Rq with d 
given by Eq. IHU (using (f) <^ 9), i.e., 

4->-^ (86) 

TT B 

For the large 9-case (fibers with crossed linkers) one finds from Eq. IHSl the condition 

^>l^ (87) 

The two conditions, Eqs. IHEl and |HZ1 shown schematically as a dotted curve in Fig. 20, 
lead to a forbidden strip in the structure diagram for small values of 0. 

Figure 20 does not show the interesting "fine-structure" of the boundary of that 
forbidden strip that is due to commensurate-incommensurate effects. I already noted 
that there are special 6'-values for which the projection of the linkers forms a regular 
polygonal star {9^ = tt (n — 1) /n) or a regular polygon {9'^ = 271 /n) (for small values 
of 0). In these cases the nucleosomes i and i + n "sit" on top of each other. On the 
other hand, for other values of 9, monomers of neighboring loops will be displaced with 
respect to each other. In this case monomers of one loop might be able to fill in gaps 
of neighboring loops so that the minimal allowed value of is smaller than estimated 
above. We are currently exploring the interesting mathematical problem of the exact 



CONTENTS 



66 



boundary line that is also sensitive to the exact nucleosome shape |192j . The dotted 
line in Fig. 20 only represents the upper envelope of the actual curve. 

The above given discussion of the two-angle model was based of the assumption of 
a perfectly homogeneous fiber where B, 9 and are constant throughout the fiber. For 
a discussion of the effect of randomness in these values on the fiber geometry I refer the 
reader to Ref . [211 ■ 

3.3. Chromation fiber: optimization of design? 

If one assumes that the chromatin fiber has a relatively regular structure and that the 
linker DNA is straight, then the two-angle model might be a good description of the fiber 
geometry. In that case the question arises where in the structure diagram, Fig. 20 the 
30-nm fiber is actually located. The diagram, by itself, does not favor any structure over 
another. However, the diagram plus the formulae given above allowed use in Ref. [3] 
to invoke the following two criteria to optimize the structure of the 30-nm fiber and to 
check a posteriori their usefulness. The two suggested criteria are: 

(i) maximum compaction 
(a) maximum accessibility 

The first criterion is obvious: inactive chromatin should be packed as dense as possible 
because of the very large ratio of DNA length to nucleus size (cf. also Ref. |193j to see 
how severe this packing problem actually is). By the second criterion we meant that a 
local accessibility mechanism is required for gene transcription and that this mechanism 
should somehow be optimized (see below). 

In order to attain maximum compaction one needs structures that lead to high 
bulk densities p = 1/ (2y/3R'^so^ (assuming that the 30-nm fibers are packed in parallel 
forming a hexagonal lattice). A comparison of the 3d densities of all possible structures 
shows that fibers with internal linkers have highest densities p, namely (cf. Eq. l82|l 

^ ~ 2^30 (tt - 6) 53 ^ ' 

In particular, the highest density is achieved for the largest possible value of 6 and 
the smallest possible value of that is still in accordance with the excluded volume 
condition. This set of angles is located at the point where the dotted curve and the 
dashed line in Fig. 20 cross each other. Apparently this also represents the only region 
in the phase diagram where excluded volume effects are operative on a short-range and a 
long-range scale at the same time, i.e., nucleosome i is in close contact with nucleosome 
i — 2 and i + 2 as well as with nucleosomes father apart along the contour length of the 
necklace. This unique set of angles is given by 9raa.x = 2arccos (Rq/B), cf. Eq. ISil and 
0^i, ~ (8/7r) (i?o/5), cf. Eq.EZI 

In order to achieve maximum accessibility we looked in Ref. [25 for structures 
that, for a given entry-exit angle vr — ^ of a highly compacted structure, achieve the 
maximum reduction in nucleosome line density pl = Sq^ for a given small change 
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of the angle 9. In other words, we looked for a maximum of (ipj^/c?^ which we called the 
"accessibility". Interestingly, the accessibility is maximized at the same unique pair of 
angles (^max, 0min)- This can be seen from its angle dependence for fibers with crossed 
linkers 

^ ^ 1 (89) 

de (j){n-efB ^ ^ 

Note that this change in with 6 is achieved by changing the number of monomers per 
vertical repeat length d. The length d itself is only weakly dependent on n according to 
Eq.lHl 

The above given formulas are now compared with experimental results. For chicken 
erythrocyte chromatin one has B ^ 20 nm (center-to-center distance of nucleosomes, 
[Ej). Together with i?o ~ 5 nm this leads to 6'max — 151°, 0min — 36° and pL — 6.9 
nucleosomes per 11 nm (using Eqs. [THllHHlandlHZl the approximate formula. Eg. IH^ gives 
Pl ~ 6.8). The theoretically derived values can now be compared with the experimental 
ones reported by Bednar et al. [22 for chicken erythrocyte chromatin fibers. From their 
table 1 one finds that for an ionic strength of 80 mM (which is close to the physiological 
value) 6 ~ 145° and pi = 6.0 nucleosomes per 11 nm. Furthermore, electron 
cryotomography-constructed stereo pair images of an oligonucleosome (cf. Fig. 3(b) 
in Ref. |2S1) indicate that the chromatin fiber might indeed have the structure of a fiber 
with crossed linkers, with n ^ 5; this would correspond to 6 = 7i {n — 1) /n ^ 144°. 

Information concerning the preferred value for (p can be obtained from the measured 
statistical distribution of the nucleosome repeat lengths. This distribution shows 
statistically preferred linker lengths equal to 10k + Ibp with k a positive integer jl89ij . 
which, in turn, indicates that the rotation angle corresponds to a change in helical 
pitch associated with 1 bp, i.e. 360°/10 = 36°. This value coincides with 0min, the 
value that we estimated for maximum compaction. However, the statistical uncertainty 
around the expectation values for the nucleosome repeat length is sufficiently large to 
make this estimate for less reliable. 

The second feature, the local accessibility, can be monitored in vitro by changing the 
salt concentration. Bednar et al. report, for example, that 6 decreases with decreasing 
ionic strength, namely 6 ^ 145° at 80 mM, 6 ^ 135° at 15 mM and 6 ^ 95° at 5 mM 
|23j . In the biochemical context the change of 6 is accomplished by other mechanisms, 
especially by the depletion of linker histones and the acetylation of core histone tails 
(cf. my discussion in Section 3.5), both of which are operative in transcriptionally active 
regions of chromatin. These mechanisms lead effectively to a decrease of 6. 

As pointed out below Eq. [HHl the decrease of 6 is accompanied by a decrease of the 
line-density pi = n/d of nucleosomes at an essentially fixed value of d. In other words, 
the number of vertices of the projected polygon decreases significantly with decreasing 
6 because 6'„ = vr (1 — 1/n). In that respect the effect of reducing 6 below the optimal 
packing value might be best viewed as an "untwisting" of the 30-nm fiber. Using the 
experimentally determined values of 6 one finds that the line density (the number of 
nucleosomes per 11 nm) is given by pi ~ 5.7 for 6 ^ 145°, ~ 4.3 for 6 ^ 135° and 
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Pl ~ 2.0 for 6 ^ 95°, which is in reasonable agreement with the experimental values 
Pl ~ 6.0, Pl ~ 3.2 and pl ~ 1-5 [5^. Furthermore, the number of polygonal vertices 
n = -K / [ti — 9) decreases as follows: n ^ 5.1 for 9 145°, n ~ 4.0 for 9 ^ 135° and 
n 2.1 for 9 95°, consistent with the stereo pair images by Bednar et al., suggesting 
n ~ 5 at an ionic strength of 80 mM and n ~ 3 at 5 mM (cf. Figs. 3(a) and (b) in |23j). 

Let me close this section with a cautionary remark The 3d density and the line 
density of the fiber can not only be changed by changing 9 ot (p but also by changing the 
linker length (in multiples of 10 bp). A variation in B changes the location of the point 
(^max) 0min) in the diagram of geometrical states, and thus the values of the maximum 
3d and line densities that can be achieved, namely 

(max) ^ ]^ ^ ^ (90) 

2VS(t>nnn{n-9,,,^,)B^ 2^RIB 

and 

(max) 4 TT^ 

- (tT - ^^ax) " 4 ^ ) 

This shows that fibers with smaller values of B can achieve higher 3d densities but have 
a smaller maximal line density (and accessibility dpL/d9 oc B"^). From this one might 
infer that active cells should have larger nucleosome repeat lengths in order to maximize 
the accessibility to their genetic material. An overview on nucleosome repeat lengths in 
different organisms and tissues is given in table 7-1 of van Holde's book ,14j. The data 
shown there do not follow this rule, unfortunately. In fact, very active cells like yeast 
cells and neuronal cells have in general short nucleosome repeat lengths while inactive 
ones like sperm cells have large ones. This shows that the optimization principle of high 
density and accessibility has to be used with caution. 



3.4- Mechanical properties of the two-angle model 

The two angle model - as discussed in the previous sections - is purely geometrical. 
Could it be useful as well for predicting physical properties of the 30-nm fiber? The 
response of the 30-nm fiber to elastic stress was indeed one of the major issues in 
our paper on the two-angle model In an independent study on the two-angle 

model by Ben-Haim, Lesne and Victor |H1] this question has been the major focus. By 
combining in this section results from both papers will allow for the first time to give 
analytical expression for the elastic properties of the two-angle model as a function of 
the underlying pair of angles 9 and 0. 

Before doing so let me remark that the elastic stress can either be of external 
or of internal origin. External stresses are exerted on the chromatin during the cell 
cycle when the mitotic spindle separates chromosome pairs |194j . The 30-nm fiber 
should be both highly flexible and extensible to survive these stresses. The in vitro 
experiments by Cui and Bustamante demonstrated that the 30-nm fiber is indeed very 
"soft" [HI]. The 30-nm fiber is also exposed to internal stresses. Attractive or repulsive 
forces between the nucleosomes will deform the linkers connecting the nucleosomes. 
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For instance, electrostatic interactions, either repulsive (due to the net charge of the 
nucleosome core particles) or attractive (bridging via the lysine-rich core histone tails 
[10]) could lead to considerable structural adjustments of the model. 

Before considering the elastic properties of the two-angle model, it is helpful 
to briefly recall some results concerning the large-scale elasticity of the DNA itself 
|1681 I195j . The measured force-extension curve of naked DNA breaks up into two 
highly distinct regimes: the "entropic" and "enthalpic" elastic regimes. For very low 
tension F (< pN), the restoring force is provided by "entropic elasticity" [ZHj- In the 
absence of any force applied to its ends, the DNA's rms end-to-end distance (chain 
length, L) is small compared to its contour length (Lq) and the chain enjoys a large 
degree of conformational disorder. Stretching DNA reduces its entropy and increases 
the free energy. The corresponding force / increases linearly with the extension L: 

F-^^, L«Lo (92) 

with lp ^ 500 A being the thermal persistence length of DNA [HH] . 

For higher forces [F > 10 pN), the end-to-end distance L is close to Lq and the 
elastic restoring force is due to distortion of the internal structure of DNA. In this 
regime, the force extension curve can be approximated by 

F ^ kpTj^^^, L>Lo (93) 

The stretching modulus 7 = (df/dL) Lo/kpT of DNA is about 300 nm-^ HHSl IW)] . 
i.e., almost four orders of magnitude larger than the corresponding value 3/lp obtained 
from Eq. EH 

In the following I shall first discuss how the mechanical properties of the linker 
backbone (modelled as a two-angle fiber) can be derived analytically from its geometry. 
Then, in Section 3.4.2 the influence of nucleosome-nucleosome interaction is considered 
before I compare in Section 3.4.3 the theoretical results with that of stretching 
experiments on chromatin fibers |8T| l82| l83] . 

3.4-1 ■ The elasticity of the linker backbone That the chromatin fiber is highly flexible 
due to the large amount of twistable and bendable linker DNA has been pointed out 
by myself, Gelbart and Bruinsma |24J. For a few special cases we were also able to 
calculate the stretching modulus of the two-linker model. A complete analysis of the 
elastic properties of the two-angle model has been given by Ben-Hai'm, Lesne and Victor 
|M] . In that paper the authors managed to relate the macroscopic mechanical properties 
of the fiber to the geometrical properties of the master solenoid (i.e. to quantities like R, 
So and ip) . Their underlying microscopic geometrical model was more complicated since 
it was assumed that the linker DNA leaves the octamer as a straight line so that entering 
and exiting strand are displaced with respect to each other. A similar arrangement has 
also been assumed in the original study by Woodcock et al. |2IIj . From cryo-EM pictures 
it is known, however, that in the presence of linker histones the entering and exiting 
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strand are glued together in a stem [221 and this is also the situation encountered in the 
mechanical stretching experiments by Cui and Bustamante [HI]. Therefore it might be 
more appropriate to model the influence of the nucleosome on the linker DNA just as 
inducing a kink on the DNA - as modelled in the above discussed variant of the two- 
angle model. Since for this case we have the exact relations between the geometrical 
parameters of the master solenoid (i?, Sq and ■?/') and the underlying two-angle geometry, 
Eqs.|7nito[7Sl the problem of calculating the mechanical properties of the two-angle fiber 
is now completely analytically solved. 

Let me sketch in the following the elegant line of arguments used by Ben-Hai'm, 
Lesne and Victor [HI] to determine the mechanical parameters. The basic idea is that 
the two-angle fiber can be described as an extensible WLC |168| I169| ll7Uj as already 
suggested in Refs. [HH IHE] . In the hnear response regime the relation is 



I F \ I knT^ ksTg \ 



/ 



u 



\ 



Mt = kBTg C 0^ 
\Mk I \ A J \ R-' J 

Here F, Mf and denote the external force and torque components: F is the force 
along the fiber axis, Mt the torsional torque (the torque component of M parallel to 
the fiber axis) and Mj, denotes the fiexural torque which is the torque component 
perpendicular to the fiber axis. These stresses are linearly related to the strains: u 
is the relative extension, Q the twist rate and is the curvature of the fiber. The 
components of the stress-strain tensor give the mechanical properties of the fiber: the 
stretching modulus 7, the bending stiffness A, the torsional stiffness C and the twist- 
stretch coupling constant g. These quantities follow from the underlying properties of 
the linker DNA that is modelled as a non-extensible WLC (like in Eq. In Ref. [HI] 
the authors wrote down the energy density as a function of the stresses (and not of the 
strains as usual). Then they compared the resulting energy per linker with the energy 
that follows from a microscopic calculation of the fiber elastic energy (again as function 
of the applied stresses). The microscopic calculation was based on the equilibrium 
condition for Kirchhoff rods (WLCs) |125j . applied to the linker DNA: The force f (s) 
on the linker at any given point r (s) (s: arclength) equals the external tension 

f(s) = F (94) 

and the local torque obeys 

m(s) = M- V AF (95) 

where v is the vector pointing from the fiber axis to the point r(s). With these 
assumptions it was possible to obtain analytical expressions for the mechanical fiber 
properties [HI]. Specifically: 

sn C + AScos^z 
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tan^ (r//2) + A - AS cos^ z] f [rj, z) 



(98) 



and 



So AS* cos z sin z 



(99) 



^ = 



keTB i?cos(r//2) 



where* 



/ (^, z) 



3A 



(100) 



+ tan2 {r]/2) {C + cos^ z) 



Note that all the parameters occurring in Eqs. IHSl to llOOl can be deduced analytically 
from the two-angle geometry. Specifically R is the fiber radius, Eq. EHl and z denotes 
the angle between the fiber axis and the linker, z = arccos (so/-B) with sq given by 
Eq. [7H| Furthermore t] = cot (ip) Sq/R is the angle between neighboring nucleosomes as 
seen when viewed down the fiber axis, i.e. rj/so is the twist rate of the unperturbed 
fiber. From Eqs. [TBI to [7H1 follows: 



The other parameters describe the mechanical properties of the DNA: the bending 
stiffness A and the torsional stiffness C as well as their difference AS* = A — C. Therefore 
we know now the macroscopic mechanical properties of the two-angle fiber as explicit 
functions of the microscopic parameters. 

These functions, 7 = 7(^,0), A = A{9,(f)), C = C{9,(j)) and g = g{9,(f)), are, 
however, rather unwieldy. To get an idea of the overall behavior one might resort to 
numerical calculations as done in Ref jHl] where the mechanical moduli of the fiber 
were calculated as a function of (j) for two values of 6, cf. Fig. 11 in that paper. It 
was found that the moduli vary strongly with (p (and thus with the linker length) and 
it was argued that this strong dependence might be used in the biological context as a 
regulatory factor. 

Having the analytical relations at hand, another approach has now become 
available, namely to look at limiting cases (solenoids, fibers with crossed linkers and 
zig-zag structures) which show simple dependences on the underlying geometry (i.e., on 
the angles 6 and 0) - as discussed in Section 3.2. Not surprisingly 7, A, C and g are 
also simple functions of these underlying angles in all the limiting cases. We will give a 
complete overview in a forthcoming publication |192j . Here I will restrict myself to two 
limits only. 

Let me start with the planar zig-zag fiber. Such a chain can be stretched via the 
bending of its linkers - maintaining the defiection angle 6 at each kink, cf. Fig. 21. It 
is also clear that the linker will not be twisted in this planar geometry. In Ref. 
we calculated the stretching modulus 7 for this special arrangement. In order to do 
so we wrote down the elastic energy of the linker (similar to Eq. and determined 

* In the original work |84j this function is caUed K (??,3). Note that there is a printing error in the 
denominator (last line of Eq. (16) in that paper). The factor 1/3A has to be removed. 
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F=0 



(a) 




(b) 



Figure 21. Stretching of a zig-zag chain, (a) The unperturbed chain, F = 0, with 
a total length has straight linker DNA. (b) The same fiber under tension F > Q 
stretches to an end-to-end distance L > Lq via the bending of its linkers. 



the deformed shape (for small perturbations) from the corresponding Euler-Lagrange 
equation taking the boundary conditions into account. We found 

^_ mcos(g/2) 

^ ksTB^ sin' (6/2) ^ ' 

This result can now also be obtained directly from the general formula, Eq. EHl by 
simply setting = vr. It is evident from Eq. 11021 that the stretching occurs via linker 
bending since 7 depends on A only whereas the DNA torsional stiffness C does not 
enter the expression. I also note that the planar zig-zag fiber shows other interesting 
features, especially two different persistence lengths for bending in plane and bending 
out of plane as discussed in Refs. [211 Ell 1197] . General features of such polymers with 
highly anisotropic bending rigidities have been considered by Nyrkova et al. |198j . 

Now I consider the case that might be of importance for 30-nm fibers: the chains 
with crossed linkers (0 -C 1, vr — ^ <^ vr). Starting from the general expressions it 
is straightforward to show that the linker geometry leads in this case to the following 
overall mechanical properties: 
3v4 _ 
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It can be seen from Eq. 11031 that stretching occurs via hnker bending (as in the case of 
zig-zag fibers) and from Eg. llOSI that also the twisting of the overall fiber is achieved via 
the bending of the linkers. The dependence of A on the DNA parameters, Eq. Il()4| shows 
that fiber bending involves both bending and twisting of the linkers, a fact that is due 
to the different orientations of individual linkers with respect to the bending direction. 
Finally, the twist-stretch coupling is very small (cf. the angle dependence in Eq. I106j) . 

The elasticity of the linker backbone is predicted to be very soft: For instance, 
the stretching modulus 7 scales for fibers with crossed linkers and zig-zag chains as 
A/ [kBTB^). This is of the order one (per nm) for an effective linker length of 20 bp 
as compared to a 7 of ~ 300 nm~^ for free DNA (see above). Of course, depending 
on the values of 6 and this value varies in a wide range. Also the other mechanical 
parameters of the two-angle fiber indicate an extremely soft structure. Because of this 
it is evident that the presence of the nucleosomes play a crucial role in determining the 
mechanical properties of the 30-nm fiber. The excluded volume will not allow a strong 
bending that would lead to overlapping nucleosomes and the nucleosome-nucleosome 
attraction counteracts the stretching of the fiber under external tension. Therefore 
before I compare in Section 3.4.3 the theoretical expressions and the results from fiber 
stretching experiments it is indispensable to discuss first how the nucleosomes modify 
the mechanical properties of the two-angle fiber. 

3.4-2. Role of the nucleosome interaction The effect of attractive interaction between 
nucleosomes is to cause a compression of the 30-nm fiber. Phase behavior studies of 
linker-free nucleosome solutions, i.e., solutions of disconnected nucleosomes |199j (cf. 
also Ref. |2()()j ) indicate that nucleosome core particles spontaneously form fiber-like 
columnar structures, presumably due to attractive nucleosome-nucleosome interaction. 
Attractive nucleosome interaction could be mediated for instance by the lysine-rich core 
histone tails [10], as mentioned above. 

Let me first discuss the role of this internucleosomal attraction on the stretching 
elasticity of a fiber. Following Ref. the special case of a planar zig-zag structure 
with elastic linkers is considered where a short-range interaction between nucleosomes 
is assumed in addition. This interaction, denoted by Uinter, is taken to be a short 
range attraction, of strength —Umin, that acts only when the nucleosomes are in close 
contact, i.e., at a distance x ~ 2i?o of the order of the hardcore diameter. For a given 
nucleosome, say the ith, the closest nucleosomes in space are number i + 2 and i — 2 (cf. 
Fig. 21). The interaction between other pairs is disregarded. The elastic interaction 
Ueiastic follows directly from Eq. 11021 applied to a trinucleosome (A^ = 2): 



elastic 
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Figure 22. (a) Intcrnucleosomal interaction potential U between nucleosome i and 
I + 2 as a function of distance x. In addition to the elastic contribution there is a 
short range attraction for nucleosome at close contact, x = 2Rq. The different curves 
correspond to different values of the angle 6. Curve "1" has the global minimum at 
large x (swollen state "S") whereas curve "3" has the minimum for nucleosomes in 
close contact (condensed state "C"). Curve "2" corresponds to the transition point. 
Also depicted is the common tangent for curve " 3" . Its slope corresponds to the critical 
stretching force fcs at which nucleosomes are transferred from the state C to state S. 
(b) Force-extension curve of a condensed fiber. For extensions L with Li < L < L2 
one finds a coexistence plateau with the restoring force fcs- 



where Xq = 2B cos (0/2) denotes the distance between nucleosome i and i-|-2 for straight 
hnkers (cf. Eq. |HH). The total internucleosomal U (x) equals Uinter (x) + UeiasUc {x). 

Fig. 22(a) shows U (x) for different values of 9. Let me assume for simplicity that 
the interaction energy Uinter remains unchanged. Curve "1" in Fig. 22(a) shows U (x) 
for a small value of 9 where the global minimum of U (x) is located at x = xq denoted 
by "S" (swollen state). Curve "2" corresponds to an intermediate value of 9 at which 
the minima at "S" and "C" have the same value. For this value of 9, 9 = 9c, the energy 
minimum shifts from " S" to a new minimum, representing the condensed state " C" . The 
change in 9 produced a structural transition from a swollen state to a condensed state. 
Finally, curve " 3" depicts U (x) for a deflection angle 9 > 9c with the minimum at " C" . 
The critical angle for the "S" to "C" transition can be determined by comparing the 
bending energy at close contact, UdasUc (2-Ro), and the strength Umin of the short range 
attraction. Equating both leads to the following condition for 9c- 
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In the swollen state the elastic properties should be in principle the one that were 
discussed in the previous section. In the condensed state, the elastic properties are 
determined by the detailed form of the nucleosome interaction potential. 

If the condensed state has a lower free energy, i.e. if ^ > 9c-, then an external 
stretching force / can induce a transition from the condensed to the swollen state. The 
transition point fcs follows from a "common-tangent" construction. The conditions are 



U' (xi) = U' (X2) = fcs and {U {x^) - U (xi)) / {x^ - x^) = fcs (cf. Fig. 22(a)) leading 



The corresponding force-extension curve has a "coexistence plateau", cf. Fig. 22(b). If 
the imposed end-to-end distance is smaller than Lq (the contour length of the condensed 
fiber) then the restoring force is entropic. For Lq < L < Li the force rises sharply with 
increasing L. This "hard elasticity" is governed by the nucleosomal interaction potential 
U inter- Then at L = Li the coexistence plateau is reached. Between L = Li and L = L2 
parts of the fiber are in the "S" state and parts are in the "C" state. For larger 
extensions, L > L2, the fiber shows soft elasticity due to the bending and twisting of 
the linkers as discussed in the previous section. 

Katritch, Bustamante and Olson |HE] presented a Monte-Carlo simulation of the 
chromatin fiber that was based on a model very similar to the two-angle model. The 
nucleosomes were modelled as spheres and attached to the kinks in the linker backbone 
via a short stem. The only difference from the above discussed two-angle model was that 
the rotational angle between each pair of nucleosomes was chosen randomly from the 
interval — vr to tt. These fibers were then stretched as in a micromanipulation experiment 
[HT] and their force-extension relationships were measured. The values for 7 were in good 
agreement with what is expected on theoretical grounds (a detailed discussion is given 
in Appendix D of Ref. ^5). What is of special interest here is that they also studied 
the effect of a short-ranged nucleosome-nucleosome attraction. Using a value Umin of 
order 2/cbT (or larger) they observed very clearly the occurrence of a pseudoplateau in 
the force-extension curve similar to Fig. 22(b). 

The nucleosomes have also a large effect on the persistence length Ip of the fiber. 
This has been demonstrated most clearly in a computer simulation by Wedemann and 
Langowski jHZj (cf. also an earlier preliminary study of this group |2nij ). Their model is 
again very closely related to the two-angle model discussed above. Differences are that 
the entering and exiting DNA at the nucleosome are slightly displaced in the direction 
of the nucleosome axis and that the screened electrostatic interaction between linkers 
was taken explicitly into account. Nevertheless, Eq. IHZI should be expected to give a 
good estimation of the contribution of the linker DNA to the fiber persistence length. 
Using the values of that simulation (6* ~ 143°, ~ 80°, B = 10 bp) gives Ip ~ 13 
nm. However, the persistence length observed in the simulation is 265 nm, i.e. 20 times 
larger! This is clearly an effect of the nucleosomes. The role of the linkers is to bring 
the nucleosomes into contact. The nucleosomes (modelled here as elhpsoids) experience 
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then in addition an attractive force where Umin has been chosen to be of order ksT. This 
leads to a very dense structure with the nucleosomes in contact so that there is hardly 
any space for fiber bending. Most clearly this is seen in Fig. 9 of that paper that shows 
a contraction of a fiber that has been stretched out first. As long as the nucleosomes 
are not in contact the fiber shows sharp bends and strong shape fluctuations. The fibers 
stiffens very strongly as soon as the dense state is reached. 

This all shows that the nucleosome interaction is a crucial element determining the 
mechanical fiber properties. Therefore a more microscopic model taking details of the 
nucleosome structure into account might be important for a theoretical prediction of the 
properties of chromatin fibers. A first step in this direction has been done by Beard and 
Schlick |2U2j . They performed a molecular dynamics simulation where the nucleosomes 
were represented by disks made of several hundred charges that were chosen to match 
the crystal structure jTHI. Di- and trinucleosomes as well as whole fibers have been 
studied. The authors demonstrated also that a fiber with a crossed linker geometry 
unfolds into an open zig-zag fiber as a result of changing ionic conditions. However, 
in their nucleosome model they neglected most of the histone tails that constitute very 
likely the crucial ingredient for the nucleosome-nucleosome attraction. 

3.4-3. Stretching chromatin It is now possible to measure the mechanical behavior of 
single chromatin fibers via micromanipulation techniques as it has been demonstrated 
in three studies jHH E21 EHj- Each of the studies focused on a different variant of the 
30-nm fiber. Cui and Bustamante [HI] stretched native chicken erythrocyte chromatin 
fibers containing linker histones and contrasted the cases of low and high ionic strength. 
For low ionic strength (5 mM NaCl) it was found that fibers are very soft. By fitting 
their data to that of an extensible WLC they found a stretching modulus of ksT'y ^ 5 
pN and a persistence length of Ip ~ 30 nm. The theoretical values are ksT^ ~ 6.3 
pN (from Eq. and Ip ~ 16 nm (from Eq. W7\ for (p = 36° and 6 = 95° (cf. Section 
3.3) and for a linker length of 20 bp (chosen on the basis of the 210 bp repeat length 
of chicken erythrocyte chromatin ^3] minus roughly 190 bp associated with the core 
and linker histones). The theoretical and experimental values are close which indicates 
that the mechanical properties of a swollen fiber at low salt concentrations are mainly 
determined by the elasticity of its linker backbone. The nucleosomes are less important 
since they are not close enough in such a swollen fiber. For high stretching forces around 
20 pN there is an irreversible change in the overall length of the fiber due to histone 
"evaporation" which has been seen more clearly in the other two stretching experiments 
(see below). 

Stress-strain curves for fibers at higher ionic strength (40 mM NaCl) are also 
reported in Ref. |HI]. In this case the fiber is much denser and nucleosomes approach 
each other closely. Attractive short-range forces and the increase of 6 associated with 
higher ionic strength should favor the condensed phase. Indeed a plateau appears at 
5 pN in the force-extension curve (cf. Fig. 4 of Ref. [HI])- From the extend of the 
plateau, 0.6 fim, its height, 5 pN, and the number of nucleosomes in the stretched fiber, 
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^ 280, it was estimated that there is an attractive interaction energy of ~ 3kBT per 
nucleosome |HI]. In Ref. [21] we used Eg. I109l to independently estimate the strength of 
the internucleosomal attraction from the value of the critical force alone. Neglecting the 
second term in that equation one finds Umin ~ fcs/ (^-^) ~ S/c^T assuming 6 = 140° 
and again B = 20 bp. Note that one finds with these values that the stretching modulus 
of the linker backbone, Eq. inSl is of order ksTj ~ 4.7 pN, i.e., even lower than the value 
6.3 pN predicted above for low ionic strength. The fact that these fibers appear much 
stiffer with respect to (small) deformations indicates that its mechanical properties are 
mainly determined by the nucleosome-nucleosome attraction and not by the backbone 
elasticity. This is also in accordance with the estimated large value of the persistence 
length of condensed chromatin fibers (~ 260 nm, cf. |203j ). 

Bennink et al. jHSj assembled chromatin fibers by exposing A-DNA to Xenopus 
egg extract before stretching them.jj This extract contains core histones but no linker 
histones. There is, however, an abundance of non-histone proteins some of which 
act close to the DNA entry-exit point similar to linker histones. For small forces 
a stretch modulus ~ 150 pN was extracted from the data. This indicates again 
that chromatin fibers are quite stiff at high salt concentrations (here 150 mM NaCl) 
compared to the pure linker backbone elasticity but still quite soft compared to naked 
DNA {ksT'y ~ 1200 pN). For stretching forces of order 20 pN irregular sawtooth-like 
fluctuations were observed, each being a result of a sudden fiber lengthening by multiples 
of ~ 65 nm. This was attributed to the unraveling of single or multiple nucleosomes. 

Finally, Brower-Toland et al. used well characterized fibers that were prepared 
from tandem repeats of the 5S rDNA positioning sequence and core histones (no linker 
histones). As shown in Ref. [22] (cf. my discussion of that paper in Section 2.4.1) most 
of the nucleosomes are localized at the preferred positions. Around 20 pN the force- 
extension curve showed a regular sawtooth pattern reminiscent of the one observed 
during the unravelling of tandem repeat domains in the protein titin |206[ I2()7j . The 
spacing of the peaks, ~ 27 nm, is indicative of the unravelling of only one turn of the 
nucleosomal DNA. The outer sections of the DNA were detected to be released at much 
smaller forces. Brower-Toland et al. jHS] gave also a theoretical explanation of this 
observation: They speculate that the first 76 bp are unwrapped much more easily due 
to weaker binding between DNA and the octamer whereas strong binding sites occur as 
soon as these first 76 bp are unwound. Only when a sufficiently large force is applied 
these binding sites are broken on the time scale of the experiment. 

However, this explanation is questionable since it does not take into account the 
actual unwinding geometry. Already Cui and Bustamante jH] pointed out that it 
requires a twisting of the core particle to unwrap the inner part of the nucleosomal 
DNA. The free DNA has to be bent strongly close the nucleosome which leads to 
a considerable barrier that has to be crossed during unwrapping. We are currently 
calculating this barrier analytically in the WLC framework |208j . This high barrier 

jj A detailed discussion of the kinetics of the chromatin assembly in this kind of experiment is given in 
Ref. |2M|; cf. also Ref. [Ml- 
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Figure 23. (a) Schematic view of a section of the 30-nm fiber (for simplicity shown 
here as a two-dimensional zig-zag), (b) Enlarged view of the stem region showing a 
speculative model of the role of the HI histone and some N-tails from the core histones, 
cf. also Fig. 7 of Ref. [THl] . 

might explain why the nucleosomes are dissociated only at surprisingly high tensions 

- even in the absence of linker histones. In fact, as Marko and Siggia pointed out in 
Ref. ^j, one would find nucleosome release (and a corresponding plateau in the force- 
distance curve) around 2 pN - if such a barrier could be neglected. This value follows 
from the comparison of adsorption energy (~ 30A;bT, cf. Section 2.1) to wrapping length 
(~ 50 nm): / ^ 30/^^^/50 nm^ 2 pN. 

3.5. Fiber swelling 

In this section I discuss how the entry-exit angle a = n—O of the DNA at the nucleosomes 
is controlled via electrostatics. A more detailed account on this subject is provided in 
Ref. ini]. As mentioned above it can be seen in the cryo EM studies that the 
fibers open up and therefore become more accessible when the ionic strength is reduced 
and that this opening is directly linked to an increase in a. It was suggested that via 
other mechanisms (for instance, the acetylation of the lysine-rich histone tails 17j, as 
explained in more detail below) the angle a and therefore the degree of swelling can be 
changed for a given section of the fiber and that this constitutes a biochemical means 
to control the transcriptional activity of genes. 

Whereas the X-ray studies of the core particle [10] allow a detailed knowledge 
of the wrapped part of the DNA it does not give insight into the conformational 
properties of the entering and exiting strands. One has therefore to refer to the electron 
cryomicrographs. In these micrographs it can be seen clearly that 10 nm stretches of 
the entering and exiting DNA strands are glued together forming a unique " stem motif" 

(cf. also Fig. 23(a)). The gluing of the two equally charged chains is accomplished 

- amongst other things - via the linker histone HI as shown schematically in Fig. 23(b). 

At physiological concentrations the electrostatics is essentially short-ranged (k^^ ~ 
10 A for 100 mM salt). It seems therefore reasonable to assume that a is set within the 
small region where the two linker DNA are in close contact, i.e., within the stem region. 
This value of a in turn controls the large-scale secondary structure of chromatin, the 
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Figure 24. Idealized model for the entry-exit region of the DNA at the nucleosome. 
The thick curved lines represent the entering and exiting DNA that enclose a well- 
defined angle a that in turn determines the overall geometry of the 30-nm fiber. 



30-nm fiber, as discussed in Section 3.2. To mimic this situation I assumed in Ref. jHI] 
a geometrical arrangement where two parallel DNA strands are hold together tightly 
at y = for X < and are free for x > 0, cf. Fig. 24. Because of their mutual 
electrostatic repulsion the two strands bend away from each other. When the two strands 
are far enough from each other their interaction is screened so that they asymptotically 
approach straight lines defining the opening angle a as indicated in Fig. 24. 

The conformation of the upper DNA chain can be described by the height function 
h {x). By symmetry the position of the lower strand is then given by —h (x). To mimic 
the stem two boundary conditions are imposed at x = 0, namely h (0) = h' (0) = 0. 
The entry-exit angle a is related to the slope of h (x) at infinity via tan (q;/2) = h' (oo). 
The two DNA chains are modelled as semiflexible polymers with persistence length Ip 
and line-charge density —e/h that interact via a screened electrostatic potential. The 
free energy of the system is then given by 

2 



J^{hix)} , 
^' ' ax 



kuT Jo 




(110) 



The first term in the integral accounts for the bending of the two DNA strands and 
the second term describes the interaction between the two chains {Kq (x) being the Qth 
order modified Bessel function). Here the interaction of a given charge on one chain 
with all the charges on the other chain is approximated by the interaction of this charge 
with a straight chain at the distance 2h.]] The conformation of the upper chain, h (x), 
is then the solution of the corresponding Euler-Lagrange equation 

- ^i^i (2-/^) = (111) 

ft This can be shown to be a good approximation as long as Zp loSF = ?b/46^k^ (the Odijk- 
Skolnick-Fixman length that describes the electrostatic stiffening of the chain, cf. Ref. |2U9| ). In fact, 
it is also precisely that limit at which the mframolecular interaction can be neglected (as done here). 
One can reformulate the above condition in the simple form a ^ 1 (cf. Eg. 11121 below): for instance it 
is found that for a — 45° the approximation is still excellent and for a — 90° the chain-chain repulsion 
is overestimated by ~20%. 
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together with four boundary conditions: Two are given at the origin (see above) and 
two follow from the condition of straight "linkers" at infinity: h" {oo) = h'" {oo) = 0. 
Defining h = 2Kh and introducing the dimensionless quantity x = {AIbh"^ /b'^lp)^^^ x, 
Eg. II 111 can be rewritten as d'^h/dx^ = Ki (fij with a solution showing the dimensionless 

asymptotic slope cq= dh/dx _ .It follows immediately that tan(Q;/2) is given by 




tan(a/2) = /.'(oo) = ^ f - (112) 



where cq is of order one (cf. Ref. jHI] for details). 

This result might help to understand better how the local electrostatics controls 
the geometry of the chromatin fiber in vitro and, on a more tentative level, in vivo. 
The in wiro-experiments show that chromatin fibers "open up" with decreasing salt 
concentrations. As already mentioned in Section 3.3 it was estimated from electron 
cryomicrographs that Uexp ~ 85° for = 5 mM and Uexp ~ 45° for Cg = 15 mM and from 
electron cryotomography that aexp ~ 35° for Cg = 80 mM [23J (cf. Refs. [T7] and |210j for 
other approaches to determine a). One expects from Eg. II 121 that a ~ 2 arctan (^Cc~^^^^ 
with C being a constant. Let me take the angle at the highest salt concentration, 
Cg = 80 mM, as the reference value. From this follows C = 0.94. With this value of 
C the prediction is a ~ 51° for Cg = 15 mM and a ~ 64° for = 5 mM. Whereas 
the theoretical value a ~ 51° at intermediate ionic strength is close to ae^p ~ 45°, the 
value a ~ 64° for low salt concentrations is noticeably too low {aexp ~ 85°). However, as 
mentioned above, for such a large value of a the chain-chain repulsion is underestimated 
by -20%. 

How can the degree of swelling of the chromatin fiber be controlled in vivol Under 
the assumption that the above mentioned geometry is valid the only parameter that 
might be under biochemical control is the linear charge density 6^^. It is known that the 
formation of a dense chromatin fiber is dependent on the presence of several components, 
especially of the cationic linker histones and of some of the lysine-rich (i.e., cationic) 
N-tails of the core histones that appear to be long, flexible polyelectrolyte chains |10j . 
In Fig. 23(b) I give a tentative picture of the conformation of two N-tails that protrude 
from the histone core. It is known that if either of these components is missing the 
fiber does not fold properly (cf. Ref. J7] and references therein). As indicated in the 
Figure the tails might form a complex with the entering and exiting linker DNA in 
such a way that they effectively reduce its linear charge density h^^. It is known that 
transcriptionally active regions in chromatin show an acetylation of the core histone 
tails (i.e., the cationic groups of the lysines are neutralized). In that tentative picture 
this acetylation mechanism would increase and according to Eg. I112l this would lead 
to an opening of the entry-exit angle a. The acetylation might therefore be the first 
step in the decondensation of a stretch of the chromatin fiber that needs to be accessed 
for transcription. 

Let me note that processes that are involved in the acetylation and deacetylation 
are guite specific and involved as, for instance, discussed in Ref. |211j . The 
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histone tail modifications might serve specific functions via the modification of their 
secondary structure that in turn modifies their interaction with certain proteins |212j . 
Recently there are even attempts to decipher a specific "language" of covalent histone 
modifications |213j . It might be that such specific processes act in concert with the 
more basic charge neutralization principle discussed here. 

4. Conclusions and outlook 

Chromatin is of fundamental importance for a host of biological processes ranging 
from gene expression to cell division. Consequently there is a huge research activity 
among biologists in this area. For physicists chromatin becomes now also of interest 
since there are more and more experiments available that work under quite well- 
defined conditions. Such experiments typically involve only a few components (DNA, 
histone proteins...) but no active protein "machines." These experiments either focus 
on elucidating properties of single nucleosomes or of beads-on-a-string complexes 
("chromatin fibers"). They study the behavior of these systems under changing ionic 
conditions and/or under an externally applied tension. Also the djTiamics of these 
systems can be investigated that is solely driven by thermal fluctuations. Theoretical 
treatments and computer simulations that capture the essential features of the chromatin 
system are now possible and thus allow to estimate the energy- and time-scales occurring 
in chromatin. Other approaches look at simplified model systems and try to identify 
general physical principles that govern complexes of charged chains and macroions. With 
the better understanding of the mechanical and dynamical properties of nucleosomes 
and chromatin fibers one hopefully gains also a deeper insight into more complicated 
questions like the working of chromatin remodeling complexes, the interaction between 
RNA polymerase and nucleosomes etc. 

To proceed in this direction it is crucial to obtain reliable numbers from 
experiments. One energy scale that dominates many processes is the adsorption energy 
of DNA on the octamer that has now been measured quite directly through stretching 
experiments. Another important feature - especially in 30-nm fibers - is the nucleosome- 
nucleosome interaction energy. Again detailed experimental studies have been performed 
and await a detailed theoretical treatment. 

Chromatin is a very active and exciting field in biology where tremendous progress 
has been made in recent years. I hope that at least a few of the ideas gained from the 
physical models will be of help to biologists to develop a clearer picture of the working 
of chromatin and to design appropriate experimental setups. 
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A. Rosette in d dimensions 



The heterogeneity found for open loops (cf . Eq. is reminiscent of a phase coexistence. 
To make clearer why the loop sizes are so sensitive to small details, especially why there 
is no "phase separation" for closed chains (cf. Eq. |^ even though <ff/dP < at 
larger separations) I will present here some unpublished results in which the free energy 
is formally recomputed for arbitrary space dimensions d. In order to do so one has to 
replace the entropy term in Eq. IHHlby {d/2) In {l/lp)- Also in that case one obtains an 
analytical expression for the partition function, namely (for an open chain) 

'2xkBT' 



d+2 

2/p* 



2-d 
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ksT 
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with Ky (x) denoting the modified Bessel function of vih order. In the case of large 
"pressure" P, IpP/ksT 3> 1, it follows from the asymptotic form of Kd/2-i{x) for 
large x, Kd/2-1 {x) ~ \Jti /2xe~^ , that the leading term of the resulting free energy 
G {P) is independent of d. Therefore one recovers the 3D case, i.e., Eqs. for 



L/M <^ Mlp. For the low "pressure" regime, IpP/ksT -C 1, I use the asymptotics 
Kd/2-iix) ~ 2-ir(|rf/2-l|)(2/x)''^/^"^' for X < 1 and ^ 2. This leads to the 
following asymptotic behavior of the partition function 
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It follows then from L = dG/dP that P is given in leading order by 
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The average leaf size can in principle be calculated, as before, from Zi. Here, however, 
it turns out to be more convenient to calculate {keaf) directly: 

/o-^n^-^/^exp(-M.-^) 
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Now using the above given power law behavior of the Bessel function together with 
Eq. IA3I the average leaf size follows: 
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First note that for d < 2 the leaf size is set by the overall length of the chain but 
does not depend on /p; on the other hand, for > 4 {lieaf) is solely determined by Ip. 
Speaking in the picture of interacting particles on a track of length L one can explain 
these two extreme cases as follows. For d < 2 the increase of the "nearest neighbor 
pair potential" beyond a distance Ip (given by {d/2)ln{l/lp)) is too small to keep the 
particles together; instead they explore all available space. For higher space dimensions 
than 4 the prefactor of the log-term is large enough to keep neighboring particles close 
to the ideal distance ~ Ip given by the shallow minimum of / (/), Eq. IHHl The case 
d = 3 which I already gave in Eq. (a result recovered in Eq. IA5|) is an intermediate 
case where (keaf) reflects the overall chain length L as well as the position ~ Zp of the 
shallow minimum. Note further that in the limit L ^ oo the average size per leaf goes 
to infinity for d < 4 (but the "particles" will only be spread out over the whole volume, 



Concerning the role of dimensionality one also gains some insight by the 
following simple argument (similar to the famous Onsager /Manning argument for 
the condensation of counterions on an infinitely long charged rod jlllj ). Consider 
a pair of two particles at distance / in one dimension that attract each other via 
/ (/) = ksT {d/2) In (l/lp). Now assume that the particles move further apart from 
the distance li to the distance I2 > h- This leads to an increase in energy by 
AE = ksT {d/ 2) In {I2/I1). On the other hand the particles gain entropy since they 
are now less confined: — /cpTAS* = fcpTln (/1//2). Hence for d < 2 the particles will 
"loose" each other since their attraction to the nearest neighbors is overruled by the 
gain in entropy - as derived rigorously in Eq. IA5I 

Finally, I mention that the same extension to arbitrary dimensions d can be 
performed for closed chains. One finds then phase separation for molten rosettes if > 4. 
More specifically, loops have a preferred spacing L/N for d < 4 and 4x1 p/ {d — 4) for 
d > A. This is different from the results on open chains, Eq. lASl in the interval 2 < d < 4; 
hence in d = 3 molten rosettes respond strongly to a cutting of the chain. 

B. Formation energy for small intranucleosomal loops 

Since the configurations of small loops are essentially planar, it is convenient to describe 
them in terms of the function r (6), cf. Fig. 12(b), where r and 6 are the polar coordinates 
of an arbitrary point on the loop (with the origin chosen on the cylinder axis and the X 
axis running through the center of the loop). In these terms the line element ds takes 
the form ds = dOJr"^ + [dr/d9)^ and the loop excess length is given by 



with 29* being the aperture angle of the bulge, 29* Rq = L*. The local curvature 1/R = 
\d'^r (s) / ds'^l of the loop takes the form 1/R = |r"r — 2r'^ — / (r^ + r'^)^^^, with 
primes and double primes denoting the first and second derivative, respectively, with 



M {keaf) ^ L, 
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respect to 6. Restricting ourselves here to small loops, one can write r {6) = Rq + u {6) 
with u <^ Rq everywhere. Keeping only quadratic terms in u and its derivatives one 
obtains ds ~ RodO (1 + w/i^o + m'V (2i?o)) and 



ds- 



de 
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The bending energy of the loop is then 

1 



U +-U 



12 



U 



+ Auu" 



u ^ 1 

Rq Rc 



E. 



elastic 



2 



de 



u 



112 



i?0 -Ro 



[B2) 



(B3) 



where the boundary conditions u iff) = du {6) /d9 = at 6 = ±6* have been used. 

The variational energy of the loop E = EeiasUc ~ TAL for a given aperture angle 6* 
and subject to the constraint of a fixed loop contour length L = L* + AL follows from 
Eqs.EUandlHHlto be: 
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Here T is the Lagrange multiplier that constrains an extra length AL to be adsorbed as 
the loop is formed (T can be interpreted as a "tension" pulling in extra length). In Eq. lB4l 
the constant terms as well as the term linear in u have been dropped since solutions to 
E with and without the linear u-term differ just by a constant. The optimal loop shape 
(with given values of AL and 9*) obeys the Euler Lagrange equation 6E/6u = 0: 
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Solutions of Eq. IB 51 are of the form u oc e*'^^ with four possible values of A 
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that show the following asymptotics: 
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One expects symmetric solutions of the form u (6) = 
The boundary conditions u = and u' = at 6 
Cicos(A+r) + C2Cos(A_r) = and CiA+sin(A+r) - 
solubility condition leads to the transcendental equation 

A_ _ tan(A+r) 

A+ ~ tan(A_6'*) 
For vanishing "tension", T 



Cicos(A+^) + C2Cos(A_^). 
= 6* have then the form 
- C2A_ sin (A_^*) = 0. The 

(B8) 

0, one finds from Eq. IB 71 that the condition IB8I is of 



the form 2tan(v^r) = tan (r/v^ that has no non-trivial solution. For large T, 



T ^ A/Rq, one obtains from Eqs. 
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The left hand side of Eq IB9I is small and hence solutions are approximately given by 



TRq/A6* ~ A+6'* ~ fcvr with k = 1,2, 3, ... In the following I consider only the k = 1 
solution which is the solution that leads to the smallest elastic energy. 

Using partial integration and Eq. IB5l the loop formation energy, Eq. IHOl can be cast 
into the form 

AU = 2kBTXRo6* + — / dOu''^ - A / d9^ (BIO) 

Rn Jo Jo Rn 



To proceed further one makes use of the explicit solution given below Eq. IB 71 Assume 
the large tension case T ^ ^/-Rq (^^^ quality of this approximation will be checked 
a posteriori); then A_ ^ 1. The condition u{9*) = takes then the form C2 — 
—Ci cos (A+^^*) which leads together with A+6'* ~ tt to C2 = Ci. The loop shape is thus 
approximately given by u (0) ^ Ci (1 + cos (ttO/O*)) from which follows its formation 
energy 

AU = 2ksTXRoe* + ^ - ^ (BID 

and excess length, Eq. IHIl AL ~ 2Ci9* + n'^Cf/ (2i?o^*) ^ 2Cie* where I used the fact 
that small loops have small amplitudes: Ci <^ Rq{9*)^. For fixed AL and 9* follows 
Ci ~ AL/ (29*). Inserting this into Eq- lBTTI and using T ^ n^A/ {Ro9*f leads to 

^^2Afi„r + ^M^-^ (B12) 

Now minimizing AU with respect to 9* (for AL fixed) gives the optimal aperture angle 

Combining Eas. lBT2l and lBT3l one arrives at the final expression for the formation energy 
of a (small) loop of excess length AL, Eq. IHTl 



References 

[1] Alberts B, Bray D, Lewis J, RafF M, Roberts K and Watson JD 1994 Molecular Biology of the 

Cell 3rd edition (Garland: New York) 
[2] Bloomfield V A 1996 Curr. Opin. Struct. Biol. 6 334 

[3] Olson N H, Gingery M, Eiserling F A and Baker T S 2001 Virology 279 385 
[4] Cerritelli M E, Cheng N, Rosenberg A H, McPherson C E, Booy F P and Steven A C 1997 Cell 
91 271 

[5] Kindt J, Tzlil S, Ben-Shaul A and Gelbart W M 2001 Proc. Natl. Acad. Sci. USA 98 13671 

[6] Kornberg R D 1974 Science 184 868 

[7] Kornberg R D and Lorch Y 1999 Cell 98 285 

[8] Khrapunov S N, Dragan A I, Sivolob A V and Zagariya A M 1997 Biochim. Biophys. Acta 1351 
213 

[9] Raspaud E I, Chaparon I, Leforestier A and Livolant F 1999 Biophys. J. 77 1547 
[10] Luger K, Mader A W, Richmond R K, Sargent D F and Richmond T J 1997 Nature 389 251 
[11] Davey C A, Sargent D F, Luger K, Maeder A W and Richmond T J 2002 J. Mol. Biol. 319 1097 
[12] Thoma F, KoUer Th and Klug A 1979 /. Cell. Biol. 83 403 



CONTENTS 



86 



[13] Widom J 1986 J. Mol. Biol. 190 411 

[14] van Holde K E 1989 Chromatin (New York: Springer) 

[15] Widom J 1989 Annu. Rev. Biophys. Chem. 18 365 

[16] van Holde K and Zlatanova J 1995 J. Biol. Chem. 93 8373 

[17] van Holde K and Zlatanova J 1996 Proc. Natl. Acad. Sci. USA 93 10548 

[18] Finch J T and Klug A 1976 Proc. Natl. Acad. Sci. USA 73 1897 

[19] Widom J and Klug A 1985 Cell 43 207 

[20] Woodcock C L, Grigoryev S A, Horowitz R A and Whitaker N 1993 Proc. Natl. Acad. Sci. USA 
90 9021 

[21] Horowitz R A, Agard D A, Sedat J W and Woodcock C L 1994 J. Cell Biol. 131 1365 

[22] Leuba S H, Yang G, Robert C, Samori B, van Holde K, Zlatanova J and Bustamante C 1994 

Proc. Natl. Acad. Sci. USA 91 11621 
[23] Bednar J, Horowitz R A, Grigoryev S A, Carruthers L M, Hansen J C, Koster A J and Woodcock 

C L 1998 Proc. Natl. Acad. Sci. USA 95 14173 
[24] Schiessel H, Gelbart W M and Bruinsma R 2001 Biophys. J. 80 1940 
[25] Workman J L and Kingston R E 1998 Annu. Rev. Biochem. 67 545 
[26] Widom J 1998 Annu. Rev. Biophys. Struct. 27 285 
[27] Polach K J and Widom J 1995 J. Mol. Biol. 254 130 
[28] Polach K J and Widom J 1996 J. Mol. Biol. 258 800 

[29] Pennings S, Meersseman G and Bradbury E M 1991 J. Mol. Biol. 220 101 

[30] Meersseman G, Pennings S and Bradbury E M 1992 EMBO J. 11 2951 

[31] Pennings S, Meersseman G and Bradbury E M 1994 Proc. Natl. Acad. Sci. 91 10275 

[32] Varga-Weisz P D and Becker P B 1998 Curr. Opin. Cell Biol. 10 346 

[33] Guschin D and WolfFe A P 1999 Current Biology 9 R742 

[34] Peterson C L 2000 FEBS Letters 476 68 

[35] Flaus A and Owcn-Hughcs 2001 Curr. Opin. Cenet. Dev. 11 148 
[36] Studitsky V M, Clark D J and Felsenfeld G 1994 Cell 76 371 
[37] Studitsky V M, Clark D J and Felsenfeld G 1995 Cell 83 19 

[38] Felsenfeld G, Boyes J, Chung J, Clark D and Studitsky V 1996 Proc. Natl. Acad. Sci. USA 93 

9384 

[39] Studitsky V M, Kassavetis G A, Geiduschek E P and Felsenfeld G 1997 Science 278 1960 

[40] Felsenfeld G, Clark D and Studitsky V 2000 Biophys. Chem. 86 231 

[41] Luger K and Richmond T J 1998 Curr. Opin. Cenet. Dev. 8 140. 

[42] Grosbcrg A Yu, Nguyen T T and Shklovskii B I 2002 Rev. Mod. Phys. 74 329 

[43] Park S Y, Bruinsma R F and Gelbart W M 1999 Europhys. Lett. 46 454 

[44] Mateescu E M, Jeppesen C and Pincus P 1999 Europhys. Lett. 46 493 

[45] Gurovitch E and Sens P 1999 Phys. Rev. Lett. 82 339 

[46] Netz R R and Joanny J-F 1999 Macromolecules 32 9026 

[47] Kunze K-K and Netz R R 2000 Phys. Rev. Lett. 85 4389 

[48] Nguyen T T and Shklovskii B I 2001 Physica A 293 324 

[49] Nguyen T T and Shklovskii B I 2001 J. Chem. Phys. 114 5905 

[50] Nguyen T T and Shklovskii B I 2001 J. Chem. Phys. 115 7298 

[51] Schiessel H, Bruinsma R F and Gelbart W M 2001 J. Chem. Phys. 115 7245 

[52] Nguyen T T, Grosbcrg A Yu and Shklovskii B I 2001 in: Electrostatic Effects in Soft Matter and 

Biophysics, eds: C Holm, P Kekicheff and R Podgornik (Kluwer: Dordrecht) p. 469 
[53] Gelbart W M 2001 in: Electrostatic Effects in Soft Matter and Biophysics, eds: C Holm, P 

Kekicheff and R Podgornik (Kluwer: Dordrecht) p. 53 
[54] Kunze K-K and Netz R R 2002 Phys. Rev. E 66 011918 
[55] Schiessel H 2002 Macromolecules in press 
[56] Odijk T 1980 Macromolecules 13 1542 

[57] Pincus P A, Sandroff C J and Witten T A 1984 J. Phys. (Paris) 45 725 



CONTENTS 



87 



[58] von Goclcr F and Muthukumar M 1994 J. Chem. Phys. 100 7796 

[59] Wallin T and Linse P 1996 Langmuir 12 305 

[60] Wallin T and Linse P 1996 J. Phys. Chem. 100 17873 

[61] Wallin T and Linse P 1997 J. Phys. Chem. B 101 5506 

[62] Wallin T and Linse P 1998 J. Chem. Phys. 109 5089 

[63] Haronska P, Vilgis T A, Grottenmiiller R and Schmidt M 1998 Macrom. Theory Simul. 7 241 
[64] Jonsson M and Linse P 2001 J. Chem. Phys. 115 3406 

[65] Chodanowski P and Stoll S 2001 Macromolecules 34 2320 

[66] Sakaue T, Yoshikawa K, Yoshimura S H and Takeyasu K 2001 Phys. Rev. Lett. 87 078105 

[67] Messina R, Holm C and Kremer K 2002 Phys. Rev. E 65 041805 

[68] Messina R, Holm C and Kremer K 2002 J. Chem. Phys. 117 2947 

[69] Jonsson M and Linse P 2001 J. Chem. Phys. 115 10975 

[70] Akinchina A and Linse P 2002 Macromolecules 35 5183 

[71] Dzubiella J, Moreira A G and Pincus P A 2003 Macromolecules36 1741 

[72] Yager T D, McMurray C T and van Holde K E 1989 Biochemistry 28 2271 

[73] Marky N L and Manning G S 1991 Biopolymers 31 1543 

[74] Schiessel H, Rudnick J, Bruinsma R and Gelbart W M 2000 Europhys. Lett. 51 237 

[75] Schiessel H, Widom J, Bruinsma R F and Gelbart W M 2001 Phys. Rev. Lett. 86 4414; Phys. 

Rev. Lett. 2002 88 129902 
[76] Kulic I and Schiessel H, Biophys. J., in press 
[77] Kulic I and Schiessel H, preprint 
[78] de Gennes P G 1971 J. Chem. Phys. 55 572 

[79] de Gcnncs P-G 1979 Scaling Concepts in Polymer Physics (Cornell University Press: Ithaca) 
[80] Doi M and Edwards S F 1986 The Theory of Polymer Dynamcis (Clerendon: Oxford) 
[81] Cui Y and Bustamante C 2000 Proc. Natl. Acad. Sci. USA 97 127 

[82] Bennink M L, Leuba S H, Leno G H, Zlatanova J, de Grooth B G and Greve J 2001 Nat. Struct. 
Biol. 8 606 

[83] Brower-Toland B D, Smith C L, Yeh R C, Lis J T, Peterson C L and Wang M D 2002 Proc. Natl. 
Acad. Sci. USA 99 1960 

[84] Ben-Haim E, Lcsnc A and Victor J-M 2001 Phys. Rev. E 64 0519219 
[85] Victor J-M, Ben-Haim E and Lesne A 2002 Phys. Rev. E 66 060901 (R) 
[86] Katritch V, Bustamante C and Olson W K 2000 J. Mol. Biol. 295 29 
[87] Wedemann G and Langowski J 2002 Biophys. J. 82 2847 

[88] Mangenot S, Lcforestior A, Vachcttc P, Dnrand D and Livolant F 2002 Biophys. J. 82 345 
[89] Mangenot S, Raspaud E, Tribet C, Belloni L and Livolant F 2002 Eur. Phys. J. E 7 221 
[90] Marko J F and Siggia E D 1997 Biophys. J. 73 2 
[91] Schiessel H 2002 Europhys. Lett. 58 140 

[92] Houchmandzadeh B, Marko JF, Chatenay D and Libchaber A 1997 J. Cell. Biol. 139 1 

[93] Poirier M, Eroglu S, Chatenay D and Marko J F 2000 Mol. Biol. Cell 11 269 

[94] Poirier M G, Ncmani A, Gupta P, Eroglu S, Marko J F 2001 Phys. Rev. Lett. 86 360 

[95] Poirier M G, Monhait T, Marko J F 2002 J. Cell. Biochem. 85 422 

[96] Marko J F and Siggia E D 1997 Mol. Biol. Cell 8 2217 

[97] Arents G, Burhngame R W, Wang B-C, Love W E and Moudrianakis E N 1991 Proc. Natl. Acad. 
Sci. USA 88 10148 

[98] Anderson J D and Widom J 2000 J. Mol. Biol. 296 979 

[99] Hagerman P J 1988 Annu. Rev. Biophys. Biophys. Chem. 17 265 
[100] Harries R A and Hearst J E 1966 J. Chem. Phys. 44 2595 
[101] Strauss J K and Mahcr L J 1994 Science 266 1829 
[102] Braunlin W H, Strick T J and Record M T 1982 Biopolymers 21 1301 
[103] Mascotti D P and Lohmann T M 1990 Proc. Natl. Acad. Sci. USA 87 3142 
[104] Radler J O, Koltover L, Salditt T and Safinya C R 1999 Science 275 810 



CONTENTS 



88 



[105] Koltover I, Salditt T and Safinya C R 1999 Biophys. J. 77 915 

[106] Record M T, Anderson C F and Lohmann T M 1978 Q. Rev. Biophys. 11 103 

[107] Bruinsma R 1998 Eur. Phys. J. B 4 75 

[108] Harries D, May S, Gelbart W M and Ben-Shaul A 1998 Biophys. J. 75 159 

[109] Sens P and Joanny J-F 2000 Phys. Rev. Lett. 84 4862 

[110] Oosawa F 1971 Poly electrolytes (Decker: New York) 

[111] Manning G S 1978 Q. Rev. Biophys. 11 179 

[112] Rouzina I and Bloomfield V A 1996 J. Phys. Chem. 100 4292 

[113] Alexander S, Chaikin P M, Grant P, Morales G J, Pincus P and Hone D 1984 J. Chem. Phys. 
80 5776 

[114] Golestanian R 1999 Phys. Rev. Lett. 83 2473; Sens P ihib. 2474 

[115] Nctz R R and Joanny J-F 1999 Macromolecules 32 9013 

[116] Schiessel H and Aranda-Espinoza H 2001 Eur. Phys. J. E 5 499 

[117] Nguyen T T, Grosberg A Yu and Shklovskii B I 2000 J. Chem. Phys. 113 1110 

[118] Ganachaud F, Elai'ssari A, Pichot F, Laayoun A and Cros P 1997 Langmuir 13 701 

[119] Sukhonikov G B, Donath E, Davis S, Lichtenfeld H, Caruso F, Popov Vi I and Mohwald H 1998 

Polym. Adv. Technol. 9 759 
[120] Kabanov A V, Sergeyev V G, Pyshkina O A, Zinchenko A A, Zezin A B, Joosten J G H, Brackman 

J, Yoshikawa K 2000 Macromolecules 33 9587 
[121] McQuigg D W, Kaplan J I and Dubin P L 1992 J. Phys. Chem. 96 1973 and references therein 
[122] Prank-Kamentskii M D 1997 Phys. Rep. 288 13 
[123] Schlick T 1995 Curr. Opin. Struct. Biol. 5 245 

[124] Love A E 1944 A Treatise of the Mathematical Theory of Elasticity 4th edition (Dover: New 

York), Chapts. 18 and 19 
[125] Nizette M and Goriely A 1999 J. Math. Phys. 40 2830 
[126] Bcnham C J 1977 Proc. Natl. Acad. Sci. USA 74 2397 
[127] Benham C J 1979 Biopolymers 18 609 
[128] Le Bret M 1979 Biopolymers 18 1709 
[129] Le Bret M 1984 Biopolymers 23 1835 
[130] Shi Y M and Hearst J E 1994 J. Chem. Phys. 101 5186 
[131] Fain B and Rudnick J 1997 Phys. Rev. E 55 7364 
[132] Swigdon D, Coleman B D and Tobias I 1998 Biophys. J. 74 2515 
[133] Fain B and Rudnick J 1999 Phys. Rev. E 60 7239 
[134] Coleman B D, Swigdon D and Tobias 1 2000 Phys. Rev. E 61 747 
[135] Zandi R and Rudnick J 2001 Phys. Rev. E 64 051918 
[136] Yamakawa H and Stockmayer W 1972 J. Chem. Phys. 57 2843 
[137] Beard P 1978 Cell 15 955 

[138] Spadafora C, Oudet P and Chambon P 1979 Eur. J. Biochem. 100 225 
[139] Watkins J F and Smerdon M J 1985 Biochemistry 24 7279 
[140] Simpson R T and Stafford D W 1983 Proc. Natl. Acad. Sci. USA 50 51 
[141] Wu H-M and Crothers D M 1984 Nature 308 509 

[142] Loomans D, Schiessel H and Blumen A 1997 J. Chem. Phys. 107 2636 
[143] Schiessel H, Sokolov I M and Blumen A 1997 Phys. Rev. E 56 R2390 

[144] Schwinefus J J and Bloomfield V A 2000 Biopolymers 54 572 
[145] Ura K, Nightingale K and Wolffe A P 1996 EMBO J 15 4959 

[146] Ura K, Kurumizaka H, Dimitrov S, Almouzni G and Wolffe A P 1997 EMBO J 16 2096 
[147] Flaus A, Luger K, Tan S and Richmond T J 1996 Proc. Natl. Acad. Sci. USA 93 1370 
[148] Flaus A and Richmond T J 1998 J. Mol. Biol. 275 427 

[149] Hamiche A, Kang J-G, Dennis C, Xiao H and Wu C 2001 Proc. Natl. Acad. Sci. USA 98 14316 
[150] van Holde K E and Yager T D 1985 in: Structure and Function of the Genetic Apparatus Nicolini 
C and Ts'o POP (Eds.), vol. 98, p. 35, NATO ASI series A 



CONTENTS 



89 



151] Sukhishvili S A, Chen Y, Miiller J D, Gratton E, Schweizer K S and Granick S 2000 Nature 
(London) 406 146 

152] Sukhishvili S A, Chen Y, Miiller J D, Gratton E, Schweizer K S and Granick S 2002 

Macromolecules 35 1776 
153] Abramowitz M and Stegun I 1972 Handbook of Mathematical Functions (Dover: New York) 
154] Bouchaud J P and Georges A 1990 Phys. Rep. 195 127 
155] Klafter J, Shlesinger M F and Zumofen G 1993 Chem. Phys. 177 821 
156] Sokolov I M, Mai J and Blumen A 1997 Phys. Rev. Lett. 79 857 
157] Mandelbrot B B 1982 The Fractal Geometry of Nature (Freeman: San Francisco) 
158] de Cloizeaux J and Jannick G 1990 Polymers in Solution (Clarendon Press: Oxford) 
159] Frenkel J and Kontorova T 1938 Zh. Eksp. Teor. Fiz. 8 1340; 1939 J. Phys. (Moscow) 1 137 
160] Floria L M and Mazo J J 1996 Adv. Phys. 45 505 

161] Braun O M, Dauxois T, Paliy M V and Peyrard M 1997 Phys. Rev. Lett. 78 1295 

162] Strunz T and Elmer F-J 1998 Phys. Rev. E 58 1601 

163] Pietroncro L, Schneider W R and Strassler S 1981 Phys. Rev. B 24 2187 

164] Aubry S 1983 J. Phys. (France) 44 147 

165] Watanabe S, van der Zant H S J, Strogatz S H and Orlando T P 1996 Physica D 97 429 
166] Uhler W and Schilling R 1988 Phys. Rev. B 37 5787 

167] Schicsscl H, Oshanin G, Cazabat A M and Moreau M 2002 Phys. Rev. E 66 056130 
168] Marko J F 1998 Phys. Rev. E 57 2134 
169] Marko J F 1997 Europhys. Lett. 38 183 

170] Kamien R D, Lubensky T V, Nelson P and O'Hern C S 1997 Europhys. Lett. 38 237 
171] Levinthal C and Crane H 1956 Proc. Natl. Acad. Sci. USA 42 436 
172] Nelson P Proc. Natl. Acad. Sci. USA 96 14342 

173] van Kampen N G 1992 Stochastic Processes in Physics and Chemistry (North-Holland: 

Amsterdam) 

174] Mattel S, Sampaolese B, De Santis P and Savino M 2002 Biophys. Chem. 97 173 

175] Anselmi C, Bocchinfuso G, De Santis P, Savino M and Scipioni A 2000 Biophys. J. 79 601 

176] Lorch Y, Zhang M and Kornberg R D 1999 Cell 96 389 

177] Bazzett- Jones D P, Cote J, Landel C C, Peterson C L and Workman J L 1999 Mol. Cell. Biol. 
19 1470 

178] Langst G and Becker P B 2001 Molecular Cell 8 1085 

179] Rees A R, Keller R W, Vesenka J P, Yang G and Bustamanto C 1993 Science 260 1646 
180] Schulz A, Muckc N, Langowski J and Rippc K 1998 ,/. Mol. Biol. 283 821 
181] Rivetti C, Guthold M and Bustamanto C 1999 EMBO J 18 4464 
182] Kulic I, private communication 

183] Woodcock C L and Dimitrov S 2001 Curr. Opin. Genet. Dev. 11 130 
184] Zlatanova J, Leuba S H and van Holde K 1998 Biophys. J. 74 2554 
185] Maret G 1976 PhD thesis Konstanz University 

186] Marct G, v. Schickfus M, Mayer A and Dransfeld K 1975 Phys. Rev. Lett. 35 397 

187] Marct G and Weill G 1983 Biopolym,ers 22 2727 

188] Martino J A, Katritch V and Olson W K 1999 Struct. Fold. Des. 7 1009 

189] Widom J 1992 Proc. Natl. Acad. Sci. USA 89 1095 

190] Plcwa J S and Wittcn T A 2000 J. Chem,. Phys. 112 10042 

191] Stryer L 1995 Biochemistry 4th Ed. (Freeman: New York) p 420 

192] Mergell B, Schiessel H and Everaers R, in preparation 

193] Daban J-R 2000 Biochemistry 39 3861 

194] Janninck G, Duplanticr B and Sikorav J L 1996 Biophys. J. 71 451 

195] Cluzel P, Lebrun A, Heller C, Lavery R,Viovy J L, Chatenay D and Caron F 1996 Science 271 
792 

[196] Smith S B, Cui Y and Bustamante C 1996 Science 271 795 



CONTENTS 



90 



[197] Mergell B, Ejtchadi M J and Everaers R 2002 Phys. R.ev. E 66 011903 

[198] Nyrkova I A, Semenov A N, Joanny J-F and Khokhlov A R 1996 J. Phys. II France 6 1411 

[199] Livolant F and Leforestier A 2000 Biophys. J. 78 2716 

[200] Fraden S and Kamien R D 2000 Biophys. J. 78 2189 

[201] Ehrlich L, Miinkel C, Chirico G and Langowski J 1997 CABIOS 13 271 

[202] Beard D A and Schlick T 2001 Structure 9 105 

[203] Miinkel C and Langowski J 1998 Phys. Rev. E 57 5888 

[204] Ladoux B, Quivy J-P, Doyle P, du Roure O, Almouzni G and Viovy J-L 2000 Proc. Natl. Acad. 

Sci. USA 97 14251 

[205] Bennink M L, Pope L H, Leuba S H, de Grooth B G and Greve J 2001 Single Mol. 2 91 
[206] Rief M, Gautel M, Oesterhelt F, Fernandez J and Gaub H 1997 Science 276 1109 
[207] Obcrliauser A, Hansma P, Carrion- Vazquez M and Fernandez J 2001 Proc. Natl. Acad. Sci. USA 
98 468 

[208] Kulic I and Schiessel H, in preparation 
[209] Barrat J-L and Joanny J-F 1996 Adv. Chem. Phys. 94 1 
[210] Toth K, Brun N and Langowski J 2001 Biochemistry 40 6921 
[211] Grunstein M 1997 Nature 389 349 

[212] Hansen J C, Tse C and Wolffe A P 1998 Biochemistry 37 17637 
[213] Strahl B D and AUis D 2000 Nature 403 41 



